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

    華南亞熱帶武夷山地區(qū)垂直植被帶表土孢粉散布規(guī)律

    2023-07-13 11:41:50龍香月黃康有謝德豪水坤春李宏衛(wèi)
    熱帶地理 2023年6期
    關(guān)鍵詞:孢粉禾本科武夷山

    龍香月,黃康有,2,陳 聰,謝德豪,水坤春,李宏衛(wèi),鄭 卓

    (1.中山大學(xué) 地球科學(xué)與工程學(xué)院//廣東省地球動(dòng)力作用與地質(zhì)災(zāi)害重點(diǎn)實(shí)驗(yàn)室,廣東 珠海 519082;2.南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室,廣東 珠海 519082;3.廣東南嶺國家級(jí)自然保護(hù)區(qū),廣東 韶關(guān) 512727;4.廣東省地質(zhì)調(diào)查院,廣州 510080)

    第四紀(jì)古植被與古氣候定量重建研究中,孢粉分析具有不可替代的作用,但由于不同植物孢粉在產(chǎn)量、傳播途徑、散布能力、搬運(yùn)方式、沉降速率和保存條件等方面都存在差異,從而導(dǎo)致孢粉組合中的孢粉種類數(shù)量、種類的含量(某種孢粉類型數(shù)目與孢粉總數(shù)的比值)與實(shí)際植被中該植物的數(shù)量和比例不一定完全一致,即孢粉與植被-氣候之間并不是簡單的線性關(guān)系,導(dǎo)致利用化石孢粉定量重建古氣候、古環(huán)境存在較大不確定性(Mauri et al.,2015)。因此正確解釋現(xiàn)代孢粉與植被-氣候之間的關(guān)系,不僅是準(zhǔn)確分析和解釋化石孢粉數(shù)據(jù)的重要依據(jù),而且為利用化石孢粉資料定量重建古氣候和古環(huán)境提供重要參考(鄭卓 等,2013)。

    近年來,現(xiàn)代孢粉過程研究越來越受重視,國內(nèi)學(xué)者對(duì)表土孢粉進(jìn)行了大量工作,研究區(qū)覆蓋中國東北(李曼玥 等,2012;余少華 等,2012;Han et al., 2020)、華北(Li et al., 2017)、青藏高原(Wei et al., 2011;李怡雯 等,2019)、西北(Qin et al., 2015; Wei et al., 2015; Huang et al., 2018;吳慧賢 等,2018)、西南(汪啟容 等,2021)、華中(Fang et al., 2015)以及華東(舒軍武 等,2010;李圓圓 等,2016;郝秀東 等,2020;全曉文 等,2022)等地區(qū)。主要探討現(xiàn)代孢粉與現(xiàn)代植被-氣候的關(guān)系、基于花粉數(shù)據(jù)下的植物群系分類研究、花粉代表性問題、利用孢粉學(xué)研究植被受到的人為影響和利用花粉數(shù)據(jù)進(jìn)行古氣候、古環(huán)境重建等內(nèi)容。表土孢粉與現(xiàn)代植被和氣候之間的關(guān)系研究揭示,現(xiàn)代表土孢粉具有較好的植被和氣候指示性,孢粉組合與現(xiàn)代植被之間具有一致性,基本可以反映現(xiàn)代植被的主要特征以及對(duì)應(yīng)的氣候類型(Wei et al., 2011;楊慶華 等,2019)。

    武夷山脈地區(qū)是中國生物多樣性的關(guān)鍵區(qū)域之一(丁暉 等,2015;張玉琴 等,2020)。前人在該地區(qū)已開展一系列沉積孢粉和古環(huán)境變化研究,戴云山泥炭孢粉分析結(jié)果揭示了該地區(qū)自中晚全新世以來的氣候變化和植被演化過程(邱宏烈,2006);此外,武夷山脈地區(qū)泥炭沼澤沉積巖芯中的花粉和碳屑記錄也揭示了該地區(qū)距今8200年以來的植被演化、氣候變化以及人類活動(dòng)對(duì)植被景觀的影響(Ma et al., 2016),而東部地區(qū)的水竹洋鉆孔孢粉記錄了該地區(qū)從末次冰消期至全新世以來的植被演變過程(Yue et al., 2012; Wan et al., 2021)。然而,武夷山地區(qū)表土孢粉散布規(guī)律與植被-氣候之間相關(guān)關(guān)系的研究成果相對(duì)較少,利用化石孢粉資料定量重建該地區(qū)古環(huán)境和古氣候仍需研究現(xiàn)代孢粉與植被之間相互關(guān)系的準(zhǔn)確性。因此,本文擬通過對(duì)武夷山脈76個(gè)不同海拔的表土孢粉樣品進(jìn)行分析,探討該地區(qū)表土孢粉的垂直分布規(guī)律及其與植被之間的關(guān)系,以期為準(zhǔn)確利用化石孢粉重建武夷山及周邊地區(qū)的古環(huán)境、古氣候提供基礎(chǔ)資料。

    1 材料與方法

    1.1 研究區(qū)概況與植被特征

    武夷山脈位于閩贛兩省之間(24°30′-28°20′N,115°33′-118°50′ E),山脈走向呈東北—西南方向,北起浙江省東南部,橫亙福建省西北部,南至廣東省南嶺地區(qū),是汀江水系、鄱陽湖水系與閩江水系的分水嶺。整體山脈地貌起伏大、地勢(shì)高差懸殊,主峰黃崗山位于北段,最高峰海拔為2 158 m。氣候上屬于中亞熱帶季風(fēng)氣候區(qū),受來自太平洋的東亞夏季風(fēng)影響,該地區(qū)年均降水量在2 000 mm左右,年均溫8.5~18 ℃,相對(duì)濕度在78%~84%之間,無霜期253~272 d,1月平均氣溫7.8℃,7月平均氣溫27.8℃(葉文 等,2007)。土壤由山麓至山頂主要有山地黃紅壤(400~600 m)、山地黃壤(600~1 300 m)、山地黃褐壤(1 300~1 900 m)和山地草甸土(>1 900 m)(何容 等,2009)。

    武夷山脈北段的國家級(jí)自然保護(hù)區(qū)植被發(fā)育狀況較為良好且保存完整,保存著中國東南地區(qū)面積最大、最完整的中亞熱帶森林生態(tài)系統(tǒng),具有世界同緯度地區(qū)現(xiàn)存大陸面積最大且保存最完整的中亞熱帶森林生態(tài)系統(tǒng),植被覆蓋率達(dá)96.3%,植被變化規(guī)律明顯。該地區(qū)典型的地帶性植被類型為亞熱帶常綠闊葉林(陳昌篤,1999)。

    武夷山脈植被的垂直地帶性分異顯著,從山麓向山頂依次分布著常綠闊葉林(200 ~1 200 m)、針闊混交林(1 000~1 350 m)、針葉林(1 350~1 750 m)、中山闊葉矮曲林(1 750~1 900 m)和山頂草甸(1 900~2 158 m)等類型(葉文 等,2007;丁暉等,2015)。此外,由于受到人類活動(dòng)的影響,武夷山山麓廣泛分布著農(nóng)作物及荒坡灌叢(權(quán)偉 等,2008)(圖1)。

    1)人工林(450 m以下)

    該植被帶分布在海拔450 m 以下,原始植被幾乎消失殆盡。以人工植被為主,如馬尾松(Pinus massoniana)、杉木(Cunninghamia lanceolata)和毛竹(Phyllostachys edulis)等,自然灌叢有烏飯(Vaccinium bracteatum)和杜鵑(Rhododendronsp.)等,草本植物主要有芒萁(Dicranopteris dichotoma);農(nóng)作物主要有水稻、玉米、小麥和茶樹等。

    2)常綠闊葉林(200~1 200 m)

    常綠闊葉林是武夷山垂直植被帶的基帶,主要由殼斗科的椎屬(Castanopsis)、石櫟屬(Lithocarpus)、青岡屬(Cyclobablanopsis)的甜櫧(C.eyrei)、米櫧(Castanopsis carlesii)、細(xì)葉青岡(Cyclobalanopsis tenuifolia)、樟科(Lauraceae)、山茶科(Theaceae)等樹種組成。群落結(jié)構(gòu)可分為喬木層、灌木層和草本層3個(gè)基本層次,其中喬木層主要有米櫧(Castanopsis carlesii)、甜櫧(Castanopsis eyrei)、毛錐(C.fordii)、南嶺栲(C.fordii)、羅浮栲(C.fabri)、少葉黃杞(Engelhardtia fenzelii)、赤楠(Syzygium buxifolium)、木荷(Schima superba)、細(xì)葉青岡(Cyclobalanopsistenuifolia)、格藥柃(Eurya muricata)等常綠闊葉樹種;灌木層主要有赤楠(Syzygium buxifolium)、微毛柃(Eurya hebeclados)、尖連蕊茶(Camellia cuspidata)、彎蒴杜鵑 (Rhododendron henryi)等;草本層主要有中華里白(Hicriopteri chinensis)、狗脊蕨(Woodwardia japonica)、芒萁(Dicranopteris dichotoma)、苔草(Carex tristachya)、蕨狀薹草(Carex filicina)和草珊瑚(Sarcandra glabra)。此外,該植被帶約海拔800 m 下分布有較大面積的人工毛竹林。

    3)針闊混交林(1 000~1 350 m)

    針闊葉林以黃山松(Pinus tanwanensis)、南方鐵杉(Tsuga chinensis)等針葉樹與殼斗科的甜櫧(Castanopsis eyrei)、木荷(Schima superba)、山茶科的格藥柃(Eurya muricata)、樟科等闊葉樹種混交為主。喬木層以黃山松(Pinus tanwanensis)、木荷(Schima superba)、甜櫧(Castanopsis eyrei)、格藥柃(Eurya muricata)、青岡(Cyclobalanopsis)、多脈青岡(Cyclobalanopsis multinervis)等植被為主;灌木層主要分布有格藥柃(Eurya muricata)、馬銀花(Rhododendron ovatum)、赤楠(Syzygium buxifolium);草本層主要分布有里白(Hicriopteris)、芒萁(Dicranopteris dichotoma)、苔草 (Carex tristachya)。

    4)針葉林(1 350~1 750 m)

    黃山松(Pinus taiwanensis)、南方鐵杉(Tsuga chinensis)和杉木(Cunninghamia lanceolata)等針葉樹構(gòu)成針葉林帶,此外主要常見薔薇科(Rosaceae)、禾本科(Poaceae)、山茶科、杜鵑花科(Ericaceae)及殼斗科等植被。喬木層主要為黃山松(Pinus taiwanensis)、杉木(Cunninghamia lanceolata)、南方鐵杉(Tsuga chinensis)、柳杉(Cryptomeria)等;灌木層分布有細(xì)齒葉柃(Eurya nitida)、小葉石楠(Photinia parvifolia)、云錦杜鵑(Rhododendron fortunei)、馬銀花(Rhododendron ovatum)、老鼠矢(Symplocos stellaris);草本層主要有芒(Miscanthus)、苔草(Carex tristachya)、青茅(Calamagrostis brachytricha)等。

    5)中山闊葉矮曲林(1 750~1 900 m)

    中山闊葉矮曲林以落葉樹種為主,苔蘚植物繁多。喬木層主要有鐵杉(Tsuga heterophylla)、木荷(Schima superba)、多脈青岡(Cyclobananopsis multinervis),灌木層有白檀(Symplocos paniculata)、云錦杜鵑(Rhododendron fortunei)等;草本層分布有龍膽(Radix gentianae)、青茅(Calamagrostis brachytricha)和野菊花(Chrysanthemum indicum)。

    6)山頂草甸(1 900~2 158 m)

    中山山頂草甸群落結(jié)構(gòu)簡單,以禾本科青茅為主。灌木層有白檀(Symplocos paniculata)、云錦杜鵑(Rhododendron fortunei)、波葉紅果樹(Stranvaesia davidiana);草本層以青茅(Calamagrostis brachytricha)為主,有少量野菊花(Chrysanthemum indicum)、龍膽(Radix gentianae)、藜蘆(Veratrum nigrum)、石松(Murraya paniculata)等;層外植物為金銀花(Lonicera japonica)。此外,草甸上還零星散布著黃山松。

    1.2 樣品采集與實(shí)驗(yàn)室處理

    本研究的表土孢粉樣品是團(tuán)隊(duì)在過去20年間多次采集積累所得,表土孢粉樣品大部分選取天然植被群落下方或者人類活動(dòng)影響較小的植被群落的表層土壤,表土樣品采集按梅花點(diǎn)狀采集方法,即在植被樣方(面積20 m×20 m)4 個(gè)角及中心位置表層1~2 cm取土,充分混合后作為1個(gè)樣品,取樣過程中記錄每個(gè)樣點(diǎn)的經(jīng)緯度和海拔高度,總共選擇76個(gè)現(xiàn)代表土孢粉樣品數(shù)據(jù)展開散布規(guī)律及氣候因子特征分析(表1)。

    表1 武夷山地區(qū)表土孢粉樣品采樣點(diǎn)信息及植被類型Table 1 Location and vegetation information of 76 surface pollen sampled in the Wuyi mountains

    表土孢粉樣品實(shí)驗(yàn)室分析采用常規(guī)的酸堿處理和重液懸浮法,每個(gè)樣品稱取約5 g放置于燒杯中,并加入1片石松孢子(27 560粒/片)以估算孢粉濃度,樣品先用10%HCl除去鈣質(zhì)膠結(jié);洗至中性后加入10%的KOH除去有機(jī)質(zhì);再一次洗至中性后,用密度為2.0 g/cm3的重液進(jìn)行浮選,將收集的孢粉加入加甘油保存,制成活動(dòng)片鑒定。孢粉形態(tài)鑒定和統(tǒng)計(jì)在蔡司光學(xué)顯微鏡Axio Scope A1 下完成,主要參考《中國熱帶亞熱帶被子植物花粉形態(tài)》(中國科學(xué)院植物研究所古植物室孢粉組,1982)和《中國植物花粉形態(tài)》(王伏雄 等,1995)等圖集,每個(gè)樣品至少鑒定統(tǒng)計(jì)孢粉400粒以上。孢粉種類百分比以陸生植物花粉總和為基數(shù)進(jìn)行計(jì)算,運(yùn)用Tilia 軟件進(jìn)行CONISS 聚類分析并繪制孢粉圖譜。

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

    冗余分析(Redundancy Analysis, RDA)和典型對(duì)應(yīng)分析(Canonical Correspondence Analysis,CCA)是物種與環(huán)境關(guān)聯(lián)分析中最常用的2種模型。為探討表土孢粉分布規(guī)律與現(xiàn)代植被以及氣候因子之間的關(guān)系,將孢粉組合中的典型代表性花粉(含量>0.2%)如柯-栲屬、松屬、櫟屬、楓香樹屬、大戟科、山茶科、薔薇科、禾本科、蒿屬等35個(gè)種類的含量作為表土孢粉數(shù)據(jù)矩陣在CANOCA 軟件中進(jìn)行降趨對(duì)應(yīng)分析(Descending Trend Correspondence Analysis, DCA),根據(jù)分析結(jié)果中梯度長度第一軸的值進(jìn)行下一步分析,當(dāng)?shù)谝惠S的值>4.0時(shí),選擇CCA;當(dāng)?shù)谝惠S的值<3.0 時(shí),選擇RDA;當(dāng)SD 值介于3.0~4.0 時(shí),選擇RDA 或CCA 進(jìn)行下一步分析均可(Lepx et al., 2003; Morris, 2015)。孢粉種類散布規(guī)律與氣候因子的相關(guān)性分析,運(yùn)用Canoco 4.5 軟件進(jìn)行典型對(duì)應(yīng)分析和相關(guān)性分析(Lepx et al., 2003)。

    本文使用的環(huán)境指標(biāo)數(shù)據(jù)有年均降雨量(Mean Annual Precipitation, MAP)、年均溫(Mean Annual Temperature, MAT)、年均相對(duì)濕度(Mean Annual Relative Humidity, MAH)、最冷月降水(Mean precipitation of coldest month, Mpco)、最冷月溫度(Mean temperature of coldest month, Mtco)和最冷月濕度(Mean humidity of coldest month,Mhco)以及最熱月降水(Mean precipitation of warmest month, Mpwa)、最熱月溫度(Mean temperature of warmest month, Mtwa)和最熱月濕度(Mean humidity of warmest month, Mhwa),氣候指標(biāo)基于中國756個(gè)氣象臺(tái)站的氣象數(shù)據(jù)運(yùn)用薄板樣條回歸插值方法進(jìn)行插值獲得。

    2 結(jié)果分析

    武夷山地區(qū)76個(gè)表土孢粉樣品總共鑒定出141個(gè)科(屬),其中喬、灌木128 科(屬),草本8 科(屬)以及蕨類孢子5科(屬)。主要孢粉類型:針葉樹花粉有松屬(Pinus)、鐵杉屬(Tsuga)、杉科(Taxodiaceae)等;闊葉樹花粉有栲屬-柯屬(Castanopsis-Lithocarpus),常綠櫟屬(Quercus-evergreen)、落葉櫟屬(Quercus-deciduous)楓香樹屬(Liquidambar)、大戟科(Euphorbiaceae)和山茶科(Theaceae)等;灌木花粉有薔薇科(Rosaceae)、杜鵑花科(Ericaceae)、楊梅屬(Myrica)和榛屬(Corylus)等;陸生草本植物花粉中禾本科(Poaceae)、蒿屬(Artemisia)、菊科(Asteraceae)、莎草科(Cyperaceae)和茜草科(Rubiaceae)等較為常見;蕨類孢子種類主要有芒萁屬(Dicranopteris)、三縫孢子(Trilete)、單縫孢子(Monolete)、石松屬(Lycopodium)和鳳尾蕨屬(Pteris)等。

    在表土孢粉組合中,喬灌木花粉種類豐富且含量高,占全部花粉總數(shù)的64.5%,其中針葉樹、闊葉樹和灌木的花粉含量分別為21.1%、40.3%和3.1%,針葉樹花粉含量波動(dòng)幅度大,介于1.1%~95.8%,以松屬為主,平均含量為19.9%;闊葉樹種類豐富,其中栲屬-柯屬含量最高,平均含量為23.2%,最高達(dá)91.4%,此外常綠櫟屬、楓香樹屬、蕓香科(Rutaceae)和木荷等植物花粉的平均含量分別為8.9%、1.7%,0.6%和0.6%;灌木植物花粉在所有植被帶中都較少出現(xiàn)(平均含量3.1%),以薔薇科和杜鵑花科為主,平均含量分別為1.1%和0.7%;陸生草本植物種類相對(duì)較少,且花粉含量較低(平均含量為8.6%),以禾本科(0.01%~58.3%)和蒿屬(0.01%~20.6%)的含量相對(duì)較高,平均為5.1%和1.1%;蕨類孢子在所有樣品中均有出現(xiàn),其中芒萁屬(0~68.1%)的含量最為豐富,平均為13.7%,三縫孢子、單縫孢子、石松屬在樣品中不同程度的出現(xiàn),平均含量分別為8.7%、2.8%和1.6%。

    2.1 垂直植被帶孢粉組合特征

    根據(jù)表土孢粉樣點(diǎn)海拔從低到高作圖,同時(shí)在結(jié)合聚類結(jié)果的基礎(chǔ)上,將武夷山地區(qū)表土孢粉組合劃分為5個(gè)孢粉帶(圖2),各孢粉組合帶的特征描述為:

    圖2 武夷山地區(qū)表土孢粉含量圖譜Fig.2 Pollen diagram of surface pollen samples in the Wuyi mountains

    2.1.1 孢粉帶I(低地常綠闊葉林帶)該孢粉帶包括15個(gè)樣品,海拔在450 m以下。孢粉組合以闊葉樹植物花粉(6.3%~97.4%,平均為49.6%)占優(yōu)勢(shì),其次為蕨類孢子(6.4%~75.5%,平均為27.3%)、針葉樹花粉(1.1%~85.9%,平均為14.1%),草本植物花粉(0.3%~44.6%,平均為7.6%)和灌木植物花粉(0.01%~7.4%,平均為1.4%)較少。闊葉樹植物花粉以殼斗科栲屬-柯屬含量最高(0.5%~90.7%),平均為28.9%,常綠櫟屬、山茶科、木荷屬和冬青屬(Ilex)等含量較高,平均分別為4.9%、3.7%、1.8%和1.3%,伴有大戟科(0.01%~8.4%,平均為0.8%)、蕓香科(0~4.3%,平均為0.8%)、漆樹科(Anacardiaceae,0.01%~5.8%,平均為0.5%)等常綠闊葉樹種,落葉闊葉樹以鵝耳櫪屬(Carpinus,0.01%~70.7%,平均為3.4%)和楓香樹屬(0.01%~9.1%,平均為0.8%)為主,混生少量的樺木(Betula)等;蕨類孢子在本帶含量較高,主要有芒萁屬(0.01%~68.5%),平均為14.7%,其次為三縫孢子、單縫孢子等,平均含量分別為7.1%和4.7%。

    2.1.2 孢粉帶II(山地常綠闊葉林帶)該孢粉帶包含海拔在450~1 000 m 之間的26個(gè)樣品,孢粉組合中闊葉樹植物花粉(2.4%~94.6%,平均為44.1%)含量略有下降,常綠櫟屬含量下降至3.7%,其余常綠闊葉樹種花粉含量也有不同程度的降低,但降幅較小,而栲屬-柯屬的含量(1.2%~91.4%,平均為31.4%)略有增加,仍為優(yōu)勢(shì)種,落葉闊葉樹花粉含量小幅度增加,如楓香樹屬和樺木屬在本帶的平均含量分別上升至3.5%和0.2%,同時(shí)落葉櫟屬花粉在本孢粉帶有一定的含量,平均約為0.4%。總體而言,孢粉帶II 闊葉樹花粉含量豐富,組成特征與帶I 相似;蕨類孢子含量(10.1%~68.6%,平均27.6%)較高,以芒萁為主,平均含量為15.1%;針葉樹花粉含量(1.5%~67.3%)顯著上升,平均為21.2%,以松屬(1.4%~91.4%,平均為20.4%)為主,含少量鐵杉,約占0.5%;灌木花粉含量(0.01%~19.0%)略增,但仍較低,平均為2.3%,主要有薔薇科(0.01%~19.0%)和杜鵑科(0.01%~5%)等,平均分別為1.0%和0.9;草本植物花粉(0.8%~19.4%,平均為4.9%)與帶I相比明顯減少,以禾本科(0.5%~9.7%,平均為2.2%)、莎草科(0.01%~4.5%,平均為0.5%)、蒿屬(0.01%~3.8%,平均為1.0%)的含量降低為主,其中禾本科降幅最為顯著;蕨類孢子中,三縫孢子(0.1%~66.8%,平均為9.9%)的含量明顯高于單縫孢子(1.0%~8.1%,平均為1.8%)。

    2.1.3 孢粉帶III(山地常綠-落葉混交林帶)該孢粉帶包括15 個(gè)樣點(diǎn),海拔為1 000~1 350 m,優(yōu)勢(shì)孢粉種類為松屬、栲屬-柯屬、櫟屬和杉科等,以柯-栲屬花粉含量顯著下降和常綠櫟屬花粉增加為特點(diǎn)區(qū)別于孢粉帶I和孢粉帶II。孢粉組合中闊葉樹花粉含量下降,平均為32.4%,比帶II 降低了11.7%,降幅明顯。主要表現(xiàn)為栲屬-柯屬由28.4%下降至16.4%,其含量變化明顯具有隨海拔升高而降低的趨勢(shì),但常綠櫟屬含量增加至7.9%,落葉櫟屬(1.0%)、樺木屬(0.9%)、金縷梅科(Hamamelidaceae,0.7%)和榿木屬(Alnus,0.5%)等落葉闊葉樹花粉含量相對(duì)于孢粉帶II也有不同程度的增加,但含量仍較低;針葉樹花粉(5.0%~77.8%,平均為22.5%)的占比保持穩(wěn)定,以松屬(21.5%)為主,含有一定量杉科(1.0%);灌木植物花粉含量小幅度增加,平均為3.6%,主要有薔薇科(1.1%)、野牡丹科(Melastomataceae,1.0%)、杜鵑花科(0.7%)等;草本植物繁盛,花粉含量(0.3%~69.3%,平均為13.0%)顯著上升,以禾本科(9.3%)為優(yōu)勢(shì)種,伴有蒿屬(1.7%)、菊科(0.6%)、莎草科(0.6%)、茜草科(0.4%)等;蕨類孢子含量(11.3%~71.3%,平均為28.5%)略有增加,以芒萁屬(13.5%)和三縫孢子(9.2%)為主,還有少量單縫孢子、石松屬,平均為3.1%和2.6%。

    2.1.4 孢粉帶IV(山地針-闊葉混交林帶)該孢粉帶海拔在1 350 m~1 850 m 之間,共采集到16 個(gè)樣品,常綠櫟屬、松屬、栲屬-柯屬、鐵杉屬、落葉櫟屬、薔薇科等為主要種屬。該帶孢粉組合以常綠櫟屬花粉含量顯著增加為明顯特征區(qū)別于其他植被帶。孢粉組合中喬灌木花粉占絕對(duì)優(yōu)勢(shì),占比明顯升高(76.2%),表現(xiàn)在闊葉樹植物花粉含量顯著上升,主要由于常綠櫟屬含量持續(xù)大幅增加,升高至32.8%,同時(shí)落葉闊葉樹中落葉櫟屬(1.5%)、楓香樹屬(0.9%)、水青岡屬(Fagus,0.7%)及鵝耳櫪屬(0.7%)等花粉含量都出現(xiàn)不同程度的上升,而栲屬-柯屬含量下降至9.0%;針葉樹中松屬的含量保持穩(wěn)定,平均為19.7%,鐵杉屬(0~7.9%,2.6%)含量在該孢粉帶達(dá)到最高值;以薔薇科(1.8%)和杜鵑花科(0.9%)為主的灌木類花粉含量較低,平均僅為3.4%。草本植物花粉(1.6%~17.9%,平均為5.0%)和蕨類孢子(12.2%~34.4%,平均為18.9%)的含量均降低,其中陸生草本植物花粉含量變化最明顯的是禾本科,平均為1.9%,蒿屬含量下降至0.9%;蕨類孢子仍以芒萁(8.3%)和三縫孢子(6.5%)占主要比例,單縫孢子(2.1%)和石松屬(1.7%)的含量略有降低。

    2.1.5 孢粉帶V(山頂灌叢植被帶)該孢粉帶分布海拔較高,約在1 850~2 200 m之間,包括4個(gè)孢粉樣品。孢粉組合以針葉樹植物花粉(20.0%~95.8%,平均31.6%)和蕨類孢子(30.5%~63.3%,平均為30.4%)占優(yōu)勢(shì),其次為草本植物花粉(0.01%~56.1%,平均為22.5%)和闊葉樹花粉(2.8%~20.9%,平均為10.9%)。闊葉樹植物花粉含量在該植被帶下降至最低,其中變化特征最明顯的是常綠櫟屬,平均僅為5.5%,其次是栲屬-柯屬,花粉含量僅為1.7%,木荷、山礬屬(Symplocos)、柳杉、柏科(Cupressaceae)等植被花粉幾乎消失;針葉樹花粉中以松屬(17.9%~95.1%,平均為30.1%)占優(yōu)勢(shì),含量最高,鐵杉含量下降至0.5%;灌木植物花粉含量(0.01%~19.0%,平均為3.6%)在該帶達(dá)到最高值,含量總體呈上升趨勢(shì),在該帶以薔薇科為主,平均為2.1%;草本植物花粉主要為禾本科(0.01%~39.1%,平均為16.7%),是武夷山脈禾本科植物最豐富的植被帶,菊科(2.1%)和毛莨科(Ranunculaceae,1.6%)均有增加,莎草科和蒿屬含量較低,平均為0.9%和0.6%;蕨類孢子在該帶占比較高,主要有芒萁屬,平均為14.4%,其次為三縫孢子(0.01%~37.1%,平均為10.1%)、單縫孢子(0.01%~8.8%,平均為3.1%)和石松(0.01%~12.5%,平均為2.5%)等。

    2.2 降趨對(duì)應(yīng)分析(DCA)與冗余分析(RDA)結(jié)果

    武夷山地區(qū)表土孢粉-氣候因子DCA排序結(jié)果顯示:第1 軸和第2 軸的特征值分別為0.369 和0.310,累計(jì)方差貢獻(xiàn)率為30.2%,而第3 軸和第4軸的特征值分別為0.158 和0.109,累計(jì)方差為11.9%,明顯低于前兩軸,由此表明表土孢粉種屬和樣品點(diǎn)的排序主要受控于第1 軸和第2 軸所代表的氣候因素。由于DCA 分析結(jié)果顯示,梯度長度值為2.247(<3.0),因此選擇冗余分析(RDA)方法進(jìn)一步探討表土孢粉數(shù)據(jù)和氣候因子之間的關(guān)系(表2)。

    表2 武夷山地區(qū)表土孢粉組合與氣候因子DCA分析結(jié)果Table 2 Results of DCA analysis between pollen assemblages and climatic factors in the Wuyi mountains

    為了進(jìn)一步分析氣候因子、植被與表土孢粉分布規(guī)律之間的關(guān)系,根據(jù)DCA 排序結(jié)果進(jìn)行冗余分析(RDA)。采用Monte Carlo置換檢驗(yàn)(n=499)評(píng)估備選環(huán)境因子的顯著性水平及其對(duì)響應(yīng)變量解釋(F值)的貢獻(xiàn),結(jié)果顯示,Mhco(F=0.668,P=0.662)顯著性水平>0.05,故除去該氣候變量。同時(shí)為了減少溫度、濕度和降水中各氣候因子之間的共線性和信息重疊,對(duì)氣候變量進(jìn)行3次方差膨脹因子檢驗(yàn)(Variance Inflation Factor, VIF)。檢驗(yàn)結(jié)果顯示,MAT (VIF=5.68)、MAH (VIF=4.34)、MAP (VIF=7.33)、Mpwa (VIF=5.32)和MPco(VIF=7.54)的VIF值較低,對(duì)武夷山脈表土孢粉組合的分布具有較強(qiáng)的影響,因此最終選擇這5個(gè)環(huán)境變量與研究區(qū)表土孢粉組合進(jìn)行分析。

    RDA 分析結(jié)果(表3)表明:第1 軸的特征值為0.129,物種與環(huán)境相關(guān)性為0.593;第2 軸的特征值為0.039,物種與環(huán)境相關(guān)性為0.457,前兩軸累計(jì)解釋了16.8%的物種信息和87.2%的物種-環(huán)境信息,第1 軸與第2 軸均表現(xiàn)出較高的物種與氣候因子的相關(guān)系數(shù),說明武夷山地區(qū)植被群落的動(dòng)態(tài)分布受到氣候因子的顯著影響。其中,Mpwa與第1軸的相關(guān)系數(shù)為0.53(0.90),MAT與第1軸的相關(guān)系數(shù)為0.53(0.89),MAH 與第1 軸的相關(guān)系數(shù)為0.45(0.75),表明這3 個(gè)氣候因子與第1 軸的相關(guān)性較高,而Mpco 與第2 軸的相關(guān)性較高,為0.34(0.74),MAP與第1軸和第2 軸的相關(guān)系數(shù)接近,分別為0.257(0.43)和0.258(0.56),與軸2 的相關(guān)性較強(qiáng)。

    表3 武夷山地區(qū)孢粉組合與氣候因子RDA分析結(jié)果Table 3 Results of RDA analysis between pollen assemblages and climatic factors in the Wuyi mountains

    從圖3 可以看出,孢粉帶IV、V 與其他孢粉帶中的樣品得到很好的區(qū)分,而孢粉帶I~I(xiàn)II 之間的樣品無法完全分開,3 個(gè)孢粉帶的成分較為相似,不存在明顯的界限。但RDA 排序圖也顯示,不同的植被帶樣品具有對(duì)應(yīng)的指示性花粉屬種,與孢粉帶I對(duì)應(yīng)的是栲屬-柯屬花粉,孢粉帶II、III的種類較多,其中與孢粉帶II對(duì)應(yīng)較好的是楓香樹屬、木荷屬等,松屬、落葉櫟屬與帶III有較好的對(duì)應(yīng),與帶IV 對(duì)應(yīng)的花粉類型主要有常綠櫟屬和鐵杉屬,孢粉帶V中的花粉種類不顯著,可能與本帶植被類型少、花粉含量低有關(guān)。

    圖3 武夷山地區(qū)主要孢粉種類和樣點(diǎn)的RDA排序及與相應(yīng)氣候因子的關(guān)系Fig.3 RDA ordination of pollen taxa and climate factors in the Wuyi mountains

    3 討論

    3.1 武夷山地區(qū)表土孢粉與現(xiàn)代植被相關(guān)關(guān)系

    武夷山地區(qū)表土孢粉組合分析結(jié)果顯示,武夷山地區(qū)孢粉種類以松屬、殼斗科的栲屬-柯屬、薔薇科、禾本科和蒿屬為優(yōu)勢(shì)類群,含較多常綠櫟屬、杉科、木荷屬、蕓香科和冬青屬等孢粉類型,同時(shí)也含有一定比例的落葉類成分,如楓香樹屬、落葉櫟屬、鵝耳櫪屬和漆樹科等,該孢粉組合與研究區(qū)亞熱帶常綠闊葉林植被的區(qū)域特征具有一致性。此外,孢粉組合中含有薔薇科、山礬屬、杜英屬等常見的喬灌木種類以及禾本科、蒿屬、莎草科、芒萁屬等草本類型。但一些在研究區(qū)常見的現(xiàn)代植物種類,如毛竹、赤楠、潤楠、樟樹、木蘭科等優(yōu)勢(shì)種類,在表土花粉中含量極低甚至沒有出現(xiàn)。主要原因可能是,竹子開花周期長而且花粉產(chǎn)量少,而樟科和木蘭科植物花粉壁極為薄弱不易保存,同時(shí)蟲媒植物花粉產(chǎn)量低,導(dǎo)致這些孢粉種類在表土孢粉譜中缺失。

    在武夷山低海拔地區(qū),花粉組合帶I 中闊葉類花粉(47.7%)占優(yōu)勢(shì),主要花粉類型有殼斗科的栲屬-柯屬、常綠櫟等,這與母體植被中殼斗科的栲屬、石櫟屬、青岡屬、樟科及山茶科等常綠闊葉樹分布范圍廣具有較好的對(duì)應(yīng)性。同時(shí),表土花粉譜中顯示孢粉帶I 松屬含量高,杉科含量也高于孢粉帶II,并含有一定量鵝耳櫟屬、楓香樹屬和落葉櫟屬等落葉類成分。而且,禾本科、芒萁屬在該帶有一定含量,最高含量達(dá)到24.1%,與山麓植被受人類活動(dòng)干擾后的表土孢粉組合特征相吻合。因此,孢粉帶I 花粉組合能較好地反映母體植被特征和優(yōu)勢(shì)種,在一定程度上可以反映人類活動(dòng)對(duì)植被群落的干擾作用。

    孢粉帶II 在常綠闊葉樹組成上與孢粉帶I 相似度較高,闊葉樹含量約為43.6%,花粉組合以栲屬-柯屬為優(yōu)勢(shì)成分,并伴有常綠櫟屬、木荷屬、蕓香科等常綠闊葉林群落,能較好地對(duì)應(yīng)武夷山地區(qū)常綠闊葉林中主要以毛椎、米櫧、毛錐、甜櫧等為優(yōu)勢(shì)種的特征。此外,花粉組合中以芒萁屬為主的蕨類孢子含量較高,與樣地草本層分布大量里白、芒萁相符合。聚類分析和RDA 分析結(jié)果表明,孢粉帶II與相鄰植被帶之間的界線不明顯。

    孢粉帶III 中喬木占優(yōu)勢(shì)(54.4%),針葉樹花粉與闊葉樹花粉含量分別為22.5%和31.9%,前者以松屬和杉科為主,后者主要有栲屬-柯屬、常綠櫟屬、落葉櫟屬、樺木屬、金縷梅科、榿木屬等,灌木花粉主要有薔薇科、野牡丹科、杜鵑花科等,草本植物中禾本科花粉含量也較高,這與植被群落以黃山松、南方鐵杉、甜櫧、木荷、山茶科、樟科等植被構(gòu)成的針闊混交林特點(diǎn)較為適應(yīng),且相較于孢粉帶II針葉樹的花粉含量(21.2%),孢粉帶III僅比前者高1.3%,但闊葉樹花粉在該帶的含量與孢粉帶II(43.6%)相比明顯降低,而陸生草本植物花粉也有顯著區(qū)別,在孢粉帶II 和III 中分別為4.9%和13.0%,相差8.1%。整體而言,孢粉帶III的花粉組合可以與其他植被帶的花粉組合相區(qū)別,能較好地反映植被群落特征。

    孢粉帶IV 中11 個(gè)樣品除了百祖山樣品外,均采于江西武夷山國家級(jí)自然保護(hù)區(qū),其中南方鐵杉天然林主要分布于海拔在1 400~2 000 m。花粉組合中針葉樹花粉含量(22.9%)相對(duì)于孢粉帶III略升高了0.5%,主要為松屬(19.7%)和鐵杉屬(2.6%),能部分反映母體植被中鐵杉的分布,但含量較低。灌木花粉含量較低(3.6%),主要有薔薇科和杜鵑花科,草本植物花粉中禾本科含量也僅為1.9%,相對(duì)于孢粉帶III,帶IV的灌木花粉(4.1%)和禾本科(9.3%)均有下降,此外有一定含量的栲屬-柯屬,并伴生落葉櫟屬、楓香樹屬、水青岡屬、鵝耳櫪屬等落葉成分,這種花粉組合在一定程度上反映針葉林母體植被中存在薔薇科、禾本科、山茶科、杜鵑花科及殼斗科等植被,但含量較低。孢粉帶IV含有松屬和鐵杉屬等種類,但其作為建群種與花粉高含量屬種的對(duì)應(yīng)性較差,故花粉組合不能代表針葉林植被群落的組成特征。

    孢粉帶V中闊葉樹花粉含量驟減,花粉組合以針葉樹(31.6%)和蕨類孢子(30.4%)占優(yōu)勢(shì),松屬花粉在黃崗山樣品達(dá)到峰值(95.1%),且該樣點(diǎn)中孢粉種類較少,僅出現(xiàn)5個(gè)孢粉種類,可能與取樣點(diǎn)植被稀疏有關(guān)。灌木和草本植物花粉總含量為26.9%,其中灌木植物花粉僅4.4%,草本植物花粉主要有禾本科(16.7%)。但該孢粉帶中,其他樣點(diǎn)孢粉譜中草本類花粉含量可達(dá)29.8%,其中以禾本科(22.2%)占絕對(duì)優(yōu)勢(shì),但灌木植物花粉(5.8%)仍較低。因此,該孢粉帶的花粉組合與灌叢草甸植被群落中主要分布青茅、白檀、云錦杜鵑優(yōu)勢(shì)植物種類有較好的對(duì)應(yīng)性。

    綜上,不同植物帶的孢粉組合具有不同的代表性成分和變化特征,通過優(yōu)勢(shì)種屬的含量可以判別不同海拔高度的植被。且聚類分析和RDA 分析結(jié)果也表明,多數(shù)植被帶可以得到區(qū)分。孢粉帶I與II與以常綠闊葉林為主的植被類型較為對(duì)應(yīng),優(yōu)勢(shì)屬種為栲屬-柯屬;孢粉帶III則以針葉樹、闊葉樹樹種為主,與針闊混交林相適應(yīng),但針葉類花粉含量與孢粉帶II相當(dāng),指示性花粉為松屬;孢粉帶IV以常綠櫟屬和鐵杉屬含量升高為特點(diǎn),與其他植被帶相區(qū)別;孢粉帶V則發(fā)育大量禾本科,對(duì)應(yīng)以青茅為主的草甸植被類型。

    總之,現(xiàn)代孢粉組合與RDA 排序結(jié)果表明,武夷山脈地區(qū)以殼斗科柯屬-栲屬、松屬為優(yōu)勢(shì)群,并含一定量常綠櫟屬、杉科、木荷屬等常綠闊葉林成分,表土孢粉組合與現(xiàn)代植被在建群、優(yōu)勢(shì)科屬上基本一致,并且不同植物帶的孢粉組合具有不同的指示性成分和變化特征,孢粉組合能較好地反映研究區(qū)對(duì)應(yīng)的現(xiàn)代植被群系及其垂直分布特點(diǎn)。

    3.2 主要花粉類型散布特征分析

    3.2.1 喬、灌木類花粉散布特征 殼斗科栲屬-柯屬是武夷山地區(qū)出現(xiàn)頻率最高的喬木花粉種屬,花粉總含量(0.2%~91.4%)也相對(duì)最高,平均為23.2%,且花粉含量變化與植被分布情況較為一致,表現(xiàn)為隨海拔升高整體呈減少趨勢(shì),高含量集中分布在以米櫧、毛錐、甜櫧等樹種為優(yōu)勢(shì)的亞熱帶常綠闊葉林群落中。在所有樣品中,栲屬-柯屬花粉含量明顯高于其植被覆蓋度,說明其具有超代表性,在亞熱帶地區(qū)其他地方表土孢粉散布研究結(jié)果(李永飛 等,2018;全曉文 等,2022)也有相似性。

    常綠櫟屬在武夷山脈表土花粉樣品中頗為常見,但含量變化大,最高為85.0%,平均為8.9%,高含量集中在海拔1 400~1 800 m 的孢粉帶IV 中,平均為31.0%,其他植被群落中平均含量低于10%。由此表明武夷山地區(qū)常綠櫟花粉具有超代表性。該結(jié)論與已有研究一致:同樣位于亞熱帶地區(qū)的云南西部、瀘沽湖、老君山和香格里拉等區(qū)域的花粉代表性問題研究表明,常綠櫟屬是超代表性花粉類型,但母體植被附近仍是花粉含量最高值的主要分布點(diǎn)(Shen et al., 2013;張華 等,2016;李永飛等,2018;譚金鳳 等,2020)。而郝秀東等(2020)在浙江西白山的研究發(fā)現(xiàn),常綠櫟屬具有中等代表性,其含量在香榧林、針闊混交林、茶園和竹林中分別為6.8%~6.9%、25.5%~26.1%、4.4%~4.5%和16.1%~16.5%,其代表性適中的原因可能是香榧(Torreya)、茶樹(Camellia)植被的大量分布致使常綠櫟生長受限。

    松屬植物廣泛分布于武夷山地區(qū),主要包括馬尾松和黃山松2個(gè)物種,其中1 100 m以下主要為人工種植的馬尾松林,而黃山松主要分布在1 100 m以上地區(qū)。本研究表土樣品的孢粉組合中均發(fā)現(xiàn)有較多松屬花粉,含量變化范圍為1.1%~95.1%,平均為19.9%,僅次于柯-栲屬花粉。其中含量>30%的樣品共27個(gè),主要集中在海拔587~1 840 m旱地、疏林地內(nèi),含量在10%~30%的樣品有29 個(gè),<10%的樣品有20個(gè),大多樣點(diǎn)采自少見或未見松屬植被分布的有林地,表明松屬花具有超代表性,花粉含量明顯高于松屬植被在群落中的比例。這可能與松屬花粉具有較發(fā)達(dá)的氣囊從而傳播能力強(qiáng)有關(guān)(唐領(lǐng)余 等,2016),且松屬花粉產(chǎn)量大,已有研究大多也揭示松屬花粉具有超代表性(李文漪 等,1990;Li et al., 2000;許清海 等,2007;李圓圓等,2016)。

    鐵杉植物在中國僅分布于橫斷山區(qū)、秦嶺—大巴山區(qū)、華東山地以及臺(tái)灣中央山脈(楊青松 等,2009)。在江西武夷山國家級(jí)自然保護(hù)區(qū)內(nèi),海拔1 400~2 000 m 范圍內(nèi)完整保存著面積達(dá)1 560 hm2的南方鐵杉天然原始森林,本文樣點(diǎn)采集于江西武夷山國家級(jí)自然保護(hù)區(qū)海拔1 288 m~2 158 m 范圍內(nèi),花粉含量范圍為0.1~7.9%,平均為2.5%,其中含量>1%的樣品集中于海拔1 422~1 843 m 以南方鐵杉、黃山松和杉木等針葉樹構(gòu)成的針葉林植被群落中,采樣點(diǎn)附近有鐵杉生長。但其他表土孢粉樣點(diǎn)幾乎不含鐵杉屬花粉。由此可見,鐵杉花粉主要分布在母體植被周圍,但花粉產(chǎn)量較低,代表性相對(duì)較低,與肖霞云等(2009)在云南玉龍雪山的表土花粉研究一致。關(guān)于鐵杉屬花粉的代表性,有學(xué)者認(rèn)為其具有超代表性,如在太白山、云南麗江-大理地區(qū)表土采樣點(diǎn)附近缺乏鐵杉植物群落的分布,但表土孢粉組合中均發(fā)現(xiàn)其花粉(李春海 等,2008;李永飛 等,2017)。

    3.2.2 草本類花粉的散布特征 在武夷山表土樣品花粉組合中,禾本科的平均含量為5.1%,但在武夷山脈高海拔地區(qū)表土樣品中禾本科花粉平均含量約為16.7%。在武夷山海拔~1 370 m 形成以青茅為建群種的中山山頂灌草叢植物群落,禾本科植物的含量高達(dá)54.6%,因此禾本科花粉的高含量與表土樣點(diǎn)附近出現(xiàn)灌草叢有關(guān),采樣點(diǎn)附近喬木層植被多樣性降低且分布稀疏,植被蓋度顯著降低,而禾本科雜草生長茂盛且花粉主要散布在母體植物周圍。在其他植被類型中,大部分表土孢粉樣品的禾本科花粉含量大多低于5%。此外,在武夷山地區(qū)大約海拔~800 m分布著大面積毛竹林,但這些樣點(diǎn)中禾本科花粉含量較低,僅為2.5%,可能與竹子開花周期長、花粉不易保存有關(guān),導(dǎo)致表土禾本科花粉含量相對(duì)較低。本研究中,禾本科花粉含量整體偏低表明該草本類花粉呈低代表性特征,大量表土花粉研究也證實(shí)禾本科植物花粉呈低代表性是一種普遍現(xiàn)象(Zheng et al., 2008;魏海成 等,2009;李圓圓 等,2016;趙楠楠 等,2019;郝秀東 等,2020)。但在青藏高原地區(qū),表土孢粉中禾本科花粉含量卻明顯高于植被蓋度,呈現(xiàn)超代表性(龐有智 等,2016),因此在不同植物群落中,禾本科植物花粉的代表性可能存在較復(fù)雜的情況(張文超等,2013)。

    蒿屬花粉在武夷山地區(qū)表土樣品中的平均含量為1.1%,但僅有1 個(gè)樣品的含量達(dá)到20.6%,而大多樣品的蒿屬花粉含量較低,與北方地區(qū)蒿屬花粉的分布特征存在一定的差異(尚雪 等,2009;張文超 等,2013;趙楠楠 等,2019)。在中國北方地區(qū),蒿屬植物花粉產(chǎn)量高、花粉較輕、易于風(fēng)力傳播而具有超代表性,因此在北方草原或者荒漠地區(qū),如果表土孢粉組合中蒿屬花粉含量較低時(shí),則不能指示周圍存在蒿屬植物,可能指示外來花粉(楊振京 等,2011)。因此,武夷山脈地區(qū)植被中,以亞熱帶常綠喬木為主,蒿屬植物相對(duì)較少,表土孢粉蒿屬植物花粉含量整體相對(duì)較低。

    3.3 人類活動(dòng)對(duì)植被的影響

    華南亞熱帶東部地區(qū)低海拔地區(qū)植被受到強(qiáng)烈人類活動(dòng)的干擾,地帶性植被受到不同程度的影響,尤其是低海拔丘陵地區(qū)在大規(guī)模人類活動(dòng)的長期作用下生境遭到嚴(yán)重破壞,人工樹種或次生樹種成為取代原始樹種的主要植被類型(丁偉 等,2011)。本研究中,武夷山山麓低海拔地區(qū)主要分布成片的茶園、毛竹林、馬尾松林和杉木林等人工植被或更新次生植被,現(xiàn)代花粉組合受到這些人工擾動(dòng)次生植被和農(nóng)作物的共同影響,主要表現(xiàn)出常綠櫟屬花粉含量極低,同時(shí)花粉種類多樣性也較低;而蕨類孢子的芒萁屬和三縫孢子含量相對(duì)較高,該孢粉組合難以反映亞熱帶常綠闊葉林的地帶性植被特征。此外,松屬花粉在低海拔地區(qū)也有較高含量(主要是人工種植的針葉類,并不是揭示自然植被演替過程中的針葉林),這與南嶺低海拔地區(qū)的花粉組合類似,即低海拔地區(qū)植被受到人類活動(dòng)的強(qiáng)烈干擾,導(dǎo)致次生林下芒箕類孢子含量相對(duì)較高(鄭卓 等,2004;全曉文 等,2022)。

    3.4 表土孢粉與環(huán)境因子的關(guān)系

    氣候因子是影響和控制植被組成及分布的主要因素之一,同時(shí)與地形、水土、海拔、土壤類型和坡度等小環(huán)境因素有關(guān)。在國內(nèi),前人已經(jīng)在大空間尺度上探討植物地理分布與孢粉散布規(guī)律及氣候因子之間的定量關(guān)系(鄭卓 等,2008),發(fā)現(xiàn)絕大多數(shù)表土孢粉含量的氣候指示性與植物的地理分布較為一致;但超代表性的孢粉種類或者R值偏高的一些孢粉種類在氣候區(qū)間上通常比植物種類分布的氣候范圍偏大,如松屬和草本類;有部分孢粉種類指示的氣候相對(duì)植物的氣候范圍存在一定的偏移,如杜鵑花科、菊科、禾本科和藜科(Chenopodiaceae)等,這些孢粉種類與植物種類在氣候特征上存在差異性,可能與孢粉在顯微鏡下分類鑒定水平和植物種類難以相對(duì)應(yīng)有關(guān);此外,由于自然植被受到人類活動(dòng)的強(qiáng)烈干擾,導(dǎo)致表土孢粉-植物-氣候之間的對(duì)應(yīng)性也存在一定的偏差(Zheng et al.,2008)。

    在區(qū)域尺度上,前人已經(jīng)開展大量關(guān)于表土孢粉散布規(guī)律與植被的相互關(guān)系研究(舒軍武 等,2010;Wei et al., 2011;李曼玥 等,2012;Fang et al., 2015; Wei et al., 2016;李圓圓 等,2016;Li et al., 2017;余少華 等,2017;吳慧賢 等,2018;李怡雯 等,2019;Han et al., 2020;郝秀東 等,2020;汪啟容 等,2021;全曉文 等,2022),主要集中在青藏高原和中國北方地區(qū)。青藏高原地區(qū)表土孢粉散布規(guī)律與氣候因子的相關(guān)性研究表明:孢粉組合主要與年降水量、年均相對(duì)濕度、年蒸發(fā)量和干燥度4個(gè)氣象因子相關(guān)性最高,揭示了降水量可能是影響干旱半干旱及高寒地區(qū)植被生長的最主要因素(王芳芳 等,2017)。華南亞熱帶南嶺地區(qū)表土孢粉研究表明,其在空間上含量的變化與年均降雨量、最暖月降水和最冷月降水有關(guān)(全曉文等,2022)。而在本研究中,武夷山地區(qū)的孢粉組合與區(qū)域溫度、降水量及海拔高度密切相關(guān)。RDA排序分析結(jié)果進(jìn)一步揭示了現(xiàn)代花粉與氣候因子之間的關(guān)系,武夷山地區(qū)花粉組合散布受控的氣候因子包括:MAT、MAP、Mpco、Mpwa 和MAH,表明武夷山地區(qū)表土孢粉種類散布規(guī)律與氣候因子之間的相關(guān)性除了與溫度和降水量有關(guān),可能還受空氣濕度控制。武夷山位于華南亞熱帶東亞夏季風(fēng)活動(dòng)的主要影響區(qū)域,表土孢粉組合與氣候因子相關(guān)性分析顯示,溫度和降水是控制武夷山地區(qū)植被孢粉分布的主要?dú)夂蛞蛩兀瑫r(shí)在垂直植被帶上孢粉散布與溫度變化緊密相關(guān)(見圖3),本研究與南嶺地區(qū)表土孢粉散布規(guī)律和氣候因子的關(guān)系相一致(李圓圓 等,2016;全曉文 等,2022)。冗余分析的橫坐標(biāo)和縱坐標(biāo)可能反映孢粉種屬對(duì)最熱月溫度和年均降雨量變化的響應(yīng)。

    4 結(jié)論與討論

    4.1 結(jié)論

    通過分析武夷山地區(qū)不同植被類型下采集的76個(gè)表土樣品孢粉組合特征,并探討其與植被群落、氣候因子之間的關(guān)系,得出以下主要結(jié)論:

    1)表土孢粉組合與現(xiàn)代植被在建群、優(yōu)勢(shì)科屬上基本一致,表現(xiàn)出明顯的亞熱帶特征,主要有松屬、殼斗科柯屬-栲屬、常綠櫟屬、杉科、木荷屬、蕓香科、冬青屬、薔薇科、禾本科、蒿屬等花粉類型,基本能反映研究區(qū)現(xiàn)生植被狀況;

    2)武夷山地區(qū)不同垂直植被帶的表土孢粉組合具有不同的指示性成分和變化特征,孢粉組合能較好地反映研究區(qū)內(nèi)對(duì)應(yīng)的現(xiàn)代植被群系及其垂直分布特點(diǎn);

    3)RDA 分析結(jié)果表明,MAT、MAH、MAP、Mpco和Mpwa是控制武夷山地區(qū)現(xiàn)代孢粉散布特征的主要?dú)夂蛞蜃?,在垂直植被帶上孢粉散布與溫度變化緊密相關(guān);

    4)花粉的代表性分析顯示,栲屬-柯屬、松屬、常綠櫟屬和蒿屬均具有超代表性,鐵杉屬和禾本科代表性較低;

    5)武夷山位于中國亞熱帶東部地區(qū),人類活動(dòng)強(qiáng)烈,現(xiàn)代孢粉譜主要表現(xiàn)在低海拔地區(qū)樣品中具有相對(duì)較高含量的禾本科、蒿屬、松屬等花粉,且花粉種類較少。因此,利用化石孢粉恢復(fù)中國亞熱帶東部地區(qū)的古環(huán)境和古氣候時(shí)不能無條件地參考現(xiàn)代植被;

    6)在亞熱帶地區(qū),由于孢粉產(chǎn)量、傳播方式、保存特性以及人類活動(dòng)的干擾,導(dǎo)致表土孢粉在組成、科屬類別上與研究區(qū)現(xiàn)生植被存在一定差異性。

    通過分析武夷山脈76個(gè)不同海拔的表土孢粉樣品,探討了該地區(qū)表土孢粉的垂直分布規(guī)律及其與母體植被之間的關(guān)系,本文可為后續(xù)在武夷山及周邊地區(qū)準(zhǔn)確利用化石孢粉重建古環(huán)境、古氣候提供基礎(chǔ)資料。當(dāng)然,本研究仍存在一些不足,首先,部分孢粉在光學(xué)顯微鏡下的鑒定精度只能達(dá)到科屬一級(jí),而同一科屬植物的種間地理空間分布的差異在一定程度上導(dǎo)致孢粉散布與氣候之間的相互關(guān)系產(chǎn)生一定的偏差;其次,現(xiàn)代表土孢粉-植物-氣候之間的相互關(guān)系仍然存在較復(fù)雜的關(guān)系,在未來仍然需要開展更多的研究。

    猜你喜歡
    孢粉禾本科武夷山
    東北地區(qū)末次盛冰期以來孢粉數(shù)據(jù)庫構(gòu)建技術(shù)與應(yīng)用
    它們的家園——武夷山
    大美武夷山
    孢粉在流域生態(tài)環(huán)境中的指示作用及研究進(jìn)展
    人民珠江(2019年10期)2019-11-11 03:22:22
    西藏仲巴縣隆格爾地區(qū)漸新統(tǒng)日貢拉組孢粉組合的發(fā)現(xiàn)及其意義
    風(fēng)雨廊橋——武夷山馀慶橋
    紅土地(2018年12期)2018-04-29 09:16:50
    武夷山
    上海師范大學(xué)揭示花粉壁孢粉素合成的調(diào)控機(jī)制
    蔬菜(2018年8期)2018-01-16 15:39:27
    禾本科植物遺傳中根癌農(nóng)桿菌介導(dǎo)法的應(yīng)用
    禾本科植物之王
    久久久久亚洲av毛片大全| 99re在线观看精品视频| 黑人操中国人逼视频| 国产欧美日韩一区二区精品| 精品久久久久久,| 亚洲天堂国产精品一区在线| or卡值多少钱| 久久人人精品亚洲av| 一边摸一边抽搐一进一出视频| 午夜老司机福利片| 久久中文字幕一级| 亚洲av成人不卡在线观看播放网| 多毛熟女@视频| 91av网站免费观看| 国产蜜桃级精品一区二区三区| 无限看片的www在线观看| 国产97色在线日韩免费| 黄色视频,在线免费观看| 国产一区在线观看成人免费| 国产一级毛片七仙女欲春2 | a在线观看视频网站| 美女午夜性视频免费| 男人的好看免费观看在线视频 | 久久久久久久精品吃奶| 97人妻精品一区二区三区麻豆 | 99精品久久久久人妻精品| 国产av精品麻豆| 国产精品久久久人人做人人爽| 久久久久久免费高清国产稀缺| 91在线观看av| 九色国产91popny在线| e午夜精品久久久久久久| 99香蕉大伊视频| 久久久久九九精品影院| 亚洲欧美一区二区三区黑人| 国产精品,欧美在线| 成人三级黄色视频| 久久精品91蜜桃| 97人妻天天添夜夜摸| 无限看片的www在线观看| 亚洲国产欧美网| av天堂久久9| 极品人妻少妇av视频| 麻豆国产av国片精品| 又大又爽又粗| 国产一区二区三区在线臀色熟女| 国产精品久久久久久精品电影 | 制服诱惑二区| 国产高清有码在线观看视频 | 美女大奶头视频| 手机成人av网站| 男男h啪啪无遮挡| 一本久久中文字幕| 俄罗斯特黄特色一大片| 久久人人97超碰香蕉20202| 一边摸一边做爽爽视频免费| 亚洲欧美激情在线| 又大又爽又粗| 久久人妻福利社区极品人妻图片| 99国产精品免费福利视频| 美女高潮喷水抽搐中文字幕| 人人妻人人爽人人添夜夜欢视频| 久久精品国产清高在天天线| 午夜免费鲁丝| 美女扒开内裤让男人捅视频| 亚洲狠狠婷婷综合久久图片| 1024香蕉在线观看| 国产不卡一卡二| 一进一出抽搐gif免费好疼| 日韩中文字幕欧美一区二区| 国产精品98久久久久久宅男小说| 国产1区2区3区精品| 看片在线看免费视频| 9191精品国产免费久久| 午夜久久久久精精品| 久久伊人香网站| 18禁裸乳无遮挡免费网站照片 | 大型黄色视频在线免费观看| 午夜福利高清视频| 在线观看66精品国产| 亚洲国产精品合色在线| 日韩精品中文字幕看吧| 老司机深夜福利视频在线观看| 午夜日韩欧美国产| 国产精品电影一区二区三区| 中文字幕色久视频| 欧美丝袜亚洲另类 | 成年版毛片免费区| 久久精品国产综合久久久| 人人澡人人妻人| 国产精品野战在线观看| 欧美一级a爱片免费观看看 | 亚洲国产毛片av蜜桃av| 日韩精品青青久久久久久| 一进一出抽搐gif免费好疼| 国产成人欧美| 亚洲电影在线观看av| 老熟妇乱子伦视频在线观看| 精品国产亚洲在线| 悠悠久久av| 国产精品免费视频内射| av免费在线观看网站| 欧美绝顶高潮抽搐喷水| 国产伦一二天堂av在线观看| 精品一区二区三区av网在线观看| 色尼玛亚洲综合影院| 老司机福利观看| 亚洲第一电影网av| 成人三级做爰电影| 欧美日韩福利视频一区二区| 亚洲五月婷婷丁香| 亚洲五月天丁香| 免费搜索国产男女视频| 中亚洲国语对白在线视频| 99riav亚洲国产免费| 国产成年人精品一区二区| 男女下面进入的视频免费午夜 | 亚洲色图 男人天堂 中文字幕| 国产精品精品国产色婷婷| 91麻豆精品激情在线观看国产| 日韩成人在线观看一区二区三区| 成人三级黄色视频| 亚洲欧美日韩无卡精品| 亚洲国产精品999在线| 自线自在国产av| 国产一区二区三区视频了| 黄网站色视频无遮挡免费观看| 激情在线观看视频在线高清| 国产精品 国内视频| 久久久久久久精品吃奶| 啪啪无遮挡十八禁网站| 99精品欧美一区二区三区四区| av免费在线观看网站| 91九色精品人成在线观看| av福利片在线| 亚洲一卡2卡3卡4卡5卡精品中文| 日本一区二区免费在线视频| 久久久久九九精品影院| 欧美色视频一区免费| 久久天躁狠狠躁夜夜2o2o| 亚洲五月天丁香| 青草久久国产| 国产精品秋霞免费鲁丝片| 免费在线观看亚洲国产| 自拍欧美九色日韩亚洲蝌蚪91| 中出人妻视频一区二区| 老司机午夜福利在线观看视频| 亚洲精品国产一区二区精华液| 美女扒开内裤让男人捅视频| 久久热在线av| 久久午夜综合久久蜜桃| 国产精品一区二区在线不卡| 午夜成年电影在线免费观看| 久久伊人香网站| 国产日韩一区二区三区精品不卡| 中文字幕久久专区| 可以在线观看毛片的网站| 中文字幕久久专区| 久久精品国产综合久久久| 国产成+人综合+亚洲专区| 亚洲av成人一区二区三| 久久青草综合色| 中文字幕高清在线视频| 国产av精品麻豆| 色播亚洲综合网| 桃红色精品国产亚洲av| 一卡2卡三卡四卡精品乱码亚洲| 热re99久久国产66热| 国产99白浆流出| av视频免费观看在线观看| 亚洲av第一区精品v没综合| 黄色片一级片一级黄色片| 亚洲 欧美一区二区三区| 日本黄色视频三级网站网址| 久久婷婷人人爽人人干人人爱 | 亚洲精品av麻豆狂野| 亚洲精品中文字幕在线视频| 精品国产美女av久久久久小说| 国产精品久久久久久人妻精品电影| 久久精品国产99精品国产亚洲性色 | 国产色视频综合| 国产国语露脸激情在线看| 欧美日本视频| 免费高清视频大片| 日本在线视频免费播放| 老汉色∧v一级毛片| 欧美日韩乱码在线| 免费一级毛片在线播放高清视频 | 国产成人啪精品午夜网站| 黄色 视频免费看| 俄罗斯特黄特色一大片| 欧美最黄视频在线播放免费| 国产免费av片在线观看野外av| 一进一出抽搐动态| 欧美日韩精品网址| 亚洲精品美女久久av网站| 女人被狂操c到高潮| 日本黄色视频三级网站网址| 亚洲九九香蕉| 国产日韩一区二区三区精品不卡| 色在线成人网| 国产熟女午夜一区二区三区| 国产精品99久久99久久久不卡| 亚洲成av片中文字幕在线观看| 精品国产一区二区久久| 亚洲熟女毛片儿| 嫩草影院精品99| 亚洲avbb在线观看| 涩涩av久久男人的天堂| 好男人在线观看高清免费视频 | 男女午夜视频在线观看| 亚洲色图 男人天堂 中文字幕| 欧美乱色亚洲激情| 中文字幕人妻丝袜一区二区| 丰满人妻熟妇乱又伦精品不卡| 97碰自拍视频| 母亲3免费完整高清在线观看| 制服人妻中文乱码| 久久国产精品影院| 淫妇啪啪啪对白视频| 一二三四在线观看免费中文在| 妹子高潮喷水视频| 亚洲男人天堂网一区| 中文字幕另类日韩欧美亚洲嫩草| 国产极品粉嫩免费观看在线| 亚洲在线自拍视频| 999久久久精品免费观看国产| 天堂√8在线中文| 真人一进一出gif抽搐免费| 窝窝影院91人妻| 18禁观看日本| 亚洲人成77777在线视频| 非洲黑人性xxxx精品又粗又长| 波多野结衣巨乳人妻| 黑人巨大精品欧美一区二区蜜桃| 成人免费观看视频高清| 日本vs欧美在线观看视频| 美国免费a级毛片| 免费少妇av软件| 精品国产乱码久久久久久男人| 国产成人精品久久二区二区91| 日韩有码中文字幕| 国产精品电影一区二区三区| 好男人电影高清在线观看| 亚洲中文av在线| av片东京热男人的天堂| 精品国产乱子伦一区二区三区| 欧美日韩黄片免| 在线免费观看的www视频| 久久精品成人免费网站| 91大片在线观看| 国产一级毛片七仙女欲春2 | 国产一区二区在线av高清观看| 久9热在线精品视频| 国产精品电影一区二区三区| 老熟妇仑乱视频hdxx| 90打野战视频偷拍视频| 两人在一起打扑克的视频| 国产熟女xx| 丁香欧美五月| 人成视频在线观看免费观看| 日本三级黄在线观看| 动漫黄色视频在线观看| 在线观看免费视频网站a站| 99国产精品一区二区蜜桃av| 国产aⅴ精品一区二区三区波| 国产精品亚洲一级av第二区| 国产一区二区激情短视频| 这个男人来自地球电影免费观看| 精品一区二区三区av网在线观看| 神马国产精品三级电影在线观看 | 精品福利观看| 一区福利在线观看| 18禁裸乳无遮挡免费网站照片 | 亚洲第一欧美日韩一区二区三区| 天天一区二区日本电影三级 | 欧洲精品卡2卡3卡4卡5卡区| 一级a爱视频在线免费观看| 少妇被粗大的猛进出69影院| 青草久久国产| 一区二区三区国产精品乱码| 美女免费视频网站| 大型av网站在线播放| 成人亚洲精品av一区二区| 国产亚洲精品久久久久久毛片| www.999成人在线观看| 精品无人区乱码1区二区| 老熟妇乱子伦视频在线观看| 亚洲精品粉嫩美女一区| 人人妻人人澡人人看| 在线观看舔阴道视频| 欧美日本视频| 日本 av在线| 1024视频免费在线观看| 国产高清videossex| 午夜精品国产一区二区电影| 国产亚洲精品久久久久5区| 老司机福利观看| 国产精品亚洲一级av第二区| 人人妻人人爽人人添夜夜欢视频| 美女 人体艺术 gogo| 手机成人av网站| 亚洲第一青青草原| 真人一进一出gif抽搐免费| 午夜日韩欧美国产| 国产av在哪里看| 久久久久久久久久久久大奶| 久久久久久久久免费视频了| 一进一出抽搐gif免费好疼| 日本免费a在线| 国产高清videossex| 欧美一区二区精品小视频在线| 久久中文字幕一级| 国产精品永久免费网站| 激情视频va一区二区三区| 国产三级黄色录像| 亚洲av电影不卡..在线观看| 欧美亚洲日本最大视频资源| 亚洲狠狠婷婷综合久久图片| 正在播放国产对白刺激| 在线观看www视频免费| 变态另类丝袜制服| 亚洲成国产人片在线观看| 变态另类成人亚洲欧美熟女 | 人人澡人人妻人| 久久狼人影院| 久久久精品国产亚洲av高清涩受| 欧美性长视频在线观看| 在线永久观看黄色视频| 女警被强在线播放| 亚洲国产精品合色在线| 国产欧美日韩精品亚洲av| 国产麻豆69| 老熟妇乱子伦视频在线观看| www.999成人在线观看| 老司机午夜福利在线观看视频| 九色亚洲精品在线播放| 久久久水蜜桃国产精品网| 桃色一区二区三区在线观看| 欧美日韩瑟瑟在线播放| 99国产精品免费福利视频| 亚洲成国产人片在线观看| 在线观看www视频免费| 嫩草影院精品99| 亚洲精品国产一区二区精华液| 黄色 视频免费看| 91av网站免费观看| 黄色丝袜av网址大全| 久久久久久人人人人人| 性色av乱码一区二区三区2| 精品久久久久久久久久免费视频| 婷婷精品国产亚洲av在线| cao死你这个sao货| 亚洲情色 制服丝袜| 亚洲美女黄片视频| 视频在线观看一区二区三区| 成人精品一区二区免费| 日韩欧美一区视频在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲av五月六月丁香网| 中文字幕人妻丝袜一区二区| 久久久国产成人精品二区| 国产午夜福利久久久久久| 丝袜美腿诱惑在线| 两个人看的免费小视频| 在线观看一区二区三区| 欧美成狂野欧美在线观看| 曰老女人黄片| 欧美日韩亚洲综合一区二区三区_| 丝袜美足系列| 国产成人精品在线电影| 一本大道久久a久久精品| 亚洲av电影在线进入| 看片在线看免费视频| 久久精品影院6| 9色porny在线观看| 亚洲av成人av| 看黄色毛片网站| 亚洲精品国产精品久久久不卡| 国产激情欧美一区二区| 欧美激情极品国产一区二区三区| 国产成人av激情在线播放| 国产精品,欧美在线| 日韩欧美在线二视频| 日韩 欧美 亚洲 中文字幕| 免费在线观看亚洲国产| 午夜视频精品福利| 国产欧美日韩一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 国产熟女午夜一区二区三区| 亚洲 国产 在线| 美女 人体艺术 gogo| 精品人妻1区二区| 久久久精品国产亚洲av高清涩受| 免费看a级黄色片| 性欧美人与动物交配| 亚洲欧美日韩无卡精品| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美一区二区三区在线观看| svipshipincom国产片| av欧美777| 伦理电影免费视频| 国产精品免费一区二区三区在线| av天堂久久9| 免费av毛片视频| 国产精品永久免费网站| 亚洲男人天堂网一区| xxx96com| 91精品三级在线观看| 久久热在线av| 亚洲视频免费观看视频| 最新美女视频免费是黄的| 99热只有精品国产| 国产激情欧美一区二区| 精品久久久久久久毛片微露脸| 美国免费a级毛片| 久久中文字幕一级| 色老头精品视频在线观看| a级毛片在线看网站| 97人妻天天添夜夜摸| 亚洲人成电影观看| 动漫黄色视频在线观看| 麻豆国产av国片精品| 纯流量卡能插随身wifi吗| 无人区码免费观看不卡| 日本五十路高清| 欧美中文日本在线观看视频| 天天躁狠狠躁夜夜躁狠狠躁| 大码成人一级视频| 日韩精品免费视频一区二区三区| 欧美中文综合在线视频| 中文字幕人妻丝袜一区二区| 日日爽夜夜爽网站| 极品教师在线免费播放| www日本在线高清视频| 天天添夜夜摸| 老司机在亚洲福利影院| 日韩高清综合在线| 少妇粗大呻吟视频| 97超级碰碰碰精品色视频在线观看| 免费看a级黄色片| 精品欧美一区二区三区在线| 欧美老熟妇乱子伦牲交| 欧美最黄视频在线播放免费| 欧美激情高清一区二区三区| 成人国产综合亚洲| 亚洲色图 男人天堂 中文字幕| 变态另类成人亚洲欧美熟女 | 国产欧美日韩一区二区三区在线| 国产1区2区3区精品| 亚洲精品美女久久av网站| 人妻丰满熟妇av一区二区三区| 久久精品影院6| 久久精品国产综合久久久| 久久中文字幕一级| 欧美乱色亚洲激情| 一级片免费观看大全| 午夜激情av网站| 免费看十八禁软件| 午夜福利高清视频| 国产区一区二久久| 美女大奶头视频| 一卡2卡三卡四卡精品乱码亚洲| 婷婷六月久久综合丁香| 久久人妻av系列| 欧美黄色片欧美黄色片| 日日爽夜夜爽网站| 狠狠狠狠99中文字幕| 国产精品99久久99久久久不卡| 久久午夜亚洲精品久久| 国产精品影院久久| 欧美另类亚洲清纯唯美| 宅男免费午夜| 一级毛片精品| av视频在线观看入口| 欧美不卡视频在线免费观看 | 免费女性裸体啪啪无遮挡网站| 又黄又爽又免费观看的视频| 午夜精品在线福利| av中文乱码字幕在线| 欧美黄色淫秽网站| 日韩欧美一区视频在线观看| 人人澡人人妻人| 国产欧美日韩精品亚洲av| av电影中文网址| 中文字幕另类日韩欧美亚洲嫩草| 亚洲专区国产一区二区| 18禁黄网站禁片午夜丰满| 国产亚洲精品第一综合不卡| 午夜久久久在线观看| 国产成年人精品一区二区| 久久久久久久久中文| 亚洲aⅴ乱码一区二区在线播放 | 变态另类成人亚洲欧美熟女 | 亚洲国产精品sss在线观看| 国产精品 欧美亚洲| 午夜免费观看网址| 9191精品国产免费久久| 99国产精品一区二区蜜桃av| 亚洲中文字幕一区二区三区有码在线看 | 91九色精品人成在线观看| 日本免费a在线| av免费在线观看网站| 久久久国产成人免费| 亚洲国产精品sss在线观看| 亚洲国产中文字幕在线视频| 黄色视频不卡| 欧美乱色亚洲激情| 十八禁人妻一区二区| 老司机靠b影院| 麻豆国产av国片精品| 午夜福利18| 亚洲第一电影网av| 美女高潮喷水抽搐中文字幕| 久久精品亚洲精品国产色婷小说| 丝袜人妻中文字幕| 国产精品久久久久久亚洲av鲁大| 精品久久久久久成人av| 在线免费观看的www视频| 久久精品国产综合久久久| 99国产精品一区二区蜜桃av| 国产高清视频在线播放一区| 97超级碰碰碰精品色视频在线观看| 日韩 欧美 亚洲 中文字幕| 黄片播放在线免费| 国产av精品麻豆| 亚洲人成网站在线播放欧美日韩| 精品欧美国产一区二区三| 日本五十路高清| 正在播放国产对白刺激| 国产精品日韩av在线免费观看 | 久久午夜亚洲精品久久| 人人妻人人澡欧美一区二区 | 久久久国产精品麻豆| 激情在线观看视频在线高清| 亚洲一卡2卡3卡4卡5卡精品中文| 搡老岳熟女国产| 91老司机精品| 男女之事视频高清在线观看| 麻豆久久精品国产亚洲av| www国产在线视频色| 天堂动漫精品| 精品国产国语对白av| 久久久久国内视频| 大型黄色视频在线免费观看| 国产成人啪精品午夜网站| 欧美日韩中文字幕国产精品一区二区三区 | 午夜免费成人在线视频| 午夜成年电影在线免费观看| 国产精品爽爽va在线观看网站 | 色哟哟哟哟哟哟| 一进一出抽搐动态| 色在线成人网| 性欧美人与动物交配| 国产又爽黄色视频| 多毛熟女@视频| 国产免费av片在线观看野外av| 亚洲国产精品久久男人天堂| 日韩欧美一区视频在线观看| 久久热在线av| 制服丝袜大香蕉在线| 动漫黄色视频在线观看| 亚洲成国产人片在线观看| 免费高清视频大片| 国产精品乱码一区二三区的特点 | 久久欧美精品欧美久久欧美| 亚洲国产中文字幕在线视频| 露出奶头的视频| 国产精品98久久久久久宅男小说| 很黄的视频免费| 国产免费av片在线观看野外av| 丝袜美足系列| 在线观看免费日韩欧美大片| 亚洲伊人色综图| 黑人欧美特级aaaaaa片| 嫩草影院精品99| 亚洲三区欧美一区| 欧美色视频一区免费| 日本免费一区二区三区高清不卡 | 国产精品影院久久| 99久久国产精品久久久| 天天躁夜夜躁狠狠躁躁| 久久久国产成人精品二区| 国产成人系列免费观看| 精品日产1卡2卡| 午夜a级毛片| av天堂久久9| АⅤ资源中文在线天堂| 热re99久久国产66热| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产精品合色在线| a级毛片在线看网站| 国产亚洲av高清不卡| 一区在线观看完整版| www.自偷自拍.com| 久久久精品国产亚洲av高清涩受| 人人妻,人人澡人人爽秒播| 国产单亲对白刺激| 久久狼人影院| 国产日韩一区二区三区精品不卡| 天堂√8在线中文| 黑人操中国人逼视频| 中文字幕另类日韩欧美亚洲嫩草| 不卡av一区二区三区| 欧美在线一区亚洲| 一区二区三区激情视频| 给我免费播放毛片高清在线观看| 亚洲国产高清在线一区二区三 | 精品久久久久久成人av| 一区福利在线观看| 狂野欧美激情性xxxx| 久久精品国产亚洲av香蕉五月| 女警被强在线播放| 啦啦啦免费观看视频1| 级片在线观看|