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

    珍稀瀕危植物江南油杉群落喬木層主要樹(shù)種種間關(guān)聯(lián)性分析

    2023-03-09 02:14:32陳緒輝葉寶鑒潘標(biāo)志林碧華羅敏賢肖麗芳劉寶何宗明鄭世群
    關(guān)鍵詞:油杉種間喬木

    陳緒輝, 葉寶鑒, 潘標(biāo)志, 林碧華, 羅敏賢, 肖麗芳, 劉寶, 何宗明, 鄭世群*

    珍稀瀕危植物江南油杉群落喬木層主要樹(shù)種種間關(guān)聯(lián)性分析

    陳緒輝1a,1b, 葉寶鑒2, 潘標(biāo)志3, 林碧華1a,1b, 羅敏賢1a,1b, 肖麗芳1a,1b, 劉寶1a, 何宗明1a, 鄭世群1a,1b*

    (1. 福建農(nóng)林大學(xué), a. 林學(xué)院;b. 海峽自然保護(hù)區(qū)研究中心,生態(tài)與資源統(tǒng)計(jì)福建省高校重點(diǎn)實(shí)驗(yàn)室,福州 350002;2.福建工程學(xué)院建筑與城鄉(xiāng)規(guī)劃學(xué)院,福州 350118;3.福建省林業(yè)調(diào)查規(guī)劃院,福州 350003)

    基于2×2列聯(lián)表,采用方差比率法(v)、2檢驗(yàn)、Pearson相關(guān)系數(shù)檢驗(yàn)和Spearman秩相關(guān)系數(shù)檢驗(yàn)等方法,對(duì)福建君子峰國(guó)家級(jí)自然保護(hù)區(qū)江南油杉(var)群落喬木層主要樹(shù)種進(jìn)行種間關(guān)聯(lián)性分析,并通過(guò)主成分分析(PCA)劃分生態(tài)種組。結(jié)果表明,21種主要樹(shù)種總體表現(xiàn)為不顯著負(fù)關(guān)聯(lián),且江南油杉與其他物種多呈不顯著相關(guān)。210個(gè)種對(duì)中,正相關(guān)的有110對(duì),負(fù)相關(guān)的有81對(duì),不相關(guān)的有19對(duì),正負(fù)比1.36。絕大多數(shù)種對(duì)間呈不顯著相關(guān),顯著率僅為0.48%,物種間相互獨(dú)立。Pearson相關(guān)系數(shù)和Spearman秩相關(guān)系數(shù)檢驗(yàn)結(jié)果具有一致性,負(fù)相關(guān)種對(duì)數(shù)大于正相關(guān)種對(duì)數(shù),絕大多數(shù)種對(duì)不顯著,顯著率分別為6.19%和7.14%。與2檢驗(yàn)相比,Pearson相關(guān)系數(shù)和Spearman秩相關(guān)系數(shù)檢驗(yàn)方式更加靈敏。運(yùn)用PCA二維排序可將主要樹(shù)種劃分為4個(gè)生態(tài)種組,同一生態(tài)種組內(nèi)物種的生態(tài)習(xí)性和資源需求相似,多呈正相關(guān),種對(duì)間相互依賴共存;不同生態(tài)種組間物種多呈負(fù)相關(guān),種對(duì)間相互競(jìng)爭(zhēng)。江南油杉群落喬木層主要樹(shù)種的種間關(guān)聯(lián)性較為松散,大多數(shù)物種為獨(dú)立分布,幼苗存在更新障礙,且幼樹(shù)個(gè)體競(jìng)爭(zhēng)壓力較大,種群規(guī)模呈逐漸縮小的趨勢(shì)。為促進(jìn)江南油杉種群正常生長(zhǎng)與更新,應(yīng)對(duì)江南油杉幼苗進(jìn)行重點(diǎn)保護(hù)。

    江南油杉;珍稀瀕危植物;種間關(guān)聯(lián);生態(tài)種組

    種間關(guān)系包括種間關(guān)聯(lián)和種間相關(guān)性,兩者都是指不同種群在空間上的相互關(guān)聯(lián)性[1],一直以來(lái)是生態(tài)學(xué)研究的熱點(diǎn)。種間關(guān)聯(lián)是植物群落中重要的結(jié)構(gòu)和數(shù)量特征之一,是不同種群間相互關(guān)系的表達(dá)形式以及群落形成和演替的基礎(chǔ)[2]。研究不同種群間的關(guān)聯(lián)性,有利于正確認(rèn)識(shí)各物種間的競(jìng)爭(zhēng)關(guān)系,客觀反映物種對(duì)環(huán)境的適應(yīng)程度,對(duì)認(rèn)識(shí)群落內(nèi)部結(jié)構(gòu)、預(yù)測(cè)種群動(dòng)態(tài)及群落演替趨勢(shì)有重要意義[3]。生態(tài)種組是群落中生活型相似物種的集合,劃分生態(tài)種組有利于探究不同種群之間的生境隔離機(jī)制,對(duì)認(rèn)識(shí)種群間的連通性和維護(hù)生態(tài)完整性具有重要指導(dǎo)作用[4]。近年來(lái),人們對(duì)閩樺()[5]、柳杉()[6]和麻櫟()[7]等植物群落的種間關(guān)聯(lián)性進(jìn)行研究,并劃分生態(tài)種組,揭示這些植物在群落中的功能地位和種群動(dòng)態(tài),為該種群保護(hù)提供了科學(xué)依據(jù)。

    江南油杉(var)系松科(Pinaceae)油杉屬常綠喬木,是油杉()的變種,為中國(guó)特有樹(shù)種,國(guó)家二級(jí)重點(diǎn)保護(hù)野生植物[8],主要分布在我國(guó)東南、華南和西南地區(qū), 大多呈散生或局部分布[9]。喜光、較耐旱,根系發(fā)達(dá),樹(shù)干通直,材質(zhì)耐腐,為建筑、橋梁和家具等優(yōu)質(zhì)用材,具有很高的經(jīng)濟(jì)價(jià)值[10]。目前對(duì)江南油杉的研究主要集中在種源[11]、幼苗[12]、細(xì)根生物量[13]和基因[14]等方面,由于野生種群數(shù)量少,生態(tài)學(xué)相關(guān)研究甚為欠缺,僅李強(qiáng)等[15]和羅金旺等[16]對(duì)其群落結(jié)構(gòu)和物種多樣性方面有初步探討,難以反映出種群間的共生關(guān)系;另一方面由于生物學(xué)特性較特殊,加上群落干擾嚴(yán)重及環(huán)境變化等多方面影響, 導(dǎo)致江南油杉處于瀕危狀態(tài),因此亟需開(kāi)展種群和群落生態(tài)學(xué)研究,對(duì)其野生種群進(jìn)行保護(hù)。福建君子峰國(guó)家級(jí)自然保護(hù)區(qū)殘存一處保存較完好的江南油杉天然群落,群落總面積約為0.6 hm2,以6株江南油杉古樹(shù)為主,冠幅大,樹(shù)干通直,平均胸徑113.3 cm,平均樹(shù)高達(dá)27.5 m,最大樹(shù)齡700 a。本文研究其喬木層主要樹(shù)種的種間關(guān)聯(lián)性,劃分生態(tài)種組,并分析江南油杉種群年齡結(jié)構(gòu),旨在探明江南油杉與其他樹(shù)種的種間關(guān)系及群落演替趨勢(shì), 從而為制定江南油杉種群保護(hù)措施提供理論依據(jù)。

    1 材料和方法

    1.1 研究區(qū)概況

    福建君子峰國(guó)家級(jí)自然保護(hù)區(qū)地處明溪縣,區(qū)內(nèi)江南油杉僅有零星幾個(gè)分布點(diǎn),種群數(shù)量稀少,多不成群落,唯有夏陽(yáng)鄉(xiāng)紫云村的羅坊水尾(26° 21′00″ N,117°28′51″ E)江南油杉群落保存較完好。該區(qū)域?yàn)榈湫偷牡蜕降孛玻0?54~756 m, 地勢(shì)較平坦。屬中亞熱帶海洋性季風(fēng)氣候,年均溫18 ℃,極端最高溫39.1 ℃,極端最低溫–8 ℃,年均降水量1 737 mm,無(wú)霜期260 d。土壤為花崗巖發(fā)育成的微酸性土壤,土層較深厚,腐殖質(zhì)層4~8 cm。群落分層明顯,江南油杉位居第一主林層。喬木層以江南油杉、杉木()、甜櫧()、馬尾松()、赤楊葉()和南酸棗()等樹(shù)種為主;灌木層種類豐富,以箬竹()為主,其蓋度達(dá)70%,還有紅皮糙果茶()、光葉山礬()、褐毛石楠()和杜莖山()等樹(shù)種。草本層較單調(diào),主要由里白()、狗脊()和薹草屬()植物等組成。

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

    根據(jù)群落范圍特點(diǎn),設(shè)置6個(gè)相鄰樣地(2排×3列),按“S”形依次編號(hào)為Ⅰ~Ⅵ,每個(gè)樣地為20 m× 20 m,總面積2 400 m2。每個(gè)樣地劃分為4個(gè)10 m× 10 m喬木樣方,記錄喬木樹(shù)種(胸徑≥5 cm)的種名、株數(shù)、胸徑、樹(shù)高、冠幅和坐標(biāo)等信息。在每個(gè)樣地四角和中心處取5個(gè)5 m×5 m的灌木樣方和5個(gè)2 m×2 m的草本樣方,調(diào)查并記錄每個(gè)小樣方中所有植物的種名、株(叢)數(shù)、高度和蓋度等。計(jì)算喬木層各物種在樣地中的重要值,選取重要值大于1的物種為研究對(duì)象。

    1.3 數(shù)據(jù)處理

    1.3.1 重要值計(jì)算

    重要值是表征物種在群落中作用和地位的指標(biāo)[17],喬木重要值(IV)=(相對(duì)多度+相對(duì)頻度+相對(duì)顯著度)/3。

    1.3.2 種群年齡結(jié)構(gòu)劃分

    根據(jù)江南油杉種群的胸徑范圍,采用“空間代替時(shí)間”的研究方法[18],用徑級(jí)代替年齡級(jí),劃分為9個(gè)齡級(jí)。Ⅰ級(jí)幼苗:DBH<5 cm,且<0.5 m; Ⅱ級(jí)幼樹(shù):DBH<5 cm,且>0.5 m;Ⅲ級(jí)小樹(shù):5 cm≤ DBH<15 cm;Ⅳ級(jí)中樹(shù):15 cm≤DBH<25 cm;Ⅴ級(jí)中樹(shù):25 cm≤DBH<35 cm;Ⅵ級(jí)大樹(shù):35 cm≤ DBH<45 cm;Ⅶ級(jí)大樹(shù):45 cm≤DBH<55 cm; Ⅷ級(jí)老樹(shù):55 cm≤DBH<65 cm;Ⅸ級(jí)老樹(shù):65 cm≤ DBH。以齡級(jí)為橫坐標(biāo),對(duì)應(yīng)株數(shù)為縱坐標(biāo),繪制江南油杉種群年齡結(jié)構(gòu)圖。

    1.3.3 總體關(guān)聯(lián)性檢驗(yàn)

    在獨(dú)立性假設(shè)條件下,v的期望值為1。若v= 1,表示群落物種間無(wú)關(guān)聯(lián);若v>1,表示正關(guān)聯(lián);若v<1,表示負(fù)關(guān)聯(lián)。采用統(tǒng)計(jì)量(=v×)來(lái)檢驗(yàn)v值偏離1的顯著程度。如20.95()≤≤20.05()的概率為90%,則種間關(guān)聯(lián)不顯著;<20.95()或>20.05(),種間關(guān)聯(lián)顯著。

    1.3.4 種間關(guān)聯(lián)性分析

    1.3.5 種間相關(guān)性分析

    1.3.6 生態(tài)種組劃分

    根據(jù)主要樹(shù)種在各樣方中的重要值,結(jié)合Pearson和Spearman相關(guān)性檢驗(yàn)結(jié)果,運(yùn)用PCA二維排序劃分生態(tài)種組[22],數(shù)據(jù)處理在Excel 2019和Origin 2022軟件中完成。

    2 結(jié)果和分析

    2.1 物種組成及重要值

    江南油杉群落喬木層共調(diào)查到植物31科50屬77種, 其中含2種以上的有15科,最多的是樟科(Lauraceae, 10種)和殼斗科(Fagaceae, 8種),其次為冬青科(Aquifoliaceae)、薔薇科(Rosaceae)和山茶科(Theaceae)各5種,漆樹(shù)科(Anacardiaceae)和山礬科(Symplocaceae)各4種,安息香科(Styracaceae)、杜鵑花科(Ericaceae)、金縷梅科(Hamamelidaceae)和五列木科(Pentaphylacaceae)各3種,杜英科(Elaeocar- paceae)、木蘭科(Magnoliaceae)、柿科(Ebenaceae)和松科(Pinaceae)各2種,其他為單科單屬單種。江南油杉的重要值最大(11.58),其次是杉木(11.05)和甜櫧(10.11),為該群落的共同優(yōu)勢(shì)種;馬尾松(6.56)、赤楊葉(5.46)和南酸棗(5.32)為亞優(yōu)勢(shì)種; 重要值大于1的還有虎皮楠(, 3.44)、薄葉山礬(, 2.40)和苦櫧(, 2.18)等15種(表1)。

    2.2 種群年齡結(jié)構(gòu)

    江南油杉種群年齡結(jié)構(gòu)見(jiàn)圖1,群落內(nèi)江南油杉共334株。第Ⅰ齡級(jí)幼苗288株,占總株數(shù)的86.23%, 林分密度為1 200 ind./hm2;第Ⅱ齡級(jí)幼樹(shù)僅17株,占總株數(shù)的5.09%,林分密度為78.8 ind./hm2;從第Ⅲ齡級(jí)后DBH≥5 cm的喬木29株,僅占總株數(shù)的8.68%,其中小樹(shù)(第Ⅲ齡級(jí)) 13株,中樹(shù)和大樹(shù)(第Ⅵ~Ⅷ齡級(jí))共6株,老樹(shù)(第Ⅷ~Ⅸ齡級(jí)) 7株,胸徑最大達(dá)159.5 cm??傮w來(lái)看,第Ⅰ到Ⅱ齡級(jí)株數(shù)驟減,幼苗死亡率極高,很難向幼樹(shù)階段發(fā)展;第Ⅲ齡級(jí)后株數(shù)減少趨勢(shì)明顯, 表明小樹(shù)難以通過(guò)環(huán)境的篩選,大多數(shù)個(gè)體在小樹(shù)階段就已淘汰。

    2.3 總體關(guān)聯(lián)性

    根據(jù)喬木層21種主要樹(shù)種在樣方中是否存在的二元矩陣,計(jì)算方差比率(V),得V=0.780<1, 表明21個(gè)種對(duì)間呈負(fù)關(guān)聯(lián)。其檢驗(yàn)統(tǒng)計(jì)量=4.683。根據(jù)自由度查2表[23],得20.95()=1.636,20.05()= 12.592,檢驗(yàn)統(tǒng)計(jì)量居于兩者之間,說(shuō)明21種主要樹(shù)種間表現(xiàn)為不顯著負(fù)關(guān)聯(lián)(>0.05),各物種間雖然有一定聯(lián)系,但尚未形成穩(wěn)定搭配,存在競(jìng)爭(zhēng)排斥關(guān)系,對(duì)環(huán)境資源的利用未達(dá)到最佳狀態(tài)。

    圖1 江南油杉種群年齡結(jié)構(gòu)。I: DBH<5 cm, H<0.5 m; II: DBH<5 cm, H>0.5 m; III: 5 cm≤DBH<15 cm; IV: 15 cm≤DBH<25 cm; V: 25 cm≤ DBH<35 cm; VI: 35 cm≤DBH<45 cm; VII: 45 cm≤DBH<55 cm; VIII: 55 cm≤DBH<65 cm; IX: 65cm≤DBH。

    表1 喬木層主要樹(shù)種的重要值(IV)

    2.4 主要樹(shù)種種間關(guān)聯(lián)性

    2檢驗(yàn)統(tǒng)計(jì)表明(表2, 圖2),在210個(gè)種對(duì)中,正相關(guān)的有110對(duì),占52.38%,其中僅有苦櫧-米櫧()為極顯著正相關(guān),占0.48%,其余均為不顯著正相關(guān)。負(fù)相關(guān)的有81對(duì),占38.57%,均為不顯著。不相關(guān)的有19對(duì),占9.05%,包括杉木-芬芳安息香()、杉木-黃丹木姜子()、甜櫧-黃丹木姜子和赤楊葉-苦櫧等。正負(fù)比1.36,顯著率為0.48%。2檢驗(yàn)結(jié)果表明江南油杉與其他種群呈不顯著相關(guān),說(shuō)明大多數(shù)物種為獨(dú)立分布,且主要樹(shù)種對(duì)環(huán)境的要求有較大差異。

    2.5 主要樹(shù)種的種間相關(guān)性

    Pearson檢驗(yàn)結(jié)果表明(表2, 圖3),正相關(guān)的種對(duì)有95對(duì),占45.24%,其中極顯著、顯著和不顯著正相關(guān)的種對(duì)分別為1、5和89對(duì)。負(fù)相關(guān)的種對(duì)有115對(duì), 占54.76%, 其中極顯著、顯著和不顯著負(fù)相關(guān)的種對(duì)分別為1、6和108對(duì)。正負(fù)比0.83,檢驗(yàn)顯著率為6.19%,說(shuō)明種間關(guān)聯(lián)極其松散,大多物種獨(dú)立分布,這與2檢驗(yàn)結(jié)果相似。

    Spearman檢驗(yàn)結(jié)果表明(表2, 圖4),喬木層正相關(guān)的種對(duì)有96對(duì),占45.72%,其中極顯著、顯著和不顯著正相關(guān)的種對(duì)分別為4、2和90對(duì)。負(fù)相關(guān)的種對(duì)有112對(duì), 占53.33%,其中極顯著、顯著和不顯著負(fù)相關(guān)的種對(duì)分別為7、2和103對(duì)。不相關(guān)的種對(duì)有2對(duì),占0.95%。正負(fù)比0.86,檢驗(yàn)顯著率為7.14%。負(fù)相關(guān)種對(duì)數(shù)大于正相關(guān)種對(duì)數(shù),表明物種間彼此存在競(jìng)爭(zhēng),生態(tài)位互相分離。這與2檢驗(yàn)和Pearson相關(guān)系數(shù)檢驗(yàn)的結(jié)果一致。

    表2 喬木層主要樹(shù)種對(duì)間的x2檢驗(yàn)、Pearson相關(guān)檢驗(yàn)和Spearman秩相關(guān)檢驗(yàn)

    括號(hào)中的數(shù)值為百分率。

    The values in brackets are percentages.

    圖2 喬木層主要樹(shù)種x2檢驗(yàn)半矩陣圖。1~21見(jiàn)表1。下同

    圖3 喬木層主要樹(shù)種Pearson相關(guān)性熱力圖。紅色: 正相關(guān); 藍(lán)色: 負(fù)相關(guān); 顏色越深, 圓形越大, 表示相關(guān)性越強(qiáng);*: P<0.05;**: P<0.01。下同

    圖4 喬木層主要樹(shù)種Spearman相關(guān)性熱力圖

    2.6 生態(tài)種組的劃分

    喬木層主要樹(shù)種大致可劃分為4個(gè)生態(tài)種組(圖5)。第Ⅰ生態(tài)種組有虎皮楠、馬尾松、秀麗錐()、杉木、苦櫧、紅皮糙果茶和米櫧,該組植物除馬尾松外普遍較耐蔭,喜暖濕氣候,適應(yīng)性較強(qiáng),杉木多與其他樹(shù)種形成針闊混交林。第Ⅱ生態(tài)種組有光葉山礬、芬芳安息香、檵木()、南酸棗和赤楊葉,該組植物個(gè)體數(shù)量最多,生態(tài)幅較廣,生境趨于一致,多生于微酸性土壤,且種間絕大多數(shù)呈正相關(guān)。第Ⅲ生態(tài)種組包括江南油杉、黃丹木姜子、尖葉川楊桐(var.)和紅楠(),該組植物喜陽(yáng)又略耐蔭,其中江南油杉重要值最大。第Ⅵ生態(tài)種組有浙江潤(rùn)楠()、柯()、薄葉山礬、樹(shù)參()和甜櫧,該組植物較耐蔭,耐寒,萌發(fā)力較強(qiáng),在群落中出現(xiàn)的樣方數(shù)較少,與其他生態(tài)種組樹(shù)種多呈負(fù)相關(guān)。

    圖5 喬木層主要樹(shù)種PCA二維排序圖

    3 結(jié)論和討論

    3.1 不同檢驗(yàn)方法的比較

    比較2檢驗(yàn)、Pearson相關(guān)系數(shù)檢驗(yàn)和Spearman秩相關(guān)系數(shù)檢驗(yàn)結(jié)果表明,3種檢驗(yàn)方法結(jié)果雖具有一致性,但也存在一定差異,本研究主要反映在以下2個(gè)方面:(1)2檢驗(yàn)中,210個(gè)種對(duì)間正相關(guān)種對(duì)數(shù)大于負(fù)相關(guān)種對(duì)數(shù),但Pearson相關(guān)系數(shù)和Spearman秩相關(guān)系數(shù)檢驗(yàn)則表現(xiàn)出相反結(jié)果;(2)2檢驗(yàn)顯著率僅為0.48%,Pearson相關(guān)系數(shù)檢驗(yàn)顯著率為6.19%,Spearman秩相關(guān)系數(shù)檢驗(yàn)顯著率為7.14%,后兩者的顯著率略高于2檢驗(yàn)。造成差異的原因在于2檢驗(yàn)是基于主要樹(shù)種在該群落樣方中是否存在的0、1二元數(shù)據(jù)作為判斷標(biāo)準(zhǔn),未考慮多度、頻度和重要值等信息,僅能判斷2物種間是否存在關(guān)聯(lián)性,無(wú)法判斷關(guān)聯(lián)程度的大小,屬于定性檢驗(yàn);而Pearson相關(guān)系數(shù)檢驗(yàn)和Spearman秩相關(guān)系數(shù)檢驗(yàn)則是以主要樹(shù)種在群落樣方中的重要值為依據(jù),考慮到了物種在樣方中的功能與地位,屬于定量檢驗(yàn),能夠較準(zhǔn)確地反映出2檢驗(yàn)中不顯著種對(duì)的關(guān)聯(lián)程度,如杉木-紅皮糙果茶、苦櫧-紅楠、樹(shù)參-芬芳安息香等。此外,Pearson相關(guān)系數(shù)檢驗(yàn)和Spearman秩相關(guān)系數(shù)檢驗(yàn)也不盡相同, 前者屬于參數(shù)檢驗(yàn),要求物種必須符合正態(tài)分布, 存在一定局限性;而后者則是一種非參數(shù)檢驗(yàn)方法,對(duì)物種的分布形式?jīng)]有要求,應(yīng)用起來(lái)相對(duì)靈活方便,靈敏度高,如甜櫧-檵木、甜櫧-樹(shù)參和薄葉山礬-光葉山礬等種對(duì),在Pearson相關(guān)系數(shù)檢驗(yàn)中不顯著, 但Spearman秩相關(guān)系數(shù)檢驗(yàn)中均表現(xiàn)為顯著相關(guān)。涂洪潤(rùn)等[24]對(duì)桂林巖溶石山青岡()群落的研究也證實(shí)了這一點(diǎn)。由此可見(jiàn),應(yīng)將Spearman秩相關(guān)系數(shù)檢驗(yàn)和Pearson相關(guān)系數(shù)檢驗(yàn)與2檢驗(yàn)充分結(jié)合,才能更全面、準(zhǔn)確反映出江南油杉群落主要樹(shù)種之間的種間關(guān)聯(lián)性。

    3.2 種間關(guān)聯(lián)與群落演替的關(guān)系

    種間關(guān)聯(lián)性與群落演替密切相關(guān),種群間的總體關(guān)聯(lián)性反映群落演替的階段及其穩(wěn)定性[25]。相關(guān)研究認(rèn)為,隨著演替的進(jìn)行,種對(duì)間的正負(fù)相關(guān)比例也伴隨著發(fā)生變化,負(fù)相關(guān)種對(duì)呈下降趨勢(shì),正相關(guān)種對(duì)則可能大幅度上升,物種間逐漸趨于獨(dú)立分布,多物種在群落中共存,群落穩(wěn)定性逐漸變強(qiáng),因此越趨向頂級(jí)群落演替[26–27]。本研究中,江南油杉群落主要樹(shù)種的總體關(guān)聯(lián)性表現(xiàn)為不顯著負(fù)關(guān)聯(lián),負(fù)相關(guān)種對(duì)數(shù)大于正相關(guān)種對(duì)數(shù),表明種間關(guān)聯(lián)性較弱,物種對(duì)生境需求具有較大差異,物種組成尚未達(dá)到穩(wěn)定和完善,這與山西五鹿山自然保護(hù)區(qū)暴馬丁香()群落木本植物種間聯(lián)結(jié)性結(jié)果相似[28]。主要有以下幾方面原因:一是各樹(shù)種在生長(zhǎng)發(fā)育和群落演替過(guò)程中,受到各方面因素的影響,包括土壤、水分、動(dòng)物和微生物等[29]; 二是喬木之間也存在激烈競(jìng)爭(zhēng),如江南油杉、馬尾松、赤楊葉和南酸棗等樹(shù)高多在15 m以上,而下木層樹(shù)高相對(duì)較低,樹(shù)種垂直分層明顯,且多以陽(yáng)性樹(shù)種為主,對(duì)環(huán)境資源尤其是光照競(jìng)爭(zhēng)激烈,各樹(shù)種在群落中占據(jù)不同的生態(tài)位。由此可見(jiàn),江南油杉群落處于演替中期。

    江南油杉與杉木和甜櫧無(wú)顯著相關(guān)性,且與其他絕大多數(shù)樹(shù)種相關(guān)性都不強(qiáng),與群落中其他樹(shù)種尚未達(dá)到穩(wěn)定共存格局,這是因?yàn)閱棠緦又械慕^對(duì)優(yōu)勢(shì)種能夠充分利用資源和環(huán)境,與其他樹(shù)種不會(huì)產(chǎn)生強(qiáng)烈的競(jìng)爭(zhēng),有其最適宜的生態(tài)位[30]。本研究絕大多數(shù)種對(duì)未達(dá)到顯著水平,各物種間互相獨(dú)立。但一些物種之間也存在顯著相關(guān),如杉木、紅皮糙果茶和米櫧彼此間呈極顯著正相關(guān)或顯著相關(guān),原因可能是3樹(shù)種都較耐蔭,喜暖濕氣候和深厚酸性土壤,生態(tài)學(xué)特性具有相似性,能夠適應(yīng)相同生境,在群落中共存。馬尾松-黃丹木姜子則表現(xiàn)為顯著負(fù)相關(guān),兩者存在競(jìng)爭(zhēng)關(guān)系。實(shí)地調(diào)查發(fā)現(xiàn),群落中馬尾松多為大型喬木,樹(shù)形高大, 位于第一喬木層,且多集中在Ⅴ和Ⅵ號(hào)樣地中;而黃丹木姜子多為中小型喬木,在Ⅴ和Ⅵ號(hào)樣地中均未出現(xiàn)。這可能與樹(shù)種本身特性的差異有關(guān),馬尾松為針葉樹(shù), 其針葉細(xì)小、稀疏,頂端優(yōu)勢(shì)明顯, 因此高生長(zhǎng)顯著; 黃丹木姜子則為闊葉樹(shù),葉片長(zhǎng)圓狀披針形, 生長(zhǎng)較緩慢,相比之下馬尾松更適宜該生境。此外,南酸棗-樹(shù)參、赤楊葉-虎皮楠也表現(xiàn)為極顯著負(fù)相關(guān)或顯著負(fù)相關(guān),原因可能是南酸棗與赤楊葉為落葉樹(shù)種, 樹(shù)參與虎皮楠為常綠樹(shù)種, 由于兩者間的生態(tài)學(xué)特性不同,對(duì)環(huán)境具有不同的生態(tài)適應(yīng)性或相互分離的生態(tài)位所致,也可能與他感作用有關(guān)[31]。其他影響種間關(guān)系的因素很多,如群落發(fā)展演替的不同階段,不同樹(shù)種個(gè)體的發(fā)育時(shí)間、功能性狀和譜系多樣性的差異等,今后可對(duì)這些因素進(jìn)行深入探討, 進(jìn)一步探究群落中物種種間關(guān)系的形成機(jī)制。

    3.3 主要樹(shù)種生態(tài)種組

    劃分生態(tài)種組有利于探討不同物種間的彼此協(xié)同適應(yīng)機(jī)制,準(zhǔn)確表達(dá)物種對(duì)環(huán)境資源和生態(tài)因子的適應(yīng)策略[32]。相同生態(tài)種組內(nèi)物種多呈正相關(guān),如杉木、紅皮糙果茶和米櫧,它們處于第Ⅰ生態(tài)種組,彼此間關(guān)聯(lián)性較強(qiáng),對(duì)資源利用方式、環(huán)境適應(yīng)能力以及植物群落結(jié)構(gòu)影響等方面具有相似性,在群落中相互促進(jìn)和依賴[24]。相反,不同生態(tài)種組間的物種則有不同的生物學(xué)習(xí)性,這些物種長(zhǎng)期適應(yīng)環(huán)境并占據(jù)不同資源空間,馬尾松-黃丹木姜子、赤楊葉-虎皮楠、南酸棗-樹(shù)參、虎皮楠-芬芳安息香等表現(xiàn)為極顯著或顯著負(fù)相關(guān),因物種生物學(xué)特性的差異和相互重疊的生態(tài)位,導(dǎo)致物種間相互排斥,在群落中表現(xiàn)為競(jìng)爭(zhēng)關(guān)系,分布在不同的生態(tài)種組。這與劉潤(rùn)紅等[33]對(duì)漓江河岸帶楓楊()群落種間關(guān)聯(lián)的研究結(jié)果一致。此外,在不同生態(tài)種組間也發(fā)現(xiàn)呈正相關(guān)的種對(duì),相同生態(tài)種組間也存在部分負(fù)相關(guān)種對(duì),說(shuō)明重要值和物種生態(tài)學(xué)特性不是影響生態(tài)種組劃分的決定因素。由于某些物種頻度較小,或因干擾導(dǎo)致在某些樣方中已淘汰,對(duì)生態(tài)種組的劃分將產(chǎn)生一定影響。此外,群落面積大小、環(huán)境變化、種群變異和演替進(jìn)程等都可能是造成差異的因素[34],只有充分考慮這些因素,與植物生態(tài)習(xí)性結(jié)合分析, 才能更加準(zhǔn)確地劃分生態(tài)種組。

    本研究中喬木層植物共有792株,而江南油杉只有29株,僅占3.67%。在母樹(shù)周圍有大量更新幼苗,但幼樹(shù)極少,說(shuō)明該群落中江南油杉幼苗死亡率高,長(zhǎng)成幼樹(shù)的可能性較低。江南油杉屬陽(yáng)性樹(shù)種,幼苗多集中生于箬竹叢下,常年處于庇蔭環(huán)境,難以獲得充足的光照,導(dǎo)致幼苗過(guò)早死亡,天然更新困難, 這是江南油杉瀕危的主要原因[15]。其次, 江南油杉種子敗育率高、質(zhì)量低[35],也可能是導(dǎo)致其瀕危的原因之一。江南油杉幼苗對(duì)環(huán)境的敏感度高,存活率低,制約了該種群的自然更新與發(fā)展, 如忽視對(duì)幼苗的保護(hù),其種群在未來(lái)一段時(shí)間內(nèi)可能被杉木這類進(jìn)展種群取代。杉木萌蘗更新能力強(qiáng), 生長(zhǎng)快,且幼、中齡林個(gè)體數(shù)量多,將逐步取代江南油杉在群落中的地位,與甜櫧、馬尾松和赤楊葉等形成針闊混交林。

    江南油杉群落喬木層主要樹(shù)種的種間關(guān)聯(lián)性較為松散,大多數(shù)物種為獨(dú)立分布,幼苗存在更新障礙,且幼樹(shù)個(gè)體競(jìng)爭(zhēng)壓力較大,種群規(guī)模呈逐漸縮小的趨勢(shì)。為更有效地保護(hù)該地江南油杉天然種群,建議人工促進(jìn)江南油杉幼苗更新,對(duì)群落內(nèi)箬竹等蓋度高的小灌木進(jìn)行適當(dāng)疏伐,提高透光率, 改善幼苗生長(zhǎng)環(huán)境;同時(shí)應(yīng)盡快建立江南油杉種質(zhì)資源庫(kù),確保種質(zhì)資源得到有效保護(hù);加強(qiáng)人工育苗工作,促進(jìn)江南油杉種群的繁育與壯大。

    致謝 在野外調(diào)查過(guò)程中,得到福建君子峰國(guó)家級(jí)自然保護(hù)區(qū)管理局的大力支持,在此深表謝意!

    [1] ZHANG J, CAIWEN D J, SUONAN C R, et al. The correlation and niche of survival community ofin Maduo County, Qinghai Province [J]. Pratacult Sci, 2019, 36(11): 2752–2765. [張靜, 才文代吉, 索南才仁, 等. 青?,敹嗨罅_以禮草生存群落種間關(guān)聯(lián)及生態(tài)位 [J]. 草業(yè)科學(xué), 2019, 36(11): 2752–2765. doi: 10. 11829/j.issn.1001-0629.2019-0310.]

    [2] WU P, PENG X Q, YANG S R, et al. Spatial distribution patterns and correlation ofpopulation in coastal wetlands of Shandong, China [J]. Chin J Plant Ecol, 2019, 43(9): 817–824. [吳盼, 彭希強(qiáng), 楊樹(shù)仁, 等. 山東省濱海濕地檉柳種群的空間分布格局及其關(guān)聯(lián)性 [J]. 植物生態(tài)學(xué)報(bào), 2019, 43(9): 817–824. doi: 10.17521/ cjpe.2018.0186.]

    [3] SU S J, LIU J F, HE Z S, et al. Ecological species groups and inter- specific association of dominant tree species in Daiyun Mountain National Nature Reserve [J]. J Mount Sci, 2015, 12(3): 637–646. doi: 10.1007/s11629-013-2935-7.

    [4] LUO Y H, WU J S, WANG X Y, et al. Understanding ecological groups under landscape fragmentation based on network theory [J]. Landscape Urban Plan, 2021, 210: 104066. doi: 10.1016/j.landurbplan.2021.104066.

    [5] GAO W, HUANG Y R, LIN J L, et al. Interspecific correlations among dominant populations of natural forest of endangered species[J]. Sci Silv Sin, 2021, 57(10): 1–14. [高偉, 黃雍容, 林建麗, 等. 瀕危樹(shù)種閩樺天然林優(yōu)勢(shì)種群種間相關(guān)性 [J]. 林業(yè)科學(xué), 2021, 57(10): 1–14. doi: 10.11707/j.1001-7488.20211001.]

    [6] WEN X H, WANG Q B, PAN H, et al. Interspecific associations of the main tree populations of thecommunity in Tianbaoyan [J]. J Forest Environ, 2022, 42(1): 1–10. [溫鑫鴻, 王其炳, 潘輝, 等. 天寶巖柳杉群落主要喬木種群的種間聯(lián)結(jié)性 [J]. 森林與環(huán)境學(xué)報(bào), 2022, 42(1): 1–10. doi: 10.13324/j.cnki.jfcf.2022.01.001.]

    [7] YE Q P, ZHANG W H, YU S C, et al. Interspecific association of the main tree populations of thecommunity in the Qiaoshan forest area [J]. Acta Ecol Sin, 2018, 38(9): 3165–3174. [葉權(quán)平, 張文輝, 于世川, 等. 橋山林區(qū)麻櫟群落主要喬木種群的種間聯(lián)結(jié)性 [J]. 生態(tài)學(xué)報(bào), 2018, 38(9): 3165–3174. doi: 10.5846/stxb20 1706131071.]

    [8] National Forestry and Grassland Administration, Ministry of Agri- culture and Rural Affairs of the People’s Republic China. Perspectives from the list of national key protected wild plant [EB/OL]. (2021- 09-08)[2022-05-13]. [國(guó)家林業(yè)和草原局, 農(nóng)業(yè)農(nóng)村部. 國(guó)家重點(diǎn)保護(hù)野生植物名錄 [EB/OL]. (2021-09-08)[2022-05-13]. http://www. forestry.gov.cn/main/5461/20210908/162515850572900.html.]

    [9] YANG M M, HE W G, CHEN W R, et al. Phenotypic traits diversity analysis of seeds of candidate superior trees of[J]. J Fujian For Sci Technol, 2020, 47(4): 18–21. [楊淼淼, 何文廣, 陳文榮,等. 江南油杉優(yōu)樹(shù)種子表型性狀的多樣性分析 [J]. 福建林業(yè)科技, 2020, 47(4): 18–21. doi: 10.13428/j.cnki.fjlk.2020.04.004.]

    [10] LIU X S, XIAO Y F, WANG Y, et al. Anatomical structures of vegetative organs of(Murr.) Carr. var.(Flous) Silba and its ecological adaptability [J]. Plant Sci J, 2020, 38(1): 39–46. [劉雄盛, 肖玉菲, 王勇, 等. 江南油杉營(yíng)養(yǎng)器官的解剖結(jié)構(gòu)及其生態(tài)適應(yīng)性 [J]. 植物科學(xué)學(xué)報(bào), 2020, 38(1): 39–46. doi: 10. 11913/PSJ.2095-0837.2020.10039.]

    [11] LIU F, ZHOU L T, JIANG Y, et al. Physiological response from different provenances ofseedlings to drought stress [J]. J CS Univ For Technol, 2018, 38(11): 35–45. [劉菲, 周隆騰, 蔣燚, 等. 不同種源江南油杉幼苗對(duì)干旱脅迫的生理響應(yīng) [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2018, 38(11): 35–45. doi: 10.14067/j.cnki.1673-923x. 2018.11.006.]

    [12] ZHANG P, PANG S J, LIU S L, et al. Effects of slow release fertilizer on growth ofseedlings cultured in container [J]. J NW Agric For Univ (Nat Sci), 2021, 49(9): 92–98. [張培, 龐圣江, 劉士玲, 等. 緩釋肥對(duì)江南油杉容器苗生長(zhǎng)的影響 [J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版), 2021, 49(9): 92–98. doi: 10.13207/j.cnki. jnwafu.2021.09.011.]

    [13] HE Y H, JIANG Y, HUANG R L, et al. Spatial structure of fine root biomass and its response to soil moisture ofFlous [J]. J CS Univ For Technol, 2021, 41(1): 151–158. [何應(yīng)會(huì), 蔣燚, 黃榮林, 等. 江南油杉細(xì)根生物量空間分布及其對(duì)土壤水分的響應(yīng) [J]. 中南林業(yè)科技大學(xué)學(xué)報(bào), 2021, 41(1): 151–158. doi: 10.14067/j. cnki.1673-923x.2021.01.016.]

    [14] LIU X S, JIANG Y, WEI S X, et al. Establishment of the ISSR-PCR reaction system and primers screening invar.[J]. Guangxi For Sci, 2017, 46(3): 266–270. [劉雄盛, 蔣燚, 韋鑠星, 等. 江南油杉ISSR-PCR反應(yīng)體系建立與引物篩選 [J]. 廣西林業(yè)科學(xué), 2017, 46(3): 266–270. doi: 10.19692/j.cnki.gfs.2017.03. 007.]

    [15] LI Q, HUANG R L, LIU X S, et al. Community structure charac- teristics ofvar.natural forest in Nandan County, Guangxi [J]. Guangxi For Sci, 2019, 48(2): 183–189. [李強(qiáng), 黃榮林, 劉雄盛, 等. 廣西南丹縣江南油杉天然林群落結(jié)構(gòu)特征 [J]. 廣西林業(yè)科學(xué), 2019, 48(2): 183–189. doi: 10.19692/j.cnki.gfs. 2019.02.009.]

    [16] LUO J W. Community structure and plant diversity ofvar.in Guangze County of Fujian Province [J]. Prot For Sci Technol, 2016(4): 11–16. [羅金旺. 福建光澤江南油杉群落結(jié)構(gòu)與植物多樣性 [J]. 防護(hù)林科技, 2016(4): 11–16. doi: 10.13601/j. issn.1005-5215.2016.04.003.]

    [17] SUN J J, SHEN A H, HUANG Y J, et al. Quantitative classification and ordination ofhabitat in Zhejiang Province [J]. J Nanjing For Univ (Nat Sci), 2019, 43(4): 85–93. [孫杰杰, 沈愛(ài)華, 黃玉潔, 等. 浙江省大葉櫸樹(shù)生境地群落數(shù)量分類與排序 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2019, 43(4): 85–93. doi: 10.3969/j. issn.1000-2006.201809027.]

    [18] WANG C H, CHEN X, WANG B Z, et al. Population structure and dynamics ofvar.in the Gaowangjie National Nature Reserve of Hunan Province [J]. J Nanjing For Univ (Nat Sci), 2022, 46(3): 57–64. [王春暉, 陳昕, 王本忠, 等. 湖南高望界國(guó)家級(jí)自然保護(hù)區(qū)雪峰山梭羅種群結(jié)構(gòu)與動(dòng)態(tài)特征 [J]. 南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2022, 46(3): 57–64. doi: 10.12302/j.issn. 1000-2006.202105012.]

    [19] SCHLUTER D. A variance test for detecting species associations, with some example applications [J]. Ecology, 1984, 65(3): 998–1005. doi: 10.2307/1938071.

    [20] WANG B S, PENG S L. Studies on the measuring techniques of interspecific association of lower-subtropical evergreen-broadleaved forests. I. The exploration and the revision on the measuring formulas of interspecific association [J]. Acta Phytoecol Geobot Sin, 1985, 9(4): 274–285. [王伯蓀, 彭少麟. 南亞熱帶常綠闊葉林種間聯(lián)結(jié)測(cè)定技術(shù)研究, Ⅰ. 種間聯(lián)結(jié)測(cè)式的探討與修正 [J]. 植物生態(tài)學(xué)與地植物學(xué)叢刊, 1985, 9(4): 274–285.]

    [21] ZHANG J T. Quantitative Ecology [M]. 2nd ed. Beijing: Science Press, 2011: 105–109. [張金屯. 數(shù)量生態(tài)學(xué) [M]. 第2版. 北京: 科學(xué)出版社, 2011: 105–109.]

    [22] JIN S, WU K S. Niche and interspecific association of dominant species in herb layer of burnedforest in the southern Taihang Mountain of northern China [J]. J Beijing For Univ, 2021, 43(4): 35–46. [金山, 武楷帥. 太行山南段油松林火燒跡地優(yōu)勢(shì)草本生態(tài)位及種間關(guān)系 [J]. 北京林業(yè)大學(xué)學(xué)報(bào), 2021, 43(4): 35–46. doi: 10.12171/j.1000-1522.20210044.]

    [23] DU R Q. Biostatistics [M]. Beijing: Higher Education Press, 1999: 277–278. [杜榮騫. 生物統(tǒng)計(jì)學(xué) [M]. 北京: 高等教育出版社, 1999: 277–278.]

    [24] TU H R, LI J F, YANG L T, et al. Interspecific associations of the main tree populations of thecommunity in Karst hills of Guilin, southwest China [J]. Chin J Appl Ecol, 2019, 30(1): 67–76. [涂洪潤(rùn), 李嬌鳳, 楊麗婷, 等. 桂林巖溶石山青岡群落主要喬木樹(shù)種的種間關(guān)聯(lián) [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2019, 30(1): 67–76. doi: 10.13287/j.1001-9332.201901.018.]

    [25] XU M H, LIU M, ZHAI D T, et al. A review of contents and methods used to analyze various aspects of plant interspecific associations [J]. Acta Ecol Sin, 2016, 36(24): 8224–8233. [徐滿厚, 劉敏, 翟大彤, 等. 植物種間聯(lián)結(jié)研究?jī)?nèi)容與方法評(píng)述 [J]. 生態(tài)學(xué)報(bào), 2016, 36(24): 8224–8233. doi: 10.5846/stxb201505311092.]

    [26] GU L, GONG Z W, LI W Z. Niches and interspecific associations of dominant populations in three changed stages of natural secondary forests on Loess Plateau, P. R. China [J]. Sci Rep, 2017, 7(1): 6604. doi: 10.1038/s41598-017-06689-9.

    [27] LIU L T, WANG X D, WEN Q, et al. Interspecific associations of plant populations in rare earth mining wasteland in southern China [J]. Int Biodeteriorat Biodegradat, 2017, 118: 82–88. doi: 10.1016/j.ibiod.2017. 01.011.

    [28] YANG Z J, ZHANG Q D, LI H, et al. Inter-specific correlations among dominant populations of woodycommunities in Wulu Mountain Nature Reserve, Shanxi Province of China [J]. Plant Sci J, 2013, 31(1): 42–48. [楊兆靜, 張欽弟, 李豪, 等. 山西五鹿山自然保護(hù)區(qū)暴馬丁香群落木本植物種間聯(lián)結(jié)性分析 [J]. 植物科學(xué)學(xué)報(bào), 2013, 31(1): 42–48. doi: 10.3724/SP.J.1142.2013.10042.]

    [29] JIN J Y, QIN W G, TAN W N, et al. Interspecific association among dominant populations in endangered plantcommunities [J]. J W China For Sci, 2013, 42(3): 86–94. [金俊彥, 覃文更, 譚衛(wèi)寧, 等. 瀕危植物單性木蘭群落主要種群種間聯(lián)結(jié)性研究 [J]. 西部林業(yè)科學(xué), 2013, 42(3): 86–94. doi: 10.16473/j.cnki.xbly kx1972.2013.03.001.]

    [30] JIAN M F, LIU Q, ZHU D, et al. Inter-specific correlations among dominant populations of tree layer species in evergreen broad-leaved forest in Jiulianshan Mountain of subtropical China [J]. Chin J Plant Ecol, 2009, 33(4): 672–680. [簡(jiǎn)敏菲, 劉琪, 朱篤, 等. 九連山常綠闊葉林喬木優(yōu)勢(shì)種群的種間關(guān)聯(lián)性分析 [J]. 植物生態(tài)學(xué)報(bào), 2009, 33(4): 672–680. doi: 10.3773/j.issn.1005-264x.2009.04.005.]

    [31] ZHU L J, SU Z X, WANG W W, et al. Numerical analysis on inter- specific relationships incommunity in Sanjiang Wolong Nature Reserve [J]. Chin J Ecol, 2005, 24(10): 1167–1171. [朱利君, 蘇智先, 王偉偉, 等. 臥龍自然保護(hù)區(qū)三江珙桐群落種間關(guān)系的數(shù)量分析 [J]. 生態(tài)學(xué)雜志, 2005, 24(10): 1167–1171.]

    [32] CHEN Q, CHEN J, ZHONG J J, et al. Interspecific association and functional group classification of the dominant populations in shrub layer in secondary forest ofin Qinling Mountain, China [J]. Chin J Appl Ecol, 2018, 29(6): 1736–1744. [陳倩, 陳杰, 鐘嬌嬌, 等. 秦嶺山地油松天然次生林灌木層主要種群種間聯(lián)結(jié)性與功能群劃分 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2018, 29(6): 1736–1744. doi: 10. 13287/j.1001-9332.201806.004.]

    [33] LIU R H, JIANG Y, CHANG B, et al. Interspecific associations and correlations among the main woody plants in acommunity in a riparian zone of Lijiang River, Guilin, southwest China [J]. Acta Ecol Sin, 2018, 38(19): 6881–6893. [劉潤(rùn)紅, 姜勇, 常斌, 等. 漓江河岸帶楓楊群落主要木本植物種間聯(lián)結(jié)與相關(guān)分析 [J]. 生態(tài)學(xué)報(bào), 2018, 38(19): 6881–6893. doi: 10.5846/stxb201802270396.]

    [34] DING W H, LI X Z, JIANG J Y, et al. Numerical analysis of inter- specific relationships in the estuary salt marsh plant community of southern Chongming Dongtan, Shanghai [J]. Chin J Appl Ecol, 2016, 27(5): 1417–1426. [丁文慧, 李秀珍, 姜俊彥, 等. 崇明東灘南部河口鹽沼植物群落種間關(guān)系的數(shù)量分析 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27 (5): 1417–1426. doi: 10.13287/j.1001-9332.201605.006.]

    [35] JIANG Y, LIU F, LIU X S, et al. Preservation and evaluation of germplasm resources of precious local tree speciesvar.and construction of strong seedling breeding system [J]. Guangxi For Sci, 2022, 51(1): 1–9. [蔣燚, 劉菲, 劉雄盛, 等. 珍貴鄉(xiāng)土樹(shù)種江南油杉種質(zhì)資源保存評(píng)價(jià)及壯苗繁育體系構(gòu)建技術(shù) [J]. 廣西林業(yè)科學(xué), 2022, 51(1): 1–9. doi: 10.19692/j.issn.1006-1126.2022 0101.]

    Interspecific Association Analysis of Dominant Tree Species of Rare and Endangered Plant invarCommunity

    CHEN Xuhui1a,1b, YE Baojian2, PAN Biaozhi3, LIN Bihua1a,1b, LUO Minxian1a,1b, XIAO Lifang1a,1b, LIU Bao1a, HE Zongming1a, ZHENG Shiqun1a,1b*

    (1a. College of forestry;1b. Cross-strait Nature Reserve Research Center, Key Laboratory of Fujian Universities for Ecology and Resource Statistics, Fujian Agriculture and Forestry University, Fuzhou 350002, China; 2. College of Architecture and Urban Planning, Fujian University of Technology, Fuzhou 350118, China; 3. Forest Investigation and Planning Institute of Fujian Province, Fuzhou 350003, China)

    The correlations of dominant tree species ofvarcommunity in Junzifeng Nature Reserve, Fujian Province were analyzed by using 2×2 contingency table, variance ratio (v),2-test, Pearson correlation coefficient test and Spearman rank correlation coefficient test, and the ecological species groups were divided based on principal component analysis (PCA). The results showed that there was not significant negative correlations existed among the 21 dominant tree species, andvar.had not significant association with other species. Among 210 species pairs, 110 species pairs were positively correlated, 81 species pairs were negatively correlated, and 19 species pairs were unrelated, with a positive-negative ratio of 1.36. Most species pairs were not significantly correlated each other, and the significantly correlated rate was only 0.48%. The results of Pearson and Spearman rank correlation coefficient test were similar. The number of negative correlation species pairs was greater than that of positive correlation species pairs, and most species pairs were not significantly, and the significantly rates were 6.19% and 7.14%, respectively. Compared with2-test, the other two test were more sensitive. The dominant tree species could be divided into four ecological species groups by PCA two-dimensional sequence. The ecological habits and resource requirements of species in the same group were similar, most of them were positively correlated, and interde- pendence coexisted. There was a negative correlation among species in different ecological species groups, and species pairs competed with each other. The interspecific association of dominant tree species ofvar.community was relatively non-significant, and most species were independently distributed. The seedlings were difficult to renew, saplings face greater competitive pressure, and population size is gradually shrinking trend. In order to promote the normal growth and regeneration ofvar.population, the seedlings should be protected.

    var.; Rare and endangered plant; Interspecific association; Ecological species group

    10.11926/jtsb.4682

    2022-05-30

    2022-09-15

    福建君子峰自然保護(hù)區(qū)維管植物資源調(diào)查項(xiàng)目(KH200237A);福建省特色林木種質(zhì)資源評(píng)價(jià)與監(jiān)測(cè)信息應(yīng)用項(xiàng)目(2021FKJ17)資助

    This work was supported by the Project for Vascular Plant Resources Survey in Junzifeng Nature Reserve of Fujian (Grant No. KH200237A), and the Project for Evaluation and Monitoring Information Application of Characteristic Forest Germplasm Resources in Fujian (Grant No. 2021FKJ17).

    陳緒輝(1996年生),男,碩士研究生,主要從事珍稀植物保護(hù)與利用研究。E-mail: fafucxh123@126.com

    E-mail: fjzsq@126.com

    猜你喜歡
    油杉種間喬木
    干旱條件對(duì)鬼針草和醉魚(yú)草種間相互作用及生長(zhǎng)的影響
    植物研究(2023年5期)2023-09-09 08:01:22
    油杉屬的地理分布及其分類學(xué)研究*
    三峽庫(kù)區(qū)支流花溪河浮游植物種間關(guān)聯(lián)及影響因子分析
    ◆ 園林苗木
    27年生江南油杉人工林天然更新及其幼苗生長(zhǎng)特征
    廣西江南油杉天然林種群分布特征
    江南油杉優(yōu)樹(shù)子代家系苗期遺傳變異與早期選擇
    江蘇省宜興市茶園秋季雜草種間生態(tài)關(guān)系及群落分類
    沈陽(yáng)市政材料價(jià)格補(bǔ)充信息
    營(yíng)養(yǎng)鹽輸入方式對(duì)3種微藻生長(zhǎng)及種間競(jìng)爭(zhēng)的影響
    午夜精品国产一区二区电影 | 免费观看精品视频网站| 99久久九九国产精品国产免费| 人妻夜夜爽99麻豆av| 亚洲欧美精品专区久久| 精品久久久久久久久av| 亚洲最大成人手机在线| 国产 一区 欧美 日韩| 人妻制服诱惑在线中文字幕| 婷婷色综合大香蕉| 亚洲欧美日韩卡通动漫| 亚洲欧美精品综合久久99| 国产成人a区在线观看| 亚洲经典国产精华液单| 国产一区二区三区av在线| 国产欧美另类精品又又久久亚洲欧美| 美女cb高潮喷水在线观看| 蜜桃久久精品国产亚洲av| 久久亚洲国产成人精品v| 成人亚洲欧美一区二区av| 97人妻精品一区二区三区麻豆| 一级毛片电影观看 | 国产精品一区二区三区四区免费观看| 精品久久久久久久久亚洲| 国产精品爽爽va在线观看网站| 欧美色视频一区免费| 国产大屁股一区二区在线视频| 亚洲国产精品合色在线| 午夜a级毛片| 精品免费久久久久久久清纯| 在线免费十八禁| 一级爰片在线观看| 国产一级毛片七仙女欲春2| 久久婷婷人人爽人人干人人爱| 在线观看美女被高潮喷水网站| 免费黄色在线免费观看| 一级毛片aaaaaa免费看小| 日本免费在线观看一区| 国产高清视频在线观看网站| 国产片特级美女逼逼视频| 久久久a久久爽久久v久久| 免费无遮挡裸体视频| av在线观看视频网站免费| 国产黄a三级三级三级人| 毛片女人毛片| 中文字幕久久专区| 午夜日本视频在线| 免费大片18禁| 国产精品久久久久久久久免| 在线播放无遮挡| 22中文网久久字幕| 噜噜噜噜噜久久久久久91| 国产美女午夜福利| 秋霞伦理黄片| 久久久亚洲精品成人影院| 国产男人的电影天堂91| 国产高清有码在线观看视频| 九九爱精品视频在线观看| 黄片无遮挡物在线观看| 激情 狠狠 欧美| 欧美xxxx性猛交bbbb| 美女被艹到高潮喷水动态| 亚洲av日韩在线播放| 久久久久久久久久黄片| 最近手机中文字幕大全| 亚洲三级黄色毛片| 熟女电影av网| 国产av不卡久久| 国产亚洲精品久久久com| 国产综合懂色| 欧美区成人在线视频| 男人舔女人下体高潮全视频| 亚洲国产最新在线播放| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线 | 国产午夜福利久久久久久| 黄片wwwwww| 日本爱情动作片www.在线观看| 久久久国产成人精品二区| 久久亚洲国产成人精品v| 2021天堂中文幕一二区在线观| 久久久精品欧美日韩精品| 免费黄色在线免费观看| 少妇熟女aⅴ在线视频| 国产真实伦视频高清在线观看| 最近最新中文字幕免费大全7| 一区二区三区免费毛片| 亚洲人成网站在线播| 少妇高潮的动态图| 免费不卡的大黄色大毛片视频在线观看 | 汤姆久久久久久久影院中文字幕 | 亚洲av.av天堂| 综合色丁香网| 国产精品久久久久久精品电影小说 | 久久精品久久久久久久性| 能在线免费看毛片的网站| 精品人妻熟女av久视频| 高清在线视频一区二区三区 | 亚洲乱码一区二区免费版| 欧美日韩在线观看h| 99久久九九国产精品国产免费| 久久久久久久久久成人| 国产av码专区亚洲av| 小说图片视频综合网站| 亚洲自拍偷在线| 日本猛色少妇xxxxx猛交久久| 69av精品久久久久久| 美女cb高潮喷水在线观看| 内地一区二区视频在线| 欧美日韩国产亚洲二区| 22中文网久久字幕| 99热网站在线观看| 国产成年人精品一区二区| 日韩三级伦理在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品嫩草影院av在线观看| 麻豆av噜噜一区二区三区| 欧美成人午夜免费资源| 观看免费一级毛片| 免费看日本二区| 精品久久久久久久久久久久久| 久久久欧美国产精品| 中文亚洲av片在线观看爽| 少妇被粗大猛烈的视频| 国产黄色小视频在线观看| 国产精品久久久久久精品电影小说 | 免费搜索国产男女视频| 在线免费十八禁| 免费看av在线观看网站| 成人欧美大片| 亚洲中文字幕一区二区三区有码在线看| 日本色播在线视频| 国产欧美另类精品又又久久亚洲欧美| 少妇的逼好多水| 日本黄大片高清| 一区二区三区高清视频在线| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲最大av| 丝袜美腿在线中文| 国产精品美女特级片免费视频播放器| 成年免费大片在线观看| 国产日韩欧美在线精品| h日本视频在线播放| 汤姆久久久久久久影院中文字幕 | 全区人妻精品视频| 嫩草影院精品99| 免费看日本二区| 夫妻性生交免费视频一级片| 国产精品一区二区三区四区久久| 麻豆乱淫一区二区| 18+在线观看网站| 成人三级黄色视频| 亚洲欧洲国产日韩| 国产精品一区www在线观看| 中文字幕av在线有码专区| 国产亚洲5aaaaa淫片| 欧美zozozo另类| av在线老鸭窝| 日本色播在线视频| 国产精品久久久久久av不卡| 午夜视频国产福利| 国产伦精品一区二区三区视频9| 春色校园在线视频观看| 青春草视频在线免费观看| 国产爱豆传媒在线观看| 视频中文字幕在线观看| 亚洲在线自拍视频| 欧美bdsm另类| 国产精品久久久久久精品电影小说 | 亚洲成人中文字幕在线播放| 97热精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 欧美激情国产日韩精品一区| 热99re8久久精品国产| 精品人妻一区二区三区麻豆| 丝袜美腿在线中文| 在线观看美女被高潮喷水网站| 欧美丝袜亚洲另类| 午夜免费男女啪啪视频观看| 日韩制服骚丝袜av| 99热6这里只有精品| 国产伦理片在线播放av一区| 国产精品,欧美在线| 水蜜桃什么品种好| 亚洲经典国产精华液单| 嫩草影院精品99| 亚洲美女搞黄在线观看| 男人舔奶头视频| 91精品伊人久久大香线蕉| 亚洲图色成人| 国语自产精品视频在线第100页| 97人妻精品一区二区三区麻豆| 国产成人a∨麻豆精品| 欧美97在线视频| 日韩人妻高清精品专区| 亚洲精品色激情综合| 国产精品一区www在线观看| 看黄色毛片网站| 九九在线视频观看精品| 国产精品精品国产色婷婷| 色尼玛亚洲综合影院| 麻豆久久精品国产亚洲av| 18禁在线播放成人免费| 男女那种视频在线观看| 91av网一区二区| 国产精品久久视频播放| 99久久精品一区二区三区| 日本色播在线视频| 狂野欧美白嫩少妇大欣赏| 黄色配什么色好看| 五月玫瑰六月丁香| 国产精品嫩草影院av在线观看| 超碰97精品在线观看| 国产欧美日韩精品一区二区| 中文在线观看免费www的网站| 最近2019中文字幕mv第一页| 麻豆成人av视频| 亚洲在久久综合| 婷婷六月久久综合丁香| 在线观看66精品国产| 丰满人妻一区二区三区视频av| 日韩强制内射视频| 久久久国产成人免费| 国产综合懂色| 一个人看的www免费观看视频| 婷婷六月久久综合丁香| 欧美激情在线99| 男女边吃奶边做爰视频| 亚洲av男天堂| 日日撸夜夜添| 丝袜美腿在线中文| 欧美成人免费av一区二区三区| 伦精品一区二区三区| av.在线天堂| 久久久久久久久久久免费av| av国产免费在线观看| 91精品伊人久久大香线蕉| 美女内射精品一级片tv| 日韩亚洲欧美综合| videos熟女内射| 亚洲av中文字字幕乱码综合| 久久99精品国语久久久| 精品免费久久久久久久清纯| 特大巨黑吊av在线直播| 成人欧美大片| 波多野结衣高清无吗| 国产一区亚洲一区在线观看| av视频在线观看入口| 亚洲欧美精品综合久久99| 九草在线视频观看| 亚洲欧洲国产日韩| 精品久久久久久久末码| 国产在视频线在精品| 亚洲中文字幕日韩| 男女边吃奶边做爰视频| 在线免费观看不下载黄p国产| 国产精品一区二区在线观看99 | 久久久久久久亚洲中文字幕| 少妇丰满av| 久久精品夜色国产| 身体一侧抽搐| 日韩高清综合在线| 欧美不卡视频在线免费观看| 久久久久精品久久久久真实原创| 午夜激情欧美在线| 国产高潮美女av| av视频在线观看入口| 波多野结衣巨乳人妻| 精品久久久久久久久亚洲| 欧美xxxx性猛交bbbb| 一级黄色大片毛片| 毛片一级片免费看久久久久| 国产又色又爽无遮挡免| 国产av不卡久久| 久久久久久久久久久丰满| 日本五十路高清| 亚洲av免费高清在线观看| 少妇熟女aⅴ在线视频| 国产国拍精品亚洲av在线观看| 老司机福利观看| 久久久久久九九精品二区国产| 男女视频在线观看网站免费| 熟妇人妻久久中文字幕3abv| 国产精品日韩av在线免费观看| 亚洲美女搞黄在线观看| 亚洲成人中文字幕在线播放| 麻豆精品久久久久久蜜桃| 在线天堂最新版资源| 岛国在线免费视频观看| 久久精品国产亚洲av涩爱| 国国产精品蜜臀av免费| 啦啦啦观看免费观看视频高清| 九九久久精品国产亚洲av麻豆| 国产麻豆成人av免费视频| 国产老妇女一区| 我要搜黄色片| 欧美高清成人免费视频www| 国产高清视频在线观看网站| 国产亚洲av片在线观看秒播厂 | 亚洲欧美成人精品一区二区| 在线观看一区二区三区| 18禁在线无遮挡免费观看视频| 日韩欧美国产在线观看| 日韩精品青青久久久久久| 青春草亚洲视频在线观看| a级毛色黄片| 黄色配什么色好看| 国产男人的电影天堂91| 色视频www国产| 亚洲色图av天堂| 亚洲成人中文字幕在线播放| 久久久色成人| 少妇熟女aⅴ在线视频| 美女国产视频在线观看| 人妻系列 视频| 国产精品国产高清国产av| 我的女老师完整版在线观看| 少妇熟女欧美另类| 一区二区三区高清视频在线| 国内精品宾馆在线| 欧美最新免费一区二区三区| 久久午夜福利片| 日韩,欧美,国产一区二区三区 | 女人久久www免费人成看片 | 大又大粗又爽又黄少妇毛片口| 成人亚洲精品av一区二区| 日本熟妇午夜| 在线观看一区二区三区| 丰满少妇做爰视频| 18+在线观看网站| 三级男女做爰猛烈吃奶摸视频| 波多野结衣高清无吗| 久久久久九九精品影院| 尾随美女入室| 精华霜和精华液先用哪个| 成人午夜高清在线视频| 久久久久久久久久久丰满| 伦精品一区二区三区| 久久99热这里只有精品18| 亚洲国产精品成人久久小说| 国产探花极品一区二区| 国产精品一及| 国产一级毛片在线| 永久网站在线| 欧美性猛交╳xxx乱大交人| 国产91av在线免费观看| 欧美不卡视频在线免费观看| 少妇的逼好多水| 91精品国产九色| 精品一区二区三区人妻视频| 国产黄色小视频在线观看| 99热精品在线国产| 国产一区二区在线观看日韩| 亚洲欧美日韩高清专用| 熟妇人妻久久中文字幕3abv| 亚洲国产欧美在线一区| 99热这里只有精品一区| 国产亚洲5aaaaa淫片| 丰满少妇做爰视频| a级一级毛片免费在线观看| 久久国产乱子免费精品| 亚洲欧美一区二区三区国产| 国产精品.久久久| 日韩强制内射视频| 久久久久久久亚洲中文字幕| 久久精品国产自在天天线| 麻豆乱淫一区二区| 欧美日韩一区二区视频在线观看视频在线 | 最近最新中文字幕免费大全7| 久久久国产成人精品二区| 国产成人精品久久久久久| 男女视频在线观看网站免费| 国产午夜精品论理片| 国产欧美另类精品又又久久亚洲欧美| 一卡2卡三卡四卡精品乱码亚洲| 小蜜桃在线观看免费完整版高清| 久久精品夜夜夜夜夜久久蜜豆| 99热全是精品| 三级男女做爰猛烈吃奶摸视频| 最近视频中文字幕2019在线8| 真实男女啪啪啪动态图| 精品国产三级普通话版| 国产亚洲一区二区精品| 日本免费一区二区三区高清不卡| 性插视频无遮挡在线免费观看| 亚洲在线观看片| 国产爱豆传媒在线观看| 一区二区三区乱码不卡18| 欧美日韩国产亚洲二区| 成年女人永久免费观看视频| av在线老鸭窝| 搡女人真爽免费视频火全软件| 欧美3d第一页| 内射极品少妇av片p| 免费电影在线观看免费观看| 亚洲欧美精品综合久久99| 久久久午夜欧美精品| 日韩高清综合在线| 女人久久www免费人成看片 | 久久这里有精品视频免费| 欧美bdsm另类| 乱系列少妇在线播放| 欧美高清成人免费视频www| 亚洲av福利一区| 蜜桃久久精品国产亚洲av| 欧美日韩在线观看h| 日韩一区二区三区影片| 日本免费一区二区三区高清不卡| 2021少妇久久久久久久久久久| 日本熟妇午夜| 国产精品日韩av在线免费观看| 国产成人精品一,二区| 日韩精品有码人妻一区| 99热这里只有是精品在线观看| 亚洲精品aⅴ在线观看| 一级毛片aaaaaa免费看小| 三级国产精品欧美在线观看| 观看美女的网站| 中文字幕av在线有码专区| 国产伦精品一区二区三区视频9| 一级黄片播放器| 久久久亚洲精品成人影院| 久久久精品欧美日韩精品| 大又大粗又爽又黄少妇毛片口| 亚洲精品乱久久久久久| 久热久热在线精品观看| 只有这里有精品99| 在线免费观看不下载黄p国产| 1000部很黄的大片| 大香蕉97超碰在线| 亚洲在久久综合| 在线观看美女被高潮喷水网站| 国产精品野战在线观看| 成人国产麻豆网| 国产精品蜜桃在线观看| 国产三级在线视频| 免费黄色在线免费观看| 中文天堂在线官网| 少妇猛男粗大的猛烈进出视频 | 3wmmmm亚洲av在线观看| 老司机影院成人| 日本黄色视频三级网站网址| 亚洲精品自拍成人| 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| 免费人成在线观看视频色| 在线免费观看的www视频| 中文在线观看免费www的网站| 亚洲怡红院男人天堂| 国产精品伦人一区二区| 久久精品人妻少妇| 国产亚洲av嫩草精品影院| 在线观看66精品国产| 亚洲高清免费不卡视频| 国产精品.久久久| 九九热线精品视视频播放| 波多野结衣巨乳人妻| a级毛片免费高清观看在线播放| 国产精品野战在线观看| 国国产精品蜜臀av免费| 国产激情偷乱视频一区二区| 国产亚洲91精品色在线| 国产色爽女视频免费观看| 国产成人免费观看mmmm| 国产精品蜜桃在线观看| 婷婷色综合大香蕉| 嘟嘟电影网在线观看| 国产精品无大码| 中文字幕精品亚洲无线码一区| 国产高清不卡午夜福利| 亚洲av男天堂| 久久人妻av系列| 亚洲欧美成人综合另类久久久 | 日韩中字成人| 国产中年淑女户外野战色| 天天一区二区日本电影三级| 午夜亚洲福利在线播放| 伊人久久精品亚洲午夜| 长腿黑丝高跟| 亚洲在线观看片| 又粗又硬又长又爽又黄的视频| 国产精品麻豆人妻色哟哟久久 | 中国国产av一级| 久久精品熟女亚洲av麻豆精品 | 99久久九九国产精品国产免费| 国产av不卡久久| 久久精品国产鲁丝片午夜精品| 国产成人免费观看mmmm| 国产又黄又爽又无遮挡在线| 国产精品久久久久久精品电影| 尤物成人国产欧美一区二区三区| 国产不卡一卡二| 亚洲国产精品久久男人天堂| 热99在线观看视频| 久久久国产成人精品二区| 少妇人妻一区二区三区视频| 成年免费大片在线观看| 人人妻人人澡欧美一区二区| 国产视频内射| 久久久亚洲精品成人影院| 国内精品宾馆在线| 精品人妻偷拍中文字幕| av女优亚洲男人天堂| 国产爱豆传媒在线观看| 亚洲成人久久爱视频| 亚洲精品久久久久久婷婷小说 | 少妇的逼水好多| 有码 亚洲区| 国产精品伦人一区二区| 亚洲精品乱码久久久久久按摩| 国内精品宾馆在线| 亚洲中文字幕一区二区三区有码在线看| 男插女下体视频免费在线播放| 插逼视频在线观看| 日韩av在线大香蕉| 黄色一级大片看看| 日韩av在线大香蕉| 欧美不卡视频在线免费观看| 国产黄片视频在线免费观看| av在线亚洲专区| 成人漫画全彩无遮挡| 看片在线看免费视频| 亚洲精品乱码久久久久久按摩| 久热久热在线精品观看| 中文精品一卡2卡3卡4更新| 特级一级黄色大片| 国模一区二区三区四区视频| 日韩av不卡免费在线播放| 蜜桃亚洲精品一区二区三区| 爱豆传媒免费全集在线观看| 国产伦理片在线播放av一区| 一本久久精品| 91aial.com中文字幕在线观看| 免费黄网站久久成人精品| 色哟哟·www| 黄色配什么色好看| 简卡轻食公司| eeuss影院久久| videos熟女内射| 欧美性感艳星| 午夜福利在线观看免费完整高清在| 水蜜桃什么品种好| 九九久久精品国产亚洲av麻豆| 少妇的逼水好多| 久久午夜福利片| 亚洲精品日韩在线中文字幕| 国产久久久一区二区三区| 97在线视频观看| 黄色日韩在线| 国产中年淑女户外野战色| 亚洲国产精品专区欧美| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看吧| 男人舔奶头视频| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久午夜电影| 中文资源天堂在线| 国模一区二区三区四区视频| 国产三级在线视频| 精品久久久噜噜| 免费观看性生交大片5| 麻豆久久精品国产亚洲av| 国产精品久久久久久久久免| 亚洲av成人av| 亚洲在久久综合| 亚洲国产精品成人久久小说| 亚洲av成人精品一二三区| 久久久欧美国产精品| 毛片一级片免费看久久久久| 日韩欧美 国产精品| 蜜桃久久精品国产亚洲av| 男女国产视频网站| 女人久久www免费人成看片 | 久久久久精品久久久久真实原创| 99久久精品一区二区三区| 美女内射精品一级片tv| 日韩一本色道免费dvd| 国产精品永久免费网站| 欧美性感艳星| av在线蜜桃| 亚洲av不卡在线观看| 亚洲av免费高清在线观看| 五月伊人婷婷丁香| 波多野结衣巨乳人妻| 插阴视频在线观看视频| 欧美日韩在线观看h| 99热精品在线国产| 亚洲五月天丁香| av播播在线观看一区| 麻豆av噜噜一区二区三区| 国产成年人精品一区二区| 免费看光身美女| 国产老妇伦熟女老妇高清| 国产不卡一卡二| 免费av观看视频| 免费观看的影片在线观看| 免费看a级黄色片| 一级爰片在线观看| 亚洲欧美日韩卡通动漫| 久久欧美精品欧美久久欧美| 国产一区有黄有色的免费视频 | 美女cb高潮喷水在线观看| 国产极品天堂在线| 99热这里只有是精品50| 久久精品久久精品一区二区三区| 22中文网久久字幕| 亚洲精品一区蜜桃| 超碰av人人做人人爽久久| 中文天堂在线官网| videos熟女内射| 亚洲av熟女| 久久欧美精品欧美久久欧美| 国产成人a区在线观看| 久久6这里有精品| 午夜日本视频在线|