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

    橋山林區(qū)麻櫟種群不同發(fā)育階段空間格局及關(guān)聯(lián)性

    2017-06-22 13:47:48薛文艷張文輝于世川
    生態(tài)學(xué)報(bào) 2017年10期
    關(guān)鍵詞:麻櫟發(fā)育階段樣地

    薛文艷, 楊 斌, 張文輝, 于世川

    西北農(nóng)林科技大學(xué)陜西省林業(yè)綜合實(shí)驗(yàn)室, 楊凌 712100

    ?

    橋山林區(qū)麻櫟種群不同發(fā)育階段空間格局及關(guān)聯(lián)性

    薛文艷, 楊 斌, 張文輝*, 于世川

    西北農(nóng)林科技大學(xué)陜西省林業(yè)綜合實(shí)驗(yàn)室, 楊凌 712100

    為探究麻櫟種群不同發(fā)育階段的空間分布特征及關(guān)聯(lián)性,以黃土高原南部橋山林區(qū)麻櫟天然次生林為研究對(duì)象,用L(r)函數(shù)對(duì)麻櫟種群不同發(fā)育階段空間分布特征及關(guān)聯(lián)性進(jìn)行分析。結(jié)果表明:(1)樣地內(nèi)麻櫟種群各發(fā)育階段數(shù)量結(jié)構(gòu)均為不規(guī)則倒“J”型,種群結(jié)構(gòu)穩(wěn)定,能實(shí)現(xiàn)持續(xù)更新。(2)麻櫟種群總體的空間分布特征受生境異質(zhì)性影響較大,純林樣地中,種群在0—10 m尺度上表現(xiàn)為高強(qiáng)度聚集分布,5—35 m尺度表現(xiàn)為低強(qiáng)度聚集;混交林樣地中,0—5 m尺度上表現(xiàn)為高強(qiáng)度聚集分布,5—35 m尺度上則表現(xiàn)為低強(qiáng)度聚集分布。(3)隨著種群的發(fā)育及尺度的增大,麻櫟種群空間聚集強(qiáng)度減弱,純林樣地內(nèi),麻櫟幼苗在0—6 m的小尺度表現(xiàn)為聚集分布,6—12 m尺度上表現(xiàn)為隨機(jī)分布;幼樹在各尺度上均表現(xiàn)為聚集分布;中齡樹與成齡樹在各尺度上均表現(xiàn)為隨機(jī)分布?;旖涣謽拥刂?僅幼苗在0—6 m尺度上表現(xiàn)為聚集分布,其余各階段在各尺度下均趨于隨機(jī)分布。(4)種群不同發(fā)育階段不同尺度空間關(guān)聯(lián)性不同,純林樣地中,幼苗與幼樹在1—2 m小尺度上呈負(fù)相關(guān),幼樹與中齡樹在8 m尺度呈正相關(guān),其余各生長(zhǎng)階段之間在任意尺度上均無(wú)關(guān)聯(lián)性;混交林樣地中,幼樹與中齡樹在2—12 m尺度上表現(xiàn)為正相關(guān),其余各生長(zhǎng)階段之間在各尺度上均無(wú)關(guān)聯(lián)。在未來(lái)的森林經(jīng)營(yíng)中,應(yīng)對(duì)不同發(fā)育階段、不同生境的種群進(jìn)行不同撫育措施,促進(jìn)種群良性發(fā)育。

    麻櫟;空間格局;空間關(guān)聯(lián);發(fā)育階段

    種群空間分布格局是指種群內(nèi)個(gè)體在生存空間的分配方式及配置特點(diǎn),是種群與環(huán)境因素相互作用的結(jié)果[1]。分析種群空間分布格局有助于了解種群生長(zhǎng)發(fā)育、生活史規(guī)律,種子傳播、種內(nèi)種間關(guān)系,以及種群與環(huán)境因素的互作效應(yīng)[2- 4]。對(duì)不同發(fā)育階段種群格局類型及其關(guān)聯(lián)性的研究可以更好的揭示種群發(fā)育過(guò)程中對(duì)環(huán)境資源的適應(yīng)和利用狀況[3]。

    研究空間格局可有多種方法,其中以點(diǎn)格局分析應(yīng)用最為廣泛[5- 6]。該方法以個(gè)體在二維空間的坐標(biāo)為基礎(chǔ),構(gòu)成點(diǎn)圖,并以此點(diǎn)圖為基礎(chǔ)進(jìn)行分析。它基于點(diǎn)之間的距離統(tǒng)計(jì),最大限度的利用了點(diǎn)之間的距離信息,可提供更全面的空間格局信息[7]。點(diǎn)格局分析可分析各尺度種群空間分布格局,同時(shí)還可分析種群間關(guān)聯(lián)性,并給出空間分布格局的最大聚集強(qiáng)度及其對(duì)應(yīng)的尺度[8]。由于空間格局具有很強(qiáng)的尺度依賴性,即空間格局的分布特性會(huì)隨尺度的不同而發(fā)生變化,因此,點(diǎn)格局分析的結(jié)果更符合生態(tài)學(xué)規(guī)律。

    麻櫟Quercusacutissima,殼斗科櫟屬,是黃土高原南部橋山地區(qū)闊葉落葉林地帶性植被的建群種,常以純林或松櫟混交林形式構(gòu)成當(dāng)?shù)刈钪匾志坝^,在當(dāng)?shù)睾B(yǎng)水源、保持水土、保護(hù)生物多樣性方面發(fā)揮著重要作用[9]。目前對(duì)櫟屬植物空間格局的研究主要集中于栓皮櫟[10]、遼東櫟[11]、銳齒櫟[12]。袁志良[10]等研究發(fā)現(xiàn)栓皮櫟種群空間格局受空間異質(zhì)性的影響并影響種群的天然更新;常偉等[12]認(rèn)為不同生境銳齒櫟空間分布格局不同,營(yíng)造適宜生境可促進(jìn)種群更新與穩(wěn)定。迄今為止,有關(guān)麻櫟種群空間分布及空間關(guān)聯(lián)鮮有報(bào)道。本研究以橋山林區(qū)麻櫟天然次生純林和油松麻櫟天然次生混交林為研究對(duì)象,用點(diǎn)格局分析法分析了不同發(fā)育階段麻櫟種群空間分布格局及空間關(guān)聯(lián),闡明了種群分布格局與發(fā)育階段的關(guān)系及不同發(fā)育階段的關(guān)聯(lián)性,揭示了種群發(fā)育動(dòng)態(tài)變化機(jī)制,為種群合理經(jīng)營(yíng)及持續(xù)穩(wěn)定更新提供依據(jù)。

    1 研究區(qū)概況

    研究區(qū)選擇在陜北黃土高原延安市橋山林業(yè)局雙龍林場(chǎng),地理位置為108°31′—109°11′ E,35°30′—35°50′ N,海拔800—1700 m,屬暖溫帶氣候區(qū),年平均氣溫9.4 ℃,最高氣溫38 ℃,最低氣溫-23 ℃,無(wú)霜期190—225 d,年均降水量677.4 mm,主要土壤類型為森林褐土和灰褐土。主要喬木樹種有遼東櫟(Quercuswutaishanica)、麻櫟(Quercusacutissima)、油松(Pinustabuliformis)、山楊 (Populusdavidiana)和白樺 (Betulaplatyphylla)等;主要灌木樹種有胡枝子(Lespedezabicolor)、狼牙刺(Sophoradavidii)、刺五加(Eleutherococcussenticosus)、水栒子(Cotoneastermultiflorus)等;主要草本植物有白茅(Imperatacylindrica)、茜草(Rubiacordifolia)等[9]。

    2 研究方法

    2.1 樣地設(shè)置與調(diào)查

    為了準(zhǔn)確揭示麻櫟種群不同發(fā)育階段的空間格局及其關(guān)聯(lián)性,于2015年7—8月經(jīng)過(guò)野外踏查后,根據(jù)點(diǎn)格局研究的取樣要求[7],在雙龍林場(chǎng)南溝,選擇一個(gè)較大的陰坡坡面由一條較大侵蝕溝分割,具有典型特征的林地:林地A為麻櫟純林,林地B為油松-麻櫟混交林。根據(jù)林場(chǎng)保留的林相圖和小班卡記錄,兩塊樣地1976年前作為薪炭林,林場(chǎng)和當(dāng)?shù)剞r(nóng)戶從中經(jīng)常性從中采樵生活用柴,形成比較殘敗的天然次生林(麻櫟純林,油松—麻櫟混交林),1976年實(shí)行封育措施,此后隨著人口遷出,林場(chǎng)職工減少,兩塊林地未經(jīng)受過(guò)較大的外界干擾。在A、B兩塊林地分別設(shè)置面積為50 m×50 m兩塊樣地(A、B樣地)。采用相鄰網(wǎng)格法,分別將兩塊樣地劃分為25個(gè)10 m×10 m的樣方;以樣地左下方為原點(diǎn),對(duì)樣方內(nèi)的麻櫟個(gè)體進(jìn)行坐標(biāo)定位;每木調(diào)查包括植株的高度、胸徑、冠幅;對(duì)于所有高度小于1 m的植株,測(cè)定其基徑、高度、冠幅;記錄樣地海拔、坡向、坡度等環(huán)境因子,并用樣點(diǎn)統(tǒng)計(jì)法調(diào)查林分郁閉度[13],在每塊樣地內(nèi)隨機(jī)設(shè)置50個(gè)樣點(diǎn),觀察樹冠是否遮蓋樣點(diǎn),統(tǒng)計(jì)被遮蓋的樣點(diǎn)數(shù),進(jìn)而得到林分郁閉度,樣地基本情況見(jiàn)表1。

    表1 橋山雙龍林場(chǎng)兩塊樣地基本情況表

    A: 麻櫟純林,Q.acutissimanaturalsecondarypureforest;B: 麻櫟-油松混交林,Q.acutissima-Pinustabuliformisnaturalsecondarymixedforest

    2.2 種群不同發(fā)育階段數(shù)量結(jié)構(gòu)

    不同發(fā)育階段種群的數(shù)量結(jié)構(gòu)在一定的程度上可以反映種群的發(fā)展趨勢(shì)[14]。櫟屬樹木通常樹體高大且木質(zhì)堅(jiān)硬,不易用生長(zhǎng)錐直接測(cè)定其年齡,通常以徑級(jí)結(jié)構(gòu)和高度級(jí)結(jié)構(gòu)代替年齡結(jié)構(gòu)進(jìn)行分析[15- 16]。本次調(diào)查徑級(jí)劃分標(biāo)準(zhǔn)為[17]:Ⅰ級(jí)(高度H ≤ 50 cm)、Ⅱ級(jí)(50 cm < H ≤ 100 cm)、Ⅲ級(jí)(H > 100 cm,DBH < 4 cm)、Ⅳ級(jí)(4 cm ≤DBH < 8 cm)、Ⅴ級(jí)(8 cm ≤ DBH < 12 cm)、Ⅵ級(jí)(12 cm ≤ DBH < 16 cm)、Ⅶ級(jí)(16 cm ≤ DBH < 20 cm)、Ⅷ級(jí)(20 cm ≤ DBH < 24 cm)、Ⅸ級(jí)(24 cm ≤ DBH < 28 cm)和Ⅹ級(jí)(DBH ≥ 28 cm)。依據(jù)實(shí)際調(diào)查結(jié)果并結(jié)合麻櫟生長(zhǎng)特性,將麻櫟種群劃分為4個(gè)生長(zhǎng)階段: s1幼苗(Ⅰ、Ⅱ級(jí))、s2幼樹(Ⅲ、Ⅳ級(jí))、s3中齡樹(Ⅴ—Ⅷ級(jí))和s4成齡樹(Ⅸ—Ⅹ級(jí))。分別統(tǒng)計(jì)各階段株樹,并以其為橫坐標(biāo),以各階段株數(shù)為縱坐標(biāo),做麻櫟種群不同發(fā)育階段數(shù)量結(jié)構(gòu)圖(圖1)。

    2.3 Morishita Index 圖

    Morishita Index(MI) 可用來(lái)表示種群聚集程度,一個(gè)MI表示固定樣方內(nèi)種群的分布格局。將固定樣方連續(xù)分割后得到不同的MI,便可將空間聚集程度與尺度聯(lián)系在一起。具體方法為將每一個(gè)大樣方逐步分割為2×2,3×3,…,n×n個(gè)小樣方,通過(guò)統(tǒng)計(jì)落入樣方內(nèi)的個(gè)體數(shù)計(jì)算一組連續(xù)的MI[18]。

    2.4 點(diǎn)格局分析

    點(diǎn)格局分析法(point pattern analysis)由Ripley最先提出[19],其主要思路為通過(guò)劃定一定面積的研究區(qū)域,以該區(qū)域內(nèi)植物個(gè)體的空間坐標(biāo)構(gòu)成的點(diǎn)圖為基礎(chǔ)數(shù)據(jù)進(jìn)行格局分析。由于其可用于分析任一尺度下種群的空間分布格局,因此在實(shí)際中得到了廣泛的應(yīng)用。Ripley′sK(r)函數(shù)公式為:

    式中,A是研究區(qū)(樣方)的面積;n為個(gè)體數(shù);uij為兩個(gè)點(diǎn)i和j之間的距離;Ir(uij)為指示函數(shù),當(dāng)uijr時(shí),Ir(uij)=0;wij為權(quán)重值,用于邊緣校正。

    通常情況下,為保持方差穩(wěn)定,通常用L(r)代替K(r)函數(shù),其表達(dá)式為:

    當(dāng)L(r)=0時(shí),表示隨機(jī)分布,位于上下包跡線之間;當(dāng)L(r)>0時(shí),為集群分布,位于上下包跡線之上;當(dāng)L(r)<0時(shí),為均勻分布,位于上下包跡線之下。

    不同生長(zhǎng)階段麻櫟個(gè)體間空間關(guān)系可以用雙變量函數(shù)分析,具體計(jì)算公式為:

    同理,可用L12(r)代替K12(r),表達(dá)式為:

    當(dāng)L12(r)=0時(shí),表示兩變量間相互獨(dú)立,位于上下包跡線之間;當(dāng)L12(r)>0時(shí),表示兩變量空間正相關(guān),位于上下包跡線之上;當(dāng)L12(r)<0 時(shí),表示兩變量間呈空間負(fù)相關(guān),位于上下包跡線之下。

    本文通過(guò)Monte-Carlo擬合檢驗(yàn)確定置信區(qū)間,選取空間完全隨機(jī)模型作為零模型[20],重復(fù)擬合99次,得到95%置信區(qū)間,即圖中上、下包跡線之間的范圍。以r為橫坐標(biāo),上下包跡線為縱坐標(biāo)作圖。本文所有數(shù)據(jù)處理、格局分析、統(tǒng)計(jì)圖制作均通過(guò)R軟件spatstat包[21]實(shí)現(xiàn)。

    3 結(jié)果與分析

    3.1 種群不同發(fā)育階段數(shù)量分布圖

    圖1 橋山雙龍林場(chǎng)兩種樣地麻櫟種群不同發(fā)育階段個(gè)體數(shù)量分布圖Fig.1 Number of Quercus acutissimain different growth

    由圖1可知,樣地A中,各階段個(gè)體數(shù)分別占個(gè)體總數(shù)的76.7%、5.4%、15.9%、2.0%。樣地B中,各階段個(gè)體數(shù)分別占個(gè)體總數(shù)的89.4%、2.1%、6.6%、1.9%。兩塊樣地中,麻櫟種群不同發(fā)育階段數(shù)量結(jié)構(gòu)均呈基部寬頂部窄的不規(guī)則倒“J”型。s1階段的麻櫟個(gè)體數(shù)最多,說(shuō)明種群更新?tīng)顟B(tài)良好。盡管幼苗數(shù)量?jī)?chǔ)備豐富,但幼樹數(shù)量不足,樣地A中,由s1階段到s2階段,種群個(gè)體數(shù)減少92.9%,樣地B相同階段種群個(gè)體數(shù)減少97.6%,在此階段競(jìng)爭(zhēng)激烈導(dǎo)致只有少數(shù)幼苗能順利進(jìn)入s2,s2階段成為制約種群發(fā)育的瓶頸期。在瓶頸期之后兩塊樣地麻櫟個(gè)體數(shù)分別增加192%,209%。在s4階段也出現(xiàn)了一個(gè)死亡高峰,但死亡率低于s2階段。

    3.2 麻櫟種群總體分布格局

    由圖2可知,0—35 m范圍內(nèi),樣地A和樣地B麻櫟種群MI始終大于1,兩麻櫟種群均表現(xiàn)為聚集分布。樣地A中,種群MI值在0—10 m尺度內(nèi)迅速減小,10—35 m尺度上MI值穩(wěn)定在1左右,說(shuō)明樣地A中麻櫟種群在10—35 m尺度上表現(xiàn)為低強(qiáng)度聚集分布;樣地B中,0—5 m尺度上種群MI值迅速降低,在5—35 m尺度上表現(xiàn)為低強(qiáng)度聚集分布,種群MI值穩(wěn)定在1左右。

    圖2 橋山雙龍林場(chǎng)兩種樣地麻櫟種群總體在不同尺度下的聚集指數(shù) MIFig.2 Morishita Index of increasing diameter of quadrat

    3.3 麻櫟種群不同發(fā)育階段的空間分布格局

    圖3、圖4表明,s1階段麻櫟種群在0—6 m的小尺度表現(xiàn)為聚集分布,稍大尺度上表現(xiàn)為隨機(jī)分布。不同的是,樣地A中,s2階段種群在0—12 m尺度上均表現(xiàn)為聚集分布;隨著種群發(fā)育至s3,s4階段,分布狀態(tài)逐漸變?yōu)殡S機(jī)分布;樣地B中,s2階段,種群基本表現(xiàn)為隨機(jī)分布,僅在1—2 m尺度上表現(xiàn)為聚集分布;s3階段,種群在3 m尺度表現(xiàn)為均勻分布,12 m尺度上表現(xiàn)為聚集分布;s4階段各尺度上種群均表現(xiàn)為隨機(jī)分布。

    圖3 橋山雙龍林場(chǎng)樣地A麻櫟種群各齡級(jí)空間分布格局Fig.3 Spatial pattern of each stage of Quercus acutissimain plot As1:幼苗,seedling stage;s2:幼樹,sapling stage;s3:中齡樹,mid-mature stages4:成齡樹,mature stage

    圖4 橋山雙龍林場(chǎng)樣地B麻櫟種群各齡級(jí)空間分布格局Fig.4 Spatial pattern of each stage of Quercus acutissimain plot B

    3.4 麻櫟種群不同發(fā)育階段的空間關(guān)聯(lián)性

    圖5、圖6結(jié)果顯示,樣地A中,s1與s2階段的個(gè)體僅在1—2 m小尺度上呈負(fù)相關(guān),s2與s3階段的個(gè)體在8 m尺度呈正相關(guān),其余各生長(zhǎng)階段之間在任意尺度上均無(wú)關(guān)聯(lián)性。樣地B中,s2—s3階段個(gè)體在2—12 m尺度上表現(xiàn)為正相關(guān),其余各生長(zhǎng)階段之間在各尺度上均無(wú)關(guān)聯(lián)。

    圖5 橋山雙龍林場(chǎng)樣地A麻櫟種群各生長(zhǎng)階段空間關(guān)系圖Fig.5 Associations of different growth of Quercus acutissimain plot A

    圖6 橋山雙龍林場(chǎng)樣地B麻櫟種群各生長(zhǎng)階段空間關(guān)系圖Fig.6 Associations of different growth of Quercus acutissimain plot B

    4 結(jié)論與討論

    植物種群的年齡結(jié)構(gòu)是指種群個(gè)體的年齡組配狀況,是植物種群生殖特性與外界環(huán)境相互作用的結(jié)果,它可以反映種群的動(dòng)態(tài)變化和群落的發(fā)展趨勢(shì)[22]。本研究中,麻櫟種群的齡級(jí)結(jié)構(gòu)是幼苗比例最大,中齡樹、成齡樹比例次之,幼樹比例最小。這說(shuō)明該種群屬進(jìn)展性種群,種群能進(jìn)行自我更新,且具有較強(qiáng)的持續(xù)發(fā)育能力。麻櫟幼苗數(shù)量多主要是由其生殖特性決定,麻櫟種群可以產(chǎn)生大量高質(zhì)量種子并萌發(fā)[23]。隨著種群發(fā)育,幼樹階段出現(xiàn)明顯的瓶頸現(xiàn)象,產(chǎn)生這種現(xiàn)象的原因是,隨著樹齡的增大,麻櫟種群不同個(gè)體對(duì)光照、水分、養(yǎng)分等環(huán)境資源需求隨之增大,由于環(huán)境資源有限,種內(nèi)競(jìng)爭(zhēng)加劇,大量的幼苗在競(jìng)爭(zhēng)中會(huì)死亡,導(dǎo)致盡管麻櫟種群幼苗充足, 卻難以轉(zhuǎn)化為中高齡個(gè)體[24]。林下經(jīng)常出現(xiàn)只見(jiàn)幼苗不見(jiàn)幼樹的現(xiàn)象。隨著種群進(jìn)一步發(fā)育,中齡樹對(duì)環(huán)境的適應(yīng)能力增強(qiáng)且中齡樹占據(jù)林冠中上層空間,光照充足,競(jìng)爭(zhēng)力提高,使得種群死亡率降低,中齡樹個(gè)體數(shù)增加[25]。成齡樹死亡率又有所提高,這主要是由于此階段個(gè)體進(jìn)入生理衰老年齡,出現(xiàn)自然死亡??偟膩?lái)說(shuō),橋山地區(qū)麻櫟種群表現(xiàn)出間歇型種群結(jié)構(gòu),但由于麻櫟種群可通過(guò)群落微結(jié)構(gòu)-林窗完成更新[26],因此這種間歇型結(jié)構(gòu)可保持相對(duì)穩(wěn)定,從而促進(jìn)群落的穩(wěn)定。

    麻櫟種群分布格局一方面受生殖特性影響,另一方面也與水分、氣候、光照、土壤等可利用生態(tài)資源有關(guān)[27]。本研究中,麻櫟種群總體上呈聚集分布。這是因?yàn)榫奂植伎梢蕴岣叻N群競(jìng)爭(zhēng)力及對(duì)外界不良環(huán)境的抵抗能力[28]。不同的是,混交林樣地喬木層有油松、白樺、茶條槭,灌木層有胡枝子、狼牙刺等。在白樺枯死木周圍有小面積林窗,麻櫟便在這些林窗聚集。隨取樣尺度增大,林窗效應(yīng)減弱,氣候,土壤水分與養(yǎng)分狀況,群落郁閉度等影響增強(qiáng),導(dǎo)致混交林樣地種群聚集程度隨尺度迅速變化。而純林樣地喬木層以麻櫟為主,在一定尺度內(nèi)生境條件相對(duì)均勻,因此樣地內(nèi)種群聚集程度變化緩慢。

    對(duì)種群在不同發(fā)育階段的空間格局變化的研究可以了解種群幼齡個(gè)體的產(chǎn)生及老齡個(gè)體的死亡等發(fā)育過(guò)程并推斷種群是否遭受人為干擾。本研究選取的兩塊樣地中,幼苗階段個(gè)體呈聚集分布,這主要是由麻櫟的生殖特性及種子散布方式?jīng)Q定。同時(shí),聚集分布使個(gè)體間能相互庇護(hù),進(jìn)而提高了種群的抗性并降低了死亡率。中齡、成齡階段個(gè)體大致呈隨機(jī)分布,這一方面是由于麻櫟幼苗耐蔭而小樹,大樹喜光,另一方面是由于隨著樹齡增大,種內(nèi)及種間競(jìng)爭(zhēng)加劇,自疏及他疏作用加強(qiáng),種群密度降低是造成聚集強(qiáng)度降低的主要因素。幼樹階段混交林樣地內(nèi)個(gè)體呈隨機(jī)分布而純林樣地內(nèi)個(gè)體呈聚集分布,是由于混交林樣地內(nèi)有枯樹形成的林窗,個(gè)體在林窗內(nèi)環(huán)境資源充足,競(jìng)爭(zhēng)力提高,因此趨于隨機(jī)分布。本研究與其他天然群落格局分析結(jié)果類似[2,29]。這是因?yàn)樽匀粻顟B(tài)下,不同發(fā)育階段的個(gè)體占據(jù)不同的空間位置,可以使個(gè)體對(duì)有限環(huán)境資源的利用達(dá)到最大化,同時(shí)還可以提高種群對(duì)環(huán)境的適應(yīng)能力,有利于其生存發(fā)展。由聚集分布到隨機(jī)分布的變化,是種群發(fā)育的一種生存策略和適應(yīng)機(jī)制[30]。

    植物種群的空間分布格局有很強(qiáng)的尺度依賴性。同一物種相同發(fā)育階段在不同尺度上表現(xiàn)出的分布狀態(tài)均有差異[31]。種群在某一尺度上的空間格局是由多種生態(tài)過(guò)程共同形成的。尺度較大時(shí),格局由物種分布、斑塊異質(zhì)性以及環(huán)境因素共同決定;尺度較小時(shí),格局與種內(nèi)種間競(jìng)爭(zhēng)、種子落地及擴(kuò)散方式有關(guān)[32]。本研究中,兩塊樣地中,幼齡個(gè)體在小尺度上呈聚集分布,這主要是由落種及散布方式?jīng)Q定。隨著尺度增大則趨于隨機(jī)分布,這種不同發(fā)育階段空間格局的鑲嵌分布有利于麻櫟種群的穩(wěn)定發(fā)育[33]。

    同一物種不同發(fā)育階段在空間分布上的相互關(guān)聯(lián)性可用雙變量種間關(guān)聯(lián)來(lái)表示,其指示了某一種群不同發(fā)育時(shí)期個(gè)體間的相互關(guān)系,這種關(guān)系既指?jìng)€(gè)體空間分布關(guān)系,也指?jìng)€(gè)體間的功能關(guān)系,可以據(jù)此了解種群的現(xiàn)狀并預(yù)測(cè)種群發(fā)展趨勢(shì)[6]。由于樹木個(gè)體之間的相互影響僅出現(xiàn)在一定范圍內(nèi)[34],因此對(duì)物種的空間分布關(guān)聯(lián)的研究多集中在較小尺度。本研究調(diào)查的兩塊樣地內(nèi),幼苗與其它階段個(gè)體在所有尺度上均為負(fù)關(guān)聯(lián)或無(wú)關(guān)聯(lián),這是由于其他個(gè)體營(yíng)造的郁閉環(huán)境不利于幼苗生長(zhǎng);幼樹與中齡樹均在一定尺度上呈正相關(guān),而與其它階段個(gè)體呈負(fù)相關(guān)或無(wú)關(guān)聯(lián),這種現(xiàn)象可能是由于種內(nèi)競(jìng)爭(zhēng)具有非對(duì)稱性[35],即當(dāng)個(gè)體大小相近時(shí),二者間關(guān)系較協(xié)調(diào);當(dāng)個(gè)體間差異較大時(shí),由于大株個(gè)體對(duì)小株個(gè)體強(qiáng)烈的排斥作用及大株個(gè)體自身較強(qiáng)的競(jìng)爭(zhēng)力,使其周圍小株個(gè)體很少[36]。因此中齡樹與成齡樹在所有尺度上均呈現(xiàn)無(wú)關(guān)聯(lián)或負(fù)相關(guān)。

    鑒于麻櫟種群現(xiàn)狀,建議有針對(duì)性地該地區(qū)麻櫟天然林進(jìn)行撫育。林下灌木和草本會(huì)影響種子萌發(fā)及幼苗存活,導(dǎo)致幼樹階段成為制約種群發(fā)育的瓶頸期。因此應(yīng)加強(qiáng)下層撫育,清除部分雜草及灌木,為幼苗生長(zhǎng)創(chuàng)造適宜環(huán)境。幼樹階段后對(duì)光照要求增加,此時(shí)應(yīng)對(duì)林冠層進(jìn)行透光伐,加快個(gè)體生長(zhǎng)發(fā)育,促進(jìn)種群良性發(fā)育。

    [1] 薛建輝. 森林生態(tài)學(xué). 北京: 中國(guó)林業(yè)出版社, 2006: 75.

    [2] Druckenbrod D L, Shugart H H, Davies I. Spatial pattern and process in forest stands within the Virginia piedmont. Journal of Vegetation Science, 2005, 16(1): 37- 48.

    [3] Nathan R. Long-distance dispersal of plants. Science, 2006, 313(5788): 786- 788.

    [4] 胡爾查, 王曉江, 張文軍, 海龍, 張雷, 張勝利, 徐鵬雁. 烏拉山自然保護(hù)區(qū)白樺種群的年齡結(jié)構(gòu)和點(diǎn)格局分析. 生態(tài)學(xué)報(bào), 2013, 33(9): 2867- 2876.

    [5] Zhang J C, Zhang Y Q, Fan D Q, Qin S G, Jia X, Wu B, Chen D, Gao H, Zhu L F. Effects of sand dune stabilization on the spatial pattern ofArtemisiaordosicapopulation in Mu Us desert, Northwest China. PLoS One, 2015, 10(6): e0129728.

    [6] Shields J M, Jenkins M A, Saunders M R, Zhang H, Jenkins L H, Parks A M. Age distribution and spatial patterning of an invasive shrub in secondary hardwood forests. Forest Science, 2014, 60(5): 830- 840.

    [7] 張金屯. 植物種群空間分布的點(diǎn)格局分析. 植物生態(tài)學(xué)報(bào), 1998, 22(4): 344- 349.

    [8] Kaito C, Dieckmann U, Sasaki A, Takasu F. Beyond pairs: definition and interpretation of third-order structure in spatial point patterns. Journal of Theoretical Biology, 2015, 372(6): 22- 38.

    [9] 曹旭平, 張文輝. 陜西橋山省級(jí)自然保護(hù)區(qū)綜合科學(xué)考察報(bào)告. 楊凌: 西北農(nóng)林科技大學(xué)出版社, 2015: 1- 5.

    [10] 袁志良, 王婷, 朱學(xué)靈, 沙迎迎, 葉永忠. 寶天曼落葉闊葉林樣地栓皮櫟種群空間格局. 生物多樣性, 2011, 19(2): 224- 231.

    [11] 王道亮, 李衛(wèi)忠, 曹鑄,王慶, 王廣儒. 黃龍山天然次生林遼東櫟種群空間格局. 生態(tài)學(xué)報(bào), 2016, 36(9): 1- 10.

    [12] 常偉, 黨坤良, 武朋輝, 李明雨. 秦嶺林區(qū)銳齒櫟次生林種群空間分布格局. 生態(tài)學(xué)報(bào), 2016, 36(4): 1021- 1029.

    [13] RautiainenM, StenbergP, NilsonT. Estimating canopy cover in Scots pine stands. Silva Fennica, 2005, 39(1): 137- 142

    [14] Lan G Y, Getzin S, Wiegand T, Hu Y H, Xie G S, Zhu H, Cao M. Spatial distribution and interspecific associations of tree species in a tropical seasonal rain forest of China. PLoS One, 2010, 7(11): e46074.

    [15] 蔡飛. 杭州西湖山區(qū)青岡種群結(jié)構(gòu)和動(dòng)態(tài)的研究. 林業(yè)科學(xué), 2000, 36(3): 67- 72.

    [16] Frost I, Rydin H. Spatial pattern and size distribution of the animal-dispersed treeQuercusroburin two spruce-dominated forests. écoscience, 2000, 7(1): 38- 44.

    [17] 周建云, 李榮, 張文輝, 何景峰. 不同間伐強(qiáng)度下遼東櫟種群結(jié)構(gòu)特征與空間分布格局. 林業(yè)科學(xué), 2012, 48(4): 149- 155.

    [18] 王曉雨, 于大炮, 周莉, 周旺明, 吳志軍, 郭焱, 包也, 孟瑩瑩, 代力民. 長(zhǎng)白山北坡林線岳樺種群空間分布格局. 生態(tài)學(xué)報(bào), 2015, 35(1): 116- 124.

    [19] Ripley B D. Modelling spatial patterns.Journal of the Royal statistical Society.Series B, 1977, 39(2): 172- 212.

    [20] Wang X T, Liang C Z, Wang W. Balance between facilitation and competition determines spatial patterns in a plant population. Chinese Science Bulletin, 2014, 59(13): 1405- 1415.

    [21] BaddeleyA, Turner B. Spatstat: an R package for analyzing spatial point patterns. Journal of Statistical Software, 2005, 12(6): 1- 42.

    [22] Benson E J, Hartnett D C. The role of seed and vegetative reproduction in plant recruitment and demography in tallgrass prairie.Plant Ecology, 2006, 187(2): 163- 178.

    [23] Wang XG, Ye J, Li BH, Zhang J, Lin F, Hao ZQ. Spatial distributions of species in an old-growth temperate forest, northeastern China. Canadian Journal of Forest Research, 2010, 40(6): 1011- 1019.

    [24] 蔡新斌, 吳俊俠. 甘家湖自然保護(hù)區(qū)白梭梭種群特征與動(dòng)態(tài)分析. 干旱區(qū)資源與環(huán)境, 2016, 30(7): 90- 94.

    [25] 閆興富, 杜茜, 石淳, 周立彪, 張靠穩(wěn). 六盤山區(qū)遼東櫟的實(shí)生苗更新及其影響因子. 植物生態(tài)學(xué)報(bào), 2011, 35(9): 914- 925.

    [26] 趙學(xué)農(nóng), 劉倫輝, 高圣義, 和愛(ài)軍. 版納青梅種群結(jié)構(gòu)動(dòng)態(tài)與分布格局. 植物學(xué)報(bào), 1993, 35(7): 552- 560.

    [27] Li X M, Zhou W Q, Ouyang Z Y. Relationship between land surface temperature and spatial pattern of greenspace: what are the effects of spatial resolution?.Landscape and Urban Planning, 2013, 114(8): 1- 8.

    [28] 龍成, 楊小波, 龍文興, 李東海. 銅鼓嶺熱帶常綠季雨矮林5種蒲桃屬植物的種群結(jié)構(gòu)及空間格局. 林業(yè)科學(xué), 2015, 51(2): 18- 27.

    [29] 劉振, 李紅麗, 董智, 李鋼鐵, 萬(wàn)玲玲, 岳永杰. 渾善達(dá)克沙地2種生境下榆樹種群空間點(diǎn)格局. 林業(yè)科學(xué), 2012, 48(1): 29- 34.

    [30] Hubbell S P. Neutral theory and the evolution of ecological equivalence. Ecology, 2006, 87(6): 1387- 1398.

    [31] Greig-Smith P. Quantitative Plant Ecology. 3rd ed. Oxford: Blackwell Scientific Publications, 1983.

    [32] Wiegand T, Gunatilleke S, Gunatilleke N, Okuda T. Analyzing the spatial structure of a Sri Lankan tree species with multiple scales of clustering. Ecology, 2007, 88(12): 3088- 3102.

    [33] 張希彪, 郭小強(qiáng), 上官周平, 劉富順. 黃土丘陵子午嶺油松天然林群落特征研究. 植物研究, 2006, 26(2): 169- 175.

    [34] 侯繼華, 黃建輝, 馬克平. 東靈山遼東櫟林主要樹種種群11年動(dòng)態(tài)變化. 植物生態(tài)學(xué)報(bào), 2004, 28(5): 609- 615.

    [35] 尤海舟, 劉興良, 繆寧, 何飛, 馬欽彥. 川滇高山櫟種群不同海拔空間格局的尺度效應(yīng)及個(gè)體間空間關(guān)聯(lián). 生態(tài)學(xué)報(bào), 2010, 30(15): 4004- 4011.

    [36] Kenkel N C. Pattern of self-thinning in Jack Pine: testing the random mortality hypothesis. Ecology, 1988, 69(4): 1017- 1024.

    Spatial pattern and spatial association ofQuercusacutissimaat different developmental stages in the Qiaoshan Mountains

    XUE Wenyan, YANG Bin, ZHANG Wenhui*, YU Shichuan

    KeyComprehensiveLaboratoryofForestryinShaanxiProvince,NorthwestA&FUniversity,Yangling712100,China

    Spatial pattern and spatial association analysis are important in the study of population characteristics and interspecies relationships. To investigate the spatial distribution pattern and association in different developmental stages in aQuercusacutissimapopulation, the natural secondaryQ.acutissimaforest of the Qiaoshan Mountains in the southern area of the Loess Plateau was used as our study area. One forest plot (50 m×50 m) was established in eachQ.acutissimanatural secondary pure forest (plot A) andQ.acutissima-Pinustabuliformisnatural secondary mixed forest (plot B). The coordinate of eachQ.acutissimawas mapped, the diameter at breast height and height of each individual were measured by plot to investigate means. We divided the individuals into four growth stages using the method of diameter class. The spatial distribution pattern and associations were analyzed using Ripley′sL(r) Function. The results showed the followings: (1)The quantity structure of each developmental stage of theQ.acutissimapopulation in the Qiaoshan Mountains exhibited an irregular inverted "J" type, indicating a stable population structure with persistent regeneration capability. (2) The spatial distribution of the two populations was greatly influenced by habitat heterogeneity. Spatial distribution exhibited high aggregation at the scale of 0—10 m in plot A, whereas it showed low aggregation at the scale of 5—35 m. In plot B, the degree of aggregation changed quickly in response to habitat heterogeneity. It exhibited high aggregation at the 0—5 m scale and low aggregation at the 5—35 m scale. (3) With the development of populations and increase in scale, the spatial aggregation of theQ.acutissimapopulation was reduced. In plot A, the seedling stage showed an aggregated distribution at the small scale of 0—6 m and a random distribution at the larger scale of 6—12 m. The sapling stage showed an aggregated distribution at all scales. The mid-mature stage and mature stage showed a random distribution at all scales. In plot B, all stages showed a random distribution at all the scales, except for the aggregated distribution of the seedling stage at the small scale of 0—6 m. (4) The spatial correlation between different growth stages and scales was discriminating in the two plots. In plot A, we found a negative correlation between the seedling and sapling stage at the 1—2 m scale. Positive correlation was found between the sapling and mid-mature stage at the 8 m scale. However, the relationship among other stages was not significant at any scale. In plot B, there was a significant positive correlation between the sapling and mid-mature stage at the 2—12 m scale. There was almost no correlation among other stages at any scales. Population distribution patterns were influenced by developmental stages, biological characteristics, and ecological factors. In future forest management, different managerial measures need to be taken to promote the regeneration of the population according to the different developmental stages and different habitats of the population.

    Quercusacutissima; spatial pattern; spatial association; developmental stages

    國(guó)家林業(yè)局全國(guó)森林經(jīng)營(yíng)基礎(chǔ)研究項(xiàng)目(1692016-03);黃土高原松櫟林多功能近自然經(jīng)營(yíng)作業(yè)法及其效益監(jiān)測(cè)研究項(xiàng)目(1692016-03)

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

    10.5846/stxb201603130447

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

    薛文艷, 楊斌, 張文輝, 于世川.橋山林區(qū)麻櫟種群不同發(fā)育階段空間格局及關(guān)聯(lián)性.生態(tài)學(xué)報(bào),2017,37(10):3375- 3384.

    Xue W Y, Yang B, Zhang W H, Yu S C.Spatial pattern and spatial association ofQuercusacutissimaat different developmental stages in the Qiaoshan Mountains.Acta Ecologica Sinica,2017,37(10):3375- 3384.

    猜你喜歡
    麻櫟發(fā)育階段樣地
    麻櫟果實(shí)的開(kāi)發(fā)利用及培育探究
    小麥生殖發(fā)育階段對(duì)低溫的敏感性鑒定
    額爾古納市興安落葉松中齡林植被碳儲(chǔ)量研究
    昆明市主要綠化樹種閾值測(cè)定與分析
    基于角尺度模型的林業(yè)樣地空間結(jié)構(gòu)分析
    麻櫟種源1年生苗生長(zhǎng)特性比較
    定向培育措施對(duì)麻櫟生長(zhǎng)的影響評(píng)價(jià)
    對(duì)森工林區(qū)在商品林基地培養(yǎng)速生楊樹探討
    大花黃牡丹葉片發(fā)育過(guò)程中氣孔密度和氣孔指數(shù)的動(dòng)態(tài)變化
    廣西石漠化地區(qū)不同植被恢復(fù)樣地的土壤微生物量比較
    狂野欧美激情性xxxx在线观看| 亚洲精品色激情综合| 噜噜噜噜噜久久久久久91| 精品久久久精品久久久| 色视频www国产| 日本一本二区三区精品| 成年版毛片免费区| 亚洲精品日韩在线中文字幕| 国产亚洲午夜精品一区二区久久 | 精品少妇黑人巨大在线播放| 国产永久视频网站| 在线观看三级黄色| 欧美xxxx黑人xx丫x性爽| 如何舔出高潮| 中国美白少妇内射xxxbb| 麻豆久久精品国产亚洲av| 久久久久久久久大av| 久久6这里有精品| 成年免费大片在线观看| 欧美高清成人免费视频www| 国产免费福利视频在线观看| 免费av毛片视频| 男男h啪啪无遮挡| 晚上一个人看的免费电影| 亚洲国产精品999| 亚洲激情五月婷婷啪啪| 亚洲精品成人久久久久久| 国产高清国产精品国产三级 | 精品人妻一区二区三区麻豆| 亚洲成色77777| 少妇丰满av| 精品少妇久久久久久888优播| 国产乱来视频区| 国产午夜精品一二区理论片| 交换朋友夫妻互换小说| 国产黄a三级三级三级人| 一本一本综合久久| 国产色婷婷99| 三级国产精品欧美在线观看| 国产成人精品婷婷| 国产淫片久久久久久久久| 国产成人精品一,二区| 国产有黄有色有爽视频| 18禁裸乳无遮挡动漫免费视频 | 22中文网久久字幕| 国产黄a三级三级三级人| 麻豆成人午夜福利视频| 国产探花在线观看一区二区| 免费大片黄手机在线观看| 小蜜桃在线观看免费完整版高清| 国产黄片美女视频| 国产探花在线观看一区二区| 亚洲最大成人中文| 国内精品美女久久久久久| 97人妻精品一区二区三区麻豆| 亚洲av免费在线观看| 99九九线精品视频在线观看视频| 亚洲高清免费不卡视频| 国产中年淑女户外野战色| 国产久久久一区二区三区| 青春草视频在线免费观看| 亚洲在久久综合| 亚洲欧美精品自产自拍| 亚洲精品成人久久久久久| 精品久久国产蜜桃| 大片免费播放器 马上看| 狠狠精品人妻久久久久久综合| 免费看光身美女| 亚洲精品日韩av片在线观看| 国精品久久久久久国模美| 欧美成人精品欧美一级黄| 国产精品99久久久久久久久| 国产亚洲5aaaaa淫片| av在线天堂中文字幕| 三级国产精品片| 国产毛片在线视频| 亚洲成人中文字幕在线播放| 久久久久久久精品精品| 亚洲国产成人一精品久久久| 人妻制服诱惑在线中文字幕| 少妇人妻 视频| 嫩草影院入口| 97在线人人人人妻| h日本视频在线播放| 精品人妻偷拍中文字幕| 国产高清不卡午夜福利| 久久精品国产亚洲网站| 欧美 日韩 精品 国产| 如何舔出高潮| 精品久久久噜噜| 欧美激情久久久久久爽电影| 男女边摸边吃奶| 一级毛片 在线播放| 免费大片18禁| 男人和女人高潮做爰伦理| 国产v大片淫在线免费观看| 国产欧美日韩精品一区二区| 久久鲁丝午夜福利片| freevideosex欧美| 色吧在线观看| 欧美少妇被猛烈插入视频| 性插视频无遮挡在线免费观看| 超碰97精品在线观看| 日韩视频在线欧美| 精品少妇黑人巨大在线播放| 日本免费在线观看一区| 国产黄频视频在线观看| 各种免费的搞黄视频| 欧美zozozo另类| 男人爽女人下面视频在线观看| 国产一级毛片在线| 联通29元200g的流量卡| 久久女婷五月综合色啪小说 | 精品国产三级普通话版| 久久久久久久午夜电影| 精品久久久精品久久久| 欧美日韩亚洲高清精品| 在线播放无遮挡| 国产永久视频网站| 观看美女的网站| 久久久久久久大尺度免费视频| 欧美最新免费一区二区三区| xxx大片免费视频| 黄色视频在线播放观看不卡| 国产毛片在线视频| 日本wwww免费看| 99久久精品国产国产毛片| 成人一区二区视频在线观看| 日产精品乱码卡一卡2卡三| 午夜亚洲福利在线播放| videos熟女内射| 午夜福利网站1000一区二区三区| 亚洲av电影在线观看一区二区三区 | 中文字幕久久专区| 亚洲av中文av极速乱| 国产综合精华液| 你懂的网址亚洲精品在线观看| 日韩大片免费观看网站| 国产美女午夜福利| 欧美激情国产日韩精品一区| 最新中文字幕久久久久| 热99国产精品久久久久久7| 可以在线观看毛片的网站| 婷婷色麻豆天堂久久| 国产精品女同一区二区软件| av线在线观看网站| av.在线天堂| 男女边吃奶边做爰视频| 国产精品久久久久久久久免| 精品国产露脸久久av麻豆| 九色成人免费人妻av| 欧美3d第一页| 大陆偷拍与自拍| 一级a做视频免费观看| 久久99热6这里只有精品| 国产探花极品一区二区| 高清日韩中文字幕在线| 欧美成人a在线观看| 久久久久久久精品精品| 亚洲av成人精品一二三区| 亚洲精品色激情综合| 舔av片在线| 国产在视频线精品| 午夜爱爱视频在线播放| 国产欧美日韩一区二区三区在线 | 久久精品综合一区二区三区| 国产精品久久久久久久久免| av.在线天堂| 欧美另类一区| 18禁动态无遮挡网站| 国产av码专区亚洲av| 最近最新中文字幕免费大全7| 亚洲欧美中文字幕日韩二区| 亚洲欧美日韩东京热| 亚洲精品乱久久久久久| 亚洲国产色片| 联通29元200g的流量卡| 日韩一本色道免费dvd| 免费av毛片视频| 六月丁香七月| 欧美3d第一页| 夜夜看夜夜爽夜夜摸| 狂野欧美激情性xxxx在线观看| 免费黄网站久久成人精品| 人妻夜夜爽99麻豆av| 青青草视频在线视频观看| 亚洲不卡免费看| 两个人的视频大全免费| 女人十人毛片免费观看3o分钟| 国产成人精品一,二区| av播播在线观看一区| av国产久精品久网站免费入址| 精品视频人人做人人爽| 99热国产这里只有精品6| 在现免费观看毛片| 99re6热这里在线精品视频| 丰满人妻一区二区三区视频av| 免费黄色在线免费观看| 成人亚洲欧美一区二区av| 亚洲国产欧美人成| 九九爱精品视频在线观看| 欧美精品人与动牲交sv欧美| 一二三四中文在线观看免费高清| 国产淫语在线视频| 视频区图区小说| 亚洲欧美清纯卡通| 欧美丝袜亚洲另类| 性色avwww在线观看| 亚洲av不卡在线观看| 三级经典国产精品| 国产精品久久久久久精品电影| 啦啦啦中文免费视频观看日本| 国产精品国产三级国产专区5o| 国产探花极品一区二区| av播播在线观看一区| 女人十人毛片免费观看3o分钟| 伊人久久国产一区二区| 91狼人影院| 中文欧美无线码| 亚洲av电影在线观看一区二区三区 | 最近的中文字幕免费完整| av专区在线播放| 网址你懂的国产日韩在线| 菩萨蛮人人尽说江南好唐韦庄| 插阴视频在线观看视频| 最近最新中文字幕大全电影3| 日本免费在线观看一区| 七月丁香在线播放| 成人漫画全彩无遮挡| 亚洲欧美日韩另类电影网站 | 国产成人精品福利久久| 国产亚洲最大av| 中文字幕人妻熟人妻熟丝袜美| 男人舔奶头视频| 成人无遮挡网站| 欧美日本视频| 各种免费的搞黄视频| 蜜桃亚洲精品一区二区三区| 日韩欧美 国产精品| 国产免费视频播放在线视频| 亚洲精华国产精华液的使用体验| av又黄又爽大尺度在线免费看| 18禁动态无遮挡网站| 久久热精品热| 久久国内精品自在自线图片| 三级经典国产精品| 在线观看美女被高潮喷水网站| 黄色欧美视频在线观看| 视频区图区小说| 久久人人爽人人片av| 欧美区成人在线视频| 色视频在线一区二区三区| 久久精品熟女亚洲av麻豆精品| 免费播放大片免费观看视频在线观看| 亚洲精品第二区| 国产精品久久久久久久电影| 久久6这里有精品| 另类亚洲欧美激情| 我的女老师完整版在线观看| 欧美日韩亚洲高清精品| 一区二区三区精品91| 好男人在线观看高清免费视频| 色综合色国产| 18禁在线播放成人免费| 久久久精品免费免费高清| 免费电影在线观看免费观看| 晚上一个人看的免费电影| 日韩av在线免费看完整版不卡| 97热精品久久久久久| av在线亚洲专区| 自拍欧美九色日韩亚洲蝌蚪91 | 欧美3d第一页| 午夜日本视频在线| 日韩亚洲欧美综合| 最近中文字幕2019免费版| 日本av手机在线免费观看| 男人添女人高潮全过程视频| 国产精品久久久久久精品电影| 亚洲精品456在线播放app| 超碰97精品在线观看| 久久热精品热| 中文乱码字字幕精品一区二区三区| 不卡视频在线观看欧美| 男女国产视频网站| 成年女人在线观看亚洲视频 | 80岁老熟妇乱子伦牲交| 亚洲精品乱久久久久久| av黄色大香蕉| 少妇丰满av| 人妻系列 视频| 身体一侧抽搐| 最近手机中文字幕大全| 国产黄片美女视频| 女人十人毛片免费观看3o分钟| 亚洲精品一二三| 亚洲内射少妇av| 国产大屁股一区二区在线视频| 亚洲熟女精品中文字幕| 建设人人有责人人尽责人人享有的 | 色5月婷婷丁香| 亚洲欧洲国产日韩| 欧美日韩视频精品一区| 国产精品秋霞免费鲁丝片| 亚洲高清免费不卡视频| 蜜桃久久精品国产亚洲av| 亚洲av二区三区四区| 最新中文字幕久久久久| 一二三四中文在线观看免费高清| 中文字幕人妻熟人妻熟丝袜美| 在线观看av片永久免费下载| 国产精品嫩草影院av在线观看| av天堂中文字幕网| 丝袜美腿在线中文| 十八禁网站网址无遮挡 | av在线观看视频网站免费| 在线 av 中文字幕| 国产午夜福利久久久久久| videossex国产| 国产成人freesex在线| 欧美另类一区| a级一级毛片免费在线观看| 国产一区亚洲一区在线观看| 中文在线观看免费www的网站| 国产一区二区在线观看日韩| 欧美成人精品欧美一级黄| 国产精品成人在线| 国产视频内射| 国产精品国产三级专区第一集| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 亚洲欧美一区二区三区国产| 日韩人妻高清精品专区| 黄色视频在线播放观看不卡| 午夜福利高清视频| 人妻少妇偷人精品九色| 边亲边吃奶的免费视频| 下体分泌物呈黄色| 午夜福利网站1000一区二区三区| 亚洲国产精品成人久久小说| 91久久精品电影网| 熟女人妻精品中文字幕| 91午夜精品亚洲一区二区三区| 国产淫片久久久久久久久| 国产成人精品久久久久久| 丰满少妇做爰视频| av福利片在线观看| 能在线免费看毛片的网站| 夫妻性生交免费视频一级片| 国模一区二区三区四区视频| 熟女电影av网| 国产亚洲av嫩草精品影院| 国产有黄有色有爽视频| 毛片女人毛片| 在线看a的网站| 2021少妇久久久久久久久久久| 全区人妻精品视频| 亚洲av免费高清在线观看| 大陆偷拍与自拍| 欧美xxxx黑人xx丫x性爽| 波野结衣二区三区在线| 最后的刺客免费高清国语| 成人亚洲精品一区在线观看 | 国产淫片久久久久久久久| 人妻夜夜爽99麻豆av| 在线 av 中文字幕| 精品人妻偷拍中文字幕| 韩国av在线不卡| 日韩av在线免费看完整版不卡| 乱系列少妇在线播放| 久久久久久九九精品二区国产| 国产老妇伦熟女老妇高清| 寂寞人妻少妇视频99o| 午夜福利视频1000在线观看| 国产午夜福利久久久久久| 国内精品宾馆在线| 22中文网久久字幕| a级毛色黄片| 只有这里有精品99| 极品教师在线视频| 男女下面进入的视频免费午夜| 午夜精品国产一区二区电影 | 欧美丝袜亚洲另类| 国产一区二区三区av在线| 建设人人有责人人尽责人人享有的 | 国产av码专区亚洲av| 最近手机中文字幕大全| 亚洲精品aⅴ在线观看| 黄色视频在线播放观看不卡| 蜜臀久久99精品久久宅男| a级毛色黄片| 亚洲成人精品中文字幕电影| 如何舔出高潮| 菩萨蛮人人尽说江南好唐韦庄| 精品人妻一区二区三区麻豆| 777米奇影视久久| 欧美潮喷喷水| 狂野欧美激情性xxxx在线观看| 中文资源天堂在线| 国产av国产精品国产| 国产精品.久久久| 91精品伊人久久大香线蕉| 亚洲经典国产精华液单| 久久久久九九精品影院| 永久免费av网站大全| 简卡轻食公司| 精品人妻偷拍中文字幕| 欧美3d第一页| 一区二区av电影网| 国产成人精品一,二区| 日本猛色少妇xxxxx猛交久久| a级一级毛片免费在线观看| 搡老乐熟女国产| 久久亚洲国产成人精品v| 黄色怎么调成土黄色| 人妻制服诱惑在线中文字幕| 一区二区三区免费毛片| 亚洲第一区二区三区不卡| 男人和女人高潮做爰伦理| 久久久久久九九精品二区国产| 秋霞在线观看毛片| 日韩中字成人| 国产 一区 欧美 日韩| 在线观看三级黄色| 国产在线男女| 少妇裸体淫交视频免费看高清| 看免费成人av毛片| 美女内射精品一级片tv| 亚洲欧美一区二区三区黑人 | 免费观看的影片在线观看| 国内精品美女久久久久久| 最近2019中文字幕mv第一页| 少妇人妻一区二区三区视频| 久久久久久久久久久丰满| 亚州av有码| 可以在线观看毛片的网站| 搡老乐熟女国产| 国产一区二区三区综合在线观看 | av黄色大香蕉| 精品久久国产蜜桃| 久久精品国产亚洲av涩爱| 一个人观看的视频www高清免费观看| 免费av不卡在线播放| 男人和女人高潮做爰伦理| 亚洲激情五月婷婷啪啪| 特大巨黑吊av在线直播| 丰满人妻一区二区三区视频av| 成年女人在线观看亚洲视频 | 老司机影院毛片| 午夜福利视频1000在线观看| 亚洲精品视频女| 亚洲精品456在线播放app| 日韩制服骚丝袜av| 欧美潮喷喷水| 久久久久久久精品精品| 国产有黄有色有爽视频| 久久久久国产精品人妻一区二区| 精品久久久久久久久亚洲| av线在线观看网站| 亚洲无线观看免费| 一区二区三区四区激情视频| 高清午夜精品一区二区三区| 男人和女人高潮做爰伦理| 国产免费一级a男人的天堂| 乱系列少妇在线播放| 亚洲自拍偷在线| 中文字幕免费在线视频6| 亚洲精品乱久久久久久| 美女被艹到高潮喷水动态| 激情五月婷婷亚洲| 日本熟妇午夜| 99久久精品国产国产毛片| 日韩av在线免费看完整版不卡| 嘟嘟电影网在线观看| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看| 卡戴珊不雅视频在线播放| 免费观看av网站的网址| 青春草国产在线视频| 又大又黄又爽视频免费| 成人综合一区亚洲| 欧美日韩视频高清一区二区三区二| 国产白丝娇喘喷水9色精品| 色婷婷久久久亚洲欧美| 交换朋友夫妻互换小说| 中文乱码字字幕精品一区二区三区| 欧美三级亚洲精品| av国产精品久久久久影院| 国产毛片在线视频| 国产精品爽爽va在线观看网站| 国产成人a区在线观看| 少妇人妻 视频| 欧美丝袜亚洲另类| 午夜福利网站1000一区二区三区| 国模一区二区三区四区视频| 五月开心婷婷网| 免费大片18禁| 91精品一卡2卡3卡4卡| 国国产精品蜜臀av免费| 可以在线观看毛片的网站| 久久人人爽人人爽人人片va| 国产免费又黄又爽又色| 国产精品人妻久久久影院| 国产一区有黄有色的免费视频| 亚洲伊人久久精品综合| 丝瓜视频免费看黄片| eeuss影院久久| 久热久热在线精品观看| 日本午夜av视频| 国产v大片淫在线免费观看| 汤姆久久久久久久影院中文字幕| 汤姆久久久久久久影院中文字幕| 麻豆成人午夜福利视频| 交换朋友夫妻互换小说| 亚洲成人久久爱视频| 日本午夜av视频| 天堂网av新在线| 欧美日韩视频高清一区二区三区二| 最近中文字幕2019免费版| 一级黄片播放器| 免费观看a级毛片全部| 国产高清三级在线| 综合色丁香网| 亚洲第一区二区三区不卡| 美女高潮的动态| 成人国产av品久久久| 高清av免费在线| 亚洲欧美日韩东京热| 在现免费观看毛片| 美女主播在线视频| 精品一区二区三卡| 日本爱情动作片www.在线观看| 国产精品久久久久久久电影| 欧美丝袜亚洲另类| 免费不卡的大黄色大毛片视频在线观看| 久久久欧美国产精品| 六月丁香七月| freevideosex欧美| 亚洲天堂国产精品一区在线| 欧美性感艳星| 永久网站在线| 国产一级毛片在线| 一级片'在线观看视频| 人妻制服诱惑在线中文字幕| 国产精品久久久久久精品电影小说 | 一级毛片 在线播放| 熟女av电影| 婷婷色综合www| 日韩欧美精品v在线| 丝袜喷水一区| 久久99热这里只频精品6学生| 身体一侧抽搐| 亚洲国产成人一精品久久久| 国产一区二区亚洲精品在线观看| 女人久久www免费人成看片| 亚洲欧美中文字幕日韩二区| 少妇人妻 视频| 看黄色毛片网站| 18禁裸乳无遮挡免费网站照片| 国产黄频视频在线观看| 国产成人a区在线观看| 久久精品久久精品一区二区三区| 91久久精品电影网| 99久久人妻综合| 久久影院123| 欧美日韩一区二区视频在线观看视频在线 | 91精品一卡2卡3卡4卡| 亚洲欧美清纯卡通| 97人妻精品一区二区三区麻豆| 青春草视频在线免费观看| 少妇猛男粗大的猛烈进出视频 | 熟妇人妻不卡中文字幕| 欧美少妇被猛烈插入视频| av福利片在线观看| 国产69精品久久久久777片| 自拍欧美九色日韩亚洲蝌蚪91 | 国产毛片在线视频| 色婷婷久久久亚洲欧美| 精品国产三级普通话版| 91精品一卡2卡3卡4卡| 国产亚洲av嫩草精品影院| 精品久久久久久久久亚洲| 午夜视频国产福利| 亚洲综合精品二区| 成人黄色视频免费在线看| 国产淫语在线视频| 好男人视频免费观看在线| 自拍偷自拍亚洲精品老妇| 黄色欧美视频在线观看| 国产精品久久久久久av不卡| 免费看a级黄色片| 夜夜爽夜夜爽视频| 性插视频无遮挡在线免费观看| 久久99精品国语久久久| 欧美97在线视频| 日韩av免费高清视频| av天堂中文字幕网| 又爽又黄无遮挡网站| 精品久久国产蜜桃| 欧美区成人在线视频| 国产v大片淫在线免费观看| 久久久久久久午夜电影| 亚洲电影在线观看av| 在线亚洲精品国产二区图片欧美 | 网址你懂的国产日韩在线| 少妇熟女欧美另类| 各种免费的搞黄视频| 白带黄色成豆腐渣| 色视频www国产| 少妇的逼好多水| 欧美区成人在线视频| 欧美 日韩 精品 国产| 国产精品爽爽va在线观看网站| 亚洲欧美精品自产自拍| 丝瓜视频免费看黄片| 综合色丁香网|