• <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é)任編輯 周翠鳴)

    日本精品一区二区三区蜜桃| 波多野结衣一区麻豆| 久久久久久亚洲精品国产蜜桃av| 99久久99久久久精品蜜桃| 午夜两性在线视频| 欧美精品高潮呻吟av久久| 亚洲av电影在线进入| 国产熟女午夜一区二区三区| 精品一区二区三区av网在线观看 | 水蜜桃什么品种好| 三上悠亚av全集在线观看| 黄片小视频在线播放| 亚洲av日韩在线播放| 十分钟在线观看高清视频www| 久久99热这里只频精品6学生| 少妇 在线观看| 波多野结衣av一区二区av| 日韩,欧美,国产一区二区三区| 18在线观看网站| 法律面前人人平等表现在哪些方面 | 亚洲第一欧美日韩一区二区三区 | 一二三四在线观看免费中文在| 巨乳人妻的诱惑在线观看| 国内毛片毛片毛片毛片毛片| 成年女人毛片免费观看观看9 | 爱豆传媒免费全集在线观看| 超碰成人久久| 女人高潮潮喷娇喘18禁视频| 色94色欧美一区二区| av福利片在线| 午夜福利视频精品| 999久久久国产精品视频| 一二三四社区在线视频社区8| 亚洲人成电影观看| 国产高清videossex| 97人妻天天添夜夜摸| 男女无遮挡免费网站观看| 国产免费福利视频在线观看| 法律面前人人平等表现在哪些方面 | 亚洲九九香蕉| 菩萨蛮人人尽说江南好唐韦庄| 黄色a级毛片大全视频| 少妇 在线观看| 在线永久观看黄色视频| 国产精品 国内视频| 色综合欧美亚洲国产小说| 我要看黄色一级片免费的| 在线观看人妻少妇| 美女脱内裤让男人舔精品视频| 日韩一卡2卡3卡4卡2021年| 久久 成人 亚洲| 亚洲五月色婷婷综合| 精品一区二区三卡| 美女脱内裤让男人舔精品视频| 99热网站在线观看| 欧美日韩国产mv在线观看视频| 嫩草影视91久久| 亚洲中文av在线| 一区二区av电影网| 久久青草综合色| 老司机亚洲免费影院| 精品亚洲成a人片在线观看| 日韩视频一区二区在线观看| 一本—道久久a久久精品蜜桃钙片| 午夜福利影视在线免费观看| 一区二区日韩欧美中文字幕| 欧美精品一区二区免费开放| 亚洲色图综合在线观看| 久久久久久久大尺度免费视频| www.999成人在线观看| 中文欧美无线码| 久久久久精品国产欧美久久久 | 亚洲七黄色美女视频| 老司机午夜福利在线观看视频 | 在线观看www视频免费| 亚洲第一青青草原| 一区二区三区精品91| 久久影院123| 又紧又爽又黄一区二区| kizo精华| 亚洲美女黄色视频免费看| 国产成人av激情在线播放| 欧美日韩黄片免| 捣出白浆h1v1| 久久久久久亚洲精品国产蜜桃av| 国产精品1区2区在线观看. | 狠狠精品人妻久久久久久综合| 美女高潮到喷水免费观看| 99久久国产精品久久久| netflix在线观看网站| 男女午夜视频在线观看| 男女之事视频高清在线观看| 高清在线国产一区| 亚洲熟女毛片儿| 人人妻人人爽人人添夜夜欢视频| 性少妇av在线| 国产精品99久久99久久久不卡| 日韩欧美一区二区三区在线观看 | 一本一本久久a久久精品综合妖精| 国产av又大| 啦啦啦中文免费视频观看日本| 亚洲专区中文字幕在线| 日韩视频一区二区在线观看| 欧美另类一区| 成人av一区二区三区在线看 | 老熟妇乱子伦视频在线观看 | 中文精品一卡2卡3卡4更新| 美女国产高潮福利片在线看| 他把我摸到了高潮在线观看 | 亚洲 欧美一区二区三区| 亚洲全国av大片| 精品久久久久久久毛片微露脸 | 亚洲 国产 在线| 黄色视频,在线免费观看| 各种免费的搞黄视频| 色老头精品视频在线观看| 精品人妻在线不人妻| 18在线观看网站| 精品国产乱子伦一区二区三区 | 国产亚洲欧美精品永久| 又紧又爽又黄一区二区| 捣出白浆h1v1| 99久久精品国产亚洲精品| 免费在线观看视频国产中文字幕亚洲 | 日本av免费视频播放| 久久国产亚洲av麻豆专区| 精品亚洲乱码少妇综合久久| a级毛片在线看网站| 天堂8中文在线网| 午夜视频精品福利| 免费高清在线观看日韩| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品免费视频内射| 日本黄色日本黄色录像| 性少妇av在线| 午夜两性在线视频| 亚洲中文av在线| 亚洲成国产人片在线观看| 国产精品1区2区在线观看. | videos熟女内射| 国产成人精品久久二区二区91| 正在播放国产对白刺激| 纯流量卡能插随身wifi吗| 国产精品自产拍在线观看55亚洲 | 国产一区二区在线观看av| 欧美 亚洲 国产 日韩一| 欧美黑人精品巨大| 欧美精品一区二区免费开放| 超碰成人久久| 日韩,欧美,国产一区二区三区| e午夜精品久久久久久久| 老司机亚洲免费影院| 久久中文看片网| 三级毛片av免费| 色婷婷久久久亚洲欧美| 国产一卡二卡三卡精品| 欧美日韩一级在线毛片| 免费在线观看视频国产中文字幕亚洲 | 操美女的视频在线观看| 亚洲一区中文字幕在线| 久久av网站| 亚洲国产av影院在线观看| 老司机亚洲免费影院| 1024香蕉在线观看| 免费观看a级毛片全部| 9191精品国产免费久久| 日韩精品免费视频一区二区三区| 国产欧美日韩综合在线一区二区| 老司机深夜福利视频在线观看 | 精品少妇内射三级| 亚洲国产看品久久| 香蕉国产在线看| 亚洲国产欧美网| 巨乳人妻的诱惑在线观看| 久久精品亚洲av国产电影网| 不卡av一区二区三区| 日韩欧美一区二区三区在线观看 | 美女高潮到喷水免费观看| 成人国产av品久久久| 人妻久久中文字幕网| 亚洲欧美激情在线| 久久女婷五月综合色啪小说| 亚洲国产中文字幕在线视频| 每晚都被弄得嗷嗷叫到高潮| av网站免费在线观看视频| 黄色a级毛片大全视频| 国产日韩欧美亚洲二区| 国产黄色免费在线视频| 免费一级毛片在线播放高清视频 | 高清黄色对白视频在线免费看| 男女国产视频网站| 在线亚洲精品国产二区图片欧美| 人成视频在线观看免费观看| 精品少妇内射三级| 考比视频在线观看| 亚洲国产av影院在线观看| 下体分泌物呈黄色| 91精品三级在线观看| 国产成人免费观看mmmm| 久久人人爽人人片av| 国产亚洲欧美在线一区二区| 男女国产视频网站| 嫩草影视91久久| 最近最新免费中文字幕在线| 国产成人a∨麻豆精品| 日韩视频在线欧美| 麻豆国产av国片精品| 搡老熟女国产l中国老女人| 人妻 亚洲 视频| 操出白浆在线播放| 97人妻天天添夜夜摸| 在线天堂中文资源库| 久久久久国产一级毛片高清牌| 一边摸一边做爽爽视频免费| 国产精品.久久久| 午夜激情久久久久久久| 午夜视频精品福利| 老鸭窝网址在线观看| 亚洲欧美日韩高清在线视频 | 成年美女黄网站色视频大全免费| 一区二区三区四区激情视频| av视频免费观看在线观看| 国产国语露脸激情在线看| 久久久久网色| 90打野战视频偷拍视频| 青春草亚洲视频在线观看| 亚洲精品国产av蜜桃| 69av精品久久久久久 | 丝袜在线中文字幕| 男女之事视频高清在线观看| 在线观看免费高清a一片| 久久狼人影院| 国产成人免费观看mmmm| 美女视频免费永久观看网站| 这个男人来自地球电影免费观看| 99国产精品一区二区三区| 黄片小视频在线播放| 国产97色在线日韩免费| 性色av一级| 久久精品国产亚洲av高清一级| xxxhd国产人妻xxx| 男人添女人高潮全过程视频| 别揉我奶头~嗯~啊~动态视频 | tocl精华| av免费在线观看网站| 亚洲av日韩精品久久久久久密| 免费久久久久久久精品成人欧美视频| 女人被躁到高潮嗷嗷叫费观| 最新的欧美精品一区二区| 一本久久精品| 91精品三级在线观看| 午夜激情久久久久久久| 欧美+亚洲+日韩+国产| 国产精品 国内视频| 午夜日韩欧美国产| 亚洲视频免费观看视频| av片东京热男人的天堂| 韩国高清视频一区二区三区| 啪啪无遮挡十八禁网站| 咕卡用的链子| 国产精品九九99| 欧美国产精品一级二级三级| 欧美精品高潮呻吟av久久| 99精品久久久久人妻精品| 老司机在亚洲福利影院| 99精国产麻豆久久婷婷| 色视频在线一区二区三区| 日本欧美视频一区| 久久精品久久久久久噜噜老黄| 国产精品免费视频内射| 香蕉丝袜av| 淫妇啪啪啪对白视频 | 免费不卡黄色视频| 超碰97精品在线观看| 在线看a的网站| 丁香六月欧美| 99精品欧美一区二区三区四区| 成人国产av品久久久| 精品人妻在线不人妻| 日韩一卡2卡3卡4卡2021年| 激情视频va一区二区三区| 狠狠婷婷综合久久久久久88av| 中文字幕av电影在线播放| 日韩欧美一区二区三区在线观看 | 香蕉丝袜av| 啦啦啦中文免费视频观看日本| 国产精品香港三级国产av潘金莲| 免费不卡黄色视频| 999久久久国产精品视频| 亚洲精品久久久久久婷婷小说| 美国免费a级毛片| 十八禁网站免费在线| 亚洲欧洲日产国产| 男女午夜视频在线观看| 一进一出抽搐动态| 国产成人欧美| 国产精品偷伦视频观看了| 欧美 亚洲 国产 日韩一| 国产精品九九99| 91大片在线观看| 美女中出高潮动态图| 女人精品久久久久毛片| 自线自在国产av| 免费高清在线观看视频在线观看| 男人爽女人下面视频在线观看| 精品久久久久久久毛片微露脸 | 另类亚洲欧美激情| 欧美一级毛片孕妇| 精品人妻一区二区三区麻豆| 人人妻人人澡人人看| 热re99久久精品国产66热6| 欧美精品人与动牲交sv欧美| 9热在线视频观看99| 大陆偷拍与自拍| 日本av手机在线免费观看| 考比视频在线观看| 人妻久久中文字幕网| 极品少妇高潮喷水抽搐| 国产欧美日韩综合在线一区二区| 亚洲精品国产av蜜桃| 香蕉国产在线看| 国产亚洲av高清不卡| 亚洲国产精品一区二区三区在线| 中文字幕精品免费在线观看视频| 最近中文字幕2019免费版| 国产成人免费无遮挡视频| 1024香蕉在线观看| 夫妻午夜视频| 国产精品香港三级国产av潘金莲| 亚洲精品乱久久久久久| 精品人妻1区二区| 日本一区二区免费在线视频| 欧美日韩亚洲高清精品| 一区二区三区精品91| 777久久人妻少妇嫩草av网站| 91av网站免费观看| 91九色精品人成在线观看| 亚洲美女黄色视频免费看| 亚洲欧美日韩另类电影网站| 欧美精品人与动牲交sv欧美| 午夜激情av网站| 最近最新中文字幕大全免费视频| 热99国产精品久久久久久7| 黄色视频,在线免费观看| 一级毛片精品| 亚洲精品国产精品久久久不卡| 淫妇啪啪啪对白视频 | 91麻豆av在线| 国产精品麻豆人妻色哟哟久久| 性少妇av在线| 新久久久久国产一级毛片| 极品人妻少妇av视频| 自线自在国产av| 别揉我奶头~嗯~啊~动态视频 | 女人精品久久久久毛片| 精品人妻一区二区三区麻豆| 成人国语在线视频| 免费黄频网站在线观看国产| 精品久久蜜臀av无| 波多野结衣av一区二区av| 最新在线观看一区二区三区| 国产亚洲精品第一综合不卡| 日本wwww免费看| 久久久久久人人人人人| 69精品国产乱码久久久| 狠狠精品人妻久久久久久综合| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜精品国产一区二区电影| 一区在线观看完整版| 美女视频免费永久观看网站| 欧美日韩福利视频一区二区| www.自偷自拍.com| 精品久久久久久电影网| 悠悠久久av| 国产精品成人在线| 国产亚洲一区二区精品| 在线观看免费日韩欧美大片| 国产亚洲av高清不卡| 国产精品久久久久久人妻精品电影 | 国产精品亚洲av一区麻豆| 飞空精品影院首页| 大香蕉久久成人网| 秋霞在线观看毛片| 麻豆av在线久日| 欧美激情久久久久久爽电影 | 在线观看www视频免费| 久久免费观看电影| 免费在线观看视频国产中文字幕亚洲 | 狠狠狠狠99中文字幕| 啦啦啦中文免费视频观看日本| 另类亚洲欧美激情| √禁漫天堂资源中文www| 免费黄频网站在线观看国产| 免费不卡黄色视频| 精品亚洲乱码少妇综合久久| 波多野结衣一区麻豆| 性色av一级| 免费高清在线观看视频在线观看| 97精品久久久久久久久久精品| 亚洲av欧美aⅴ国产| 王馨瑶露胸无遮挡在线观看| 国产精品一区二区在线观看99| 亚洲欧美日韩高清在线视频 | 亚洲精品国产精品久久久不卡| 777米奇影视久久| 无遮挡黄片免费观看| 亚洲国产av新网站| 伊人久久大香线蕉亚洲五| 制服诱惑二区| 精品国产超薄肉色丝袜足j| 麻豆av在线久日| av免费在线观看网站| 男女免费视频国产| 国产精品熟女久久久久浪| 日韩视频一区二区在线观看| 亚洲欧美一区二区三区久久| 亚洲九九香蕉| 国产在视频线精品| 51午夜福利影视在线观看| 热99国产精品久久久久久7| 啪啪无遮挡十八禁网站| 秋霞在线观看毛片| 大香蕉久久成人网| 老汉色∧v一级毛片| 国产极品粉嫩免费观看在线| 丁香六月天网| 亚洲国产中文字幕在线视频| 欧美黄色片欧美黄色片| 国产av国产精品国产| 19禁男女啪啪无遮挡网站| 美女福利国产在线| 俄罗斯特黄特色一大片| 岛国在线观看网站| 久久久久国内视频| 黄色 视频免费看| 国产成人av教育| av在线app专区| 一本色道久久久久久精品综合| 成年动漫av网址| 搡老乐熟女国产| 99re6热这里在线精品视频| 亚洲精品av麻豆狂野| 日本精品一区二区三区蜜桃| 9191精品国产免费久久| 高清黄色对白视频在线免费看| 9热在线视频观看99| 精品视频人人做人人爽| 精品久久蜜臀av无| 久久久精品免费免费高清| 欧美精品av麻豆av| 亚洲国产日韩一区二区| 国产成+人综合+亚洲专区| 最近最新中文字幕大全免费视频| 国产精品国产三级国产专区5o| 久9热在线精品视频| 乱人伦中国视频| 一本—道久久a久久精品蜜桃钙片| 精品视频人人做人人爽| 青青草视频在线视频观看| 国产精品九九99| 女人爽到高潮嗷嗷叫在线视频| 51午夜福利影视在线观看| 国产精品免费视频内射| 中文字幕人妻熟女乱码| 无遮挡黄片免费观看| 日本猛色少妇xxxxx猛交久久| 黄片小视频在线播放| 国产免费视频播放在线视频| 91av网站免费观看| 狂野欧美激情性bbbbbb| 久9热在线精品视频| 新久久久久国产一级毛片| 亚洲avbb在线观看| 国产91精品成人一区二区三区 | 亚洲国产成人一精品久久久| 日韩大片免费观看网站| 久久毛片免费看一区二区三区| 精品国内亚洲2022精品成人 | 国产一区有黄有色的免费视频| 国产亚洲欧美精品永久| 亚洲男人天堂网一区| 久久 成人 亚洲| 免费女性裸体啪啪无遮挡网站| 精品少妇一区二区三区视频日本电影| 欧美人与性动交α欧美精品济南到| 成年av动漫网址| 大码成人一级视频| 欧美黄色片欧美黄色片| 王馨瑶露胸无遮挡在线观看| 天天操日日干夜夜撸| 啦啦啦啦在线视频资源| 久久久久久亚洲精品国产蜜桃av| 国产av一区二区精品久久| 91老司机精品| 成在线人永久免费视频| 国产精品一区二区精品视频观看| 69精品国产乱码久久久| 波多野结衣av一区二区av| 欧美精品av麻豆av| 女人久久www免费人成看片| 欧美 亚洲 国产 日韩一| 亚洲激情五月婷婷啪啪| 欧美+亚洲+日韩+国产| 亚洲国产欧美在线一区| 亚洲成人国产一区在线观看| 一区二区三区四区激情视频| 欧美xxⅹ黑人| 久久久精品区二区三区| 少妇人妻久久综合中文| 日韩欧美一区视频在线观看| 91精品国产国语对白视频| 少妇裸体淫交视频免费看高清 | 国产男女超爽视频在线观看| 亚洲欧美精品综合一区二区三区| 视频区图区小说| 亚洲精品国产精品久久久不卡| 在线十欧美十亚洲十日本专区| 国产色视频综合| 亚洲成人免费av在线播放| 高清在线国产一区| 亚洲av电影在线进入| 丝袜脚勾引网站| 嫁个100分男人电影在线观看| 亚洲专区中文字幕在线| 午夜久久久在线观看| 十八禁网站免费在线| 一本久久精品| 国产黄频视频在线观看| 亚洲成人免费av在线播放| 制服人妻中文乱码| 久久精品人人爽人人爽视色| 国产精品1区2区在线观看. | 久久中文字幕一级| 一级黄色大片毛片| 午夜老司机福利片| 亚洲精品国产区一区二| 国产成人一区二区三区免费视频网站| 国产主播在线观看一区二区| 精品少妇内射三级| 亚洲精华国产精华精| 成年av动漫网址| svipshipincom国产片| 日韩三级视频一区二区三区| 美国免费a级毛片| 免费人妻精品一区二区三区视频| 亚洲伊人色综图| 婷婷丁香在线五月| 99国产精品免费福利视频| 超碰成人久久| 两性夫妻黄色片| 免费黄频网站在线观看国产| 人人澡人人妻人| 久久女婷五月综合色啪小说| 国产精品99久久99久久久不卡| 精品熟女少妇八av免费久了| 美女中出高潮动态图| 国产成人免费无遮挡视频| 欧美在线一区亚洲| 在线亚洲精品国产二区图片欧美| √禁漫天堂资源中文www| 亚洲中文av在线| 欧美日韩精品网址| 午夜老司机福利片| 国产一卡二卡三卡精品| 亚洲欧美一区二区三区久久| 日韩大码丰满熟妇| 久久久国产欧美日韩av| 人人澡人人妻人| 天堂中文最新版在线下载| 国产一级毛片在线| 国产精品免费大片| 国产三级黄色录像| 欧美精品亚洲一区二区| 日本vs欧美在线观看视频| 考比视频在线观看| 9191精品国产免费久久| 中文欧美无线码| 国产又爽黄色视频| 国产av国产精品国产| 一级片免费观看大全| 午夜免费成人在线视频| 搡老乐熟女国产| 丝袜美腿诱惑在线| 丰满饥渴人妻一区二区三| 午夜日韩欧美国产| 12—13女人毛片做爰片一| 色婷婷av一区二区三区视频| kizo精华| 中文字幕精品免费在线观看视频| 一区在线观看完整版| 日韩欧美一区视频在线观看| 精品国产国语对白av| 青草久久国产| 青春草视频在线免费观看| 男男h啪啪无遮挡| 又大又爽又粗| 少妇精品久久久久久久| 丁香六月欧美| 欧美日韩成人在线一区二区| 一个人免费看片子| 国精品久久久久久国模美| 麻豆av在线久日| av在线app专区| 精品一区二区三区av网在线观看 | 在线看a的网站| 丝瓜视频免费看黄片| 欧美久久黑人一区二区| 精品免费久久久久久久清纯 | 激情视频va一区二区三区| 国产男人的电影天堂91| 90打野战视频偷拍视频| 欧美日韩黄片免| 夫妻午夜视频| 少妇裸体淫交视频免费看高清 | 国产黄色免费在线视频|