龍香月，黃康有,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)系，在未來仍然需要開展更多的研究。