劉穎 田斌 歐光龍

摘 要：為揭示濕潤(rùn)常綠闊葉林和半濕潤(rùn)常綠闊葉林替代分布的氣候制約變量，該研究選擇其代表性優(yōu)勢(shì)樹(shù)種青岡（Cyclobalanopsis glauca）和滇青岡（C. glaucoides）為研究對(duì)象，收集兩個(gè)物種的標(biāo)本分布點(diǎn)數(shù)據(jù)和19個(gè)生物氣候變量圖層數(shù)據(jù)，運(yùn)用MaxEnt模型，模擬其潛在分布區(qū)，通過(guò)判別分析（DFA）、方差分析（ANOVA）和核密度分析三種方法量化兩個(gè)樹(shù)種的生態(tài)位差異，解釋兩類植被地理替代分布規(guī)律及其主導(dǎo)氣候變量。結(jié)果表明：（1）青岡和與滇青岡MaxEnt擬合的AUC值分別為0.995和0.986，準(zhǔn)確預(yù)測(cè)了兩物種的潛在分布;青岡適宜于20°～30° N的亞熱帶及北熱帶地區(qū)，滇青岡集中分布于云貴高原亞熱帶區(qū)域。（2）DFA和ANOVA的結(jié)果一致，表明晝夜溫差月均值（bio2）、等溫性（bio3）、降水季節(jié)性變化（bio15）是造成二者生態(tài)位分異的主導(dǎo)氣候變量。（3）核密度分析顯示晝夜溫差月均值（bio2）、等溫性（bio3）、溫度季節(jié)性變動(dòng)系數(shù)（bio4）是影響兩物種生態(tài)位分化的主導(dǎo)溫度變量;降水季節(jié)性變化（bio15）、年降水量（bio12）、最干月降水量（bio14）、最干季降水量（bio17）、最冷季降水量（bio19）是兩物種生態(tài)位分化的主導(dǎo)降水變量。上述結(jié)果說(shuō)明等溫性、溫度季節(jié)性變化、降水季節(jié)性變化等表征氣候季節(jié)性差異的變量是造成濕潤(rùn)與半濕潤(rùn)常綠闊葉林出現(xiàn)替代分布的主要原因。

關(guān)鍵詞： 青岡， 滇青岡， 濕潤(rùn)常綠闊葉林， 半濕潤(rùn)常綠闊葉林， 生態(tài)位模擬， 替代分布

中圖分類號(hào)：Q948

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-3142（2022）03-0460-10

Displacement distribution and climate explanation on humid? and semi-humid evergreen broadleaved

forests using niche model of Cyclobalanopsis glauca and C. glaucoides in China

LIU Ying， TIAN Bin， OU Guanglong*

（ 1. Key Laboratory of State Forestry Administration on Biodiversity Conservation in Southwest

China， Southwest Forestry University， Kunming 650224， China ）

Abstract：To explain the geographical displacement law and the dominant climate variables between humid evergreen broadleaved forests （HEBF） and semi-humid evergreen broadleaved forests（SHEBF）， Cyclobalanopsis glauca and C. glaucoides were selected as research objects because the both species are the typical dominant tree species of the HEBF and SHEBF， then the specimens data of the both species and the 19 bio-climate variables map layer were collected to simulate their distribution area using MaxEnt model. Moreover， the niche differences were quantitatively analyzed by the kernel density analysis， discriminant function analysis （DFA） and analysis of variance （ANOVA） to reveal the dominant climate variables which affect the niche differentiation between the both species， then the geographical displacement law and the dominant climate variables between the HEBF and SHEBF were analyzed. The results were as follows： （1） The AUC values of C. glauca and C. glaucoides were 0.995 and 0.986 respectively， which accurately predicted the potential distribution of the both species; C. glauca was suitable for the subtropical and north-tropical regions of 20°to 30° N， but C. glaucoides was mainly distributed in subtropical area of Yunnan-Guizhou Plateau. （2）The results of DFA and ANOVA were consistent， and the mean diurnal range （bio2）， isothermality （bio3） and the precipitation seasonality （bio15） were the climate variables that promoted the niche differentiation of the two species. （3） Kernel density analysis showed the mean diurnal range （bio2）， isothermality （bio3） and temperature seasonality （bio4） were the dominant temperature variables affecting the niche differentiation of the both species; and the precipitation seasonality （bio15）， annual precipitation （bio12）， precipitation of the driest month （bio14）， precipitation of the driest quarter （bio17） and precipitation of the coldest quarter （bio19） also significantly affected the niche differentiation of the both species. All the above results indicate that isothermality， the seasonal differences of temperature and precipitation are the main variables determining geographical displacement distribution of the HEBF and SHEBF from the niche simulation of C. glauca and C. glaucoides.

Key words： Cyclobalanopsis glauca， C. glaucoides， humid evergreen broadleaved forest， semi-humid evergreen broadleaved forest， niche model， displacement distribution

植被地理替代分布呈現(xiàn)出水平和垂直變化特征（王鐵娟，2004;李靜等，2006;劉彬彬等，2013;于明茜，2015;李緩，2019），造成地理替代現(xiàn)象的原因比較復(fù)雜，包括地理阻隔、遺傳分化、氣候以及地史等多方面因素（羅艷和周浙昆，2001;李緩，2019），其中氣候是預(yù)測(cè)樹(shù)木及植被類型分布的主要因子，并與其他驅(qū)動(dòng)因素相互作用影響森林的功能和動(dòng)態(tài)過(guò)程，從而成為植被替代分布形成的主導(dǎo)因子之一（Boisvert-Marshet et al.， 2014;Grünig et al.， 2017;鄭維艷和曹坤芳，2020;Zou et al.， 2020），明晰植被分布與氣候的關(guān)系并揭示其影響機(jī)制對(duì)于準(zhǔn)確分析植被分布規(guī)律、植被功能過(guò)程及其生態(tài)恢復(fù)等具有重要理論和實(shí)踐意義（Fang et al.，2002;Buitenwerf & Higgins，2016）。

亞熱帶常綠闊葉林作為森林生態(tài)系統(tǒng)中的重要組成部分，在維持生態(tài)系統(tǒng)平衡與穩(wěn)定方面發(fā)揮著不可替代的作用（丁圣彥和宋永昌，2004;Wang et al.， 2007），中國(guó)作為常綠闊葉林分布最為廣泛的國(guó)家，在秦嶺-淮河線以南，青藏高原以東，云南、福建、東海岸沿線以內(nèi)的廣大區(qū)域成片分布（吳征鎰，1980;Dunmei et al.， 2012;Zhu，2017），其中東部濕潤(rùn)常綠闊葉林和西部半濕潤(rùn)常綠闊葉林是主要的兩個(gè)地帶性植被亞型，且存在明顯的替代分布特征（吳征鎰，1980;吳征鎰和朱彥丞，1987;李昌華，1997;宋永昌等，2005;曾覺(jué)民，2018）。氣候和海拔的梯度變化會(huì)造成植被的替代分布差異（Ge & Xie， 2017;Elsen et al.， 2018），F(xiàn)ang & Yoda（1991）揭示了中國(guó)常綠闊葉林分布上限與氣候因子的關(guān)系，倪健和宋永昌（1997）通過(guò)青岡的地理分布與氣候因子的關(guān)系得出常綠闊葉林對(duì)氣候變化的反應(yīng)不十分敏感和劇烈，且寒冷指數(shù)等反映整體年際氣候變化的指數(shù)不能說(shuō)明其對(duì)分布界限的影響（Chiu et al.， 2014），而反映氣候的季節(jié)變化甚至日變化特征的氣候指標(biāo)則對(duì)植被分布變異具有較大影響（Ernakovich et al.， 2014;Chan et al.， 2016;Allen et al.， 2017;Ge et al.， 2019）。由此可見(jiàn)，學(xué)者們已經(jīng)對(duì)亞熱帶常綠闊葉林的替代分布現(xiàn)象以及常綠闊葉林這一植被型與氣候變量的關(guān)系進(jìn)行了較為詳細(xì)的闡述，但對(duì)于植被亞型間替代分布的主導(dǎo)變量及其替代分布影響機(jī)制，尤其是植被亞型替代分布的季節(jié)性氣候變化影響的研究卻鮮有報(bào)道。

物種分布模型因其可以準(zhǔn)確預(yù)測(cè)和分析物種的適宜生境范圍，模型運(yùn)行穩(wěn)定，還可揭示植被分布的影響氣候變量，所以被廣泛應(yīng)用于生態(tài)學(xué)、保護(hù)生物學(xué)和生物地理學(xué)等領(lǐng)域（Elith & Leathwick， 2009）。目前，最大熵模型（MaxEnt）、基于遺傳算法的規(guī)則組合模型（GARP）、CLIMEX模型、生態(tài)位因子分析模型（ENFA）以及生物氣候分析和預(yù)測(cè)系統(tǒng)模型（BIOCLIM）是最常用的五種物種分布模型。其中，MaxEnt模型可以通過(guò)物種的已知分布數(shù)據(jù)和氣候變量建立物種地理分布與氣候之間的聯(lián)系（Phillips et al.， 2008），且能在大尺度空間上獲取更多的物種信息，為建立系統(tǒng)全面的物種與環(huán)境的關(guān)系提供更為豐富的環(huán)境信息（莊鴻飛等，2018），是目前物種分布模擬中認(rèn)可度最高、擬合度較好的模型（Phillips et al.， 2006;曹向鋒等，2010）。

基于優(yōu)勢(shì)種的生態(tài)位模擬及差異分析，揭示造成其分布差異的氣候制約因素，尤其是氣候季節(jié)性變化的影響機(jī)制對(duì)于闡明中國(guó)東西部?jī)深惖貛猿＞G闊葉林植被亞型的替代分布規(guī)律具有重要理論和實(shí)踐意義。鑒于此，本研究選擇濕潤(rùn)和半濕潤(rùn)常綠闊葉林中具有代表性的優(yōu)勢(shì)樹(shù)種青岡（Cyclobalanopsis glauca）和滇青岡（C. glaucoides）為研究對(duì)象，收集標(biāo)本分布點(diǎn)數(shù)據(jù)，結(jié)合氣候變量，運(yùn)用MaxEnt模型，模擬預(yù)測(cè)兩者的適生區(qū)范圍，分析造成兩個(gè)樹(shù)種生態(tài)位分化的氣候因素，進(jìn)而探討濕潤(rùn)與半濕潤(rùn)常綠闊葉林的替代分布格局及其與氣候條件的相互關(guān)系，揭示其替代分布的主導(dǎo)因素，從而為亞熱帶常綠闊葉林分布和影響機(jī)制研究提供重要的科學(xué)依據(jù)。

1 材料與方法

1.1 標(biāo)本點(diǎn)數(shù)據(jù)

收集中國(guó)數(shù)字植物標(biāo)本數(shù)據(jù)庫(kù)（CVH，http：//www.cvh.org.cn）中青岡與滇青岡的標(biāo)本信息，篩除樣本采集地存疑和出現(xiàn)重復(fù)的樣點(diǎn)，獲取精確樣本記錄地點(diǎn)的經(jīng)緯度坐標(biāo)，作為物種生境模擬的基礎(chǔ)數(shù)據(jù)。共獲取兩個(gè)物種603個(gè)標(biāo)本點(diǎn)，其中包括432個(gè)青岡標(biāo)本點(diǎn)和171個(gè)滇青岡標(biāo)本點(diǎn)（圖1：D）。

1.2 氣候變量數(shù)據(jù)

通過(guò)Worldclim（http：//www.worldclim.org/）下載收集19個(gè)生物氣候圖層（表1），數(shù)據(jù)空間分辨率為30″（Hijmans et al.， 2005）。

1.3 生態(tài)位模型

利用物種分布模型，以物種標(biāo)本點(diǎn)和19個(gè)氣候變量為數(shù)據(jù)，擬合青岡與滇青岡的適宜生境。在MaxEnt 3.3.1中使用默認(rèn)設(shè)置為評(píng)估模型的有效性（Phillips et al.， 2008），設(shè)置75%的物種分布點(diǎn)作為訓(xùn)練數(shù)據(jù)集，剩余25%作為交互檢驗(yàn)的數(shù)據(jù)集，運(yùn)行10次重復(fù)進(jìn)行交互驗(yàn)證（Pearson et al.， 2007）。MaxEnt模型提供刀切法（Jackknife）檢驗(yàn)氣候變量對(duì)模型構(gòu)建的貢獻(xiàn)率（Miguel et al.， 2010），并將模擬結(jié)果進(jìn)行數(shù)據(jù)轉(zhuǎn)換后生成物種適宜生境分布圖，在ArcGIS 10.5中利用分類功能進(jìn)行適生等級(jí)分類，共劃分為5個(gè)等級(jí)：非適生區(qū)0.00～0.02，低適生區(qū)0.02～0.25，中適生區(qū)0.25～0.50，高適生區(qū)0.50～0.75，最適生區(qū)>0.75。

1.4 生態(tài)位差異分析

本研究采用三種統(tǒng)計(jì)方法分析青岡與滇青岡的生態(tài)位差異。首先，利用SPSS 22軟件用判別分析（DFA）評(píng)估各氣候變量對(duì)物種分布的貢獻(xiàn)度;其次，運(yùn)用SPSS 22軟件，用方差分析（ANOVA）定量描述兩物種在各生物氣候因子的種內(nèi)及種間差異;最后，基于Eviews 10軟件，采用核密度分析（kernel density analysis）將樣本數(shù)據(jù)進(jìn)行擬合，得到整個(gè)樣本集的概率密度函數(shù)，通過(guò)核密度曲線描述19個(gè)氣候變量的生態(tài)位變化幅度，由此確定兩物種在各氣候變量中的差異。

2 結(jié)果與分析

2.1 分布區(qū)預(yù)測(cè)

通過(guò)生態(tài)位模擬，兩個(gè)物種平均模型的測(cè)試AUC值均高于0.9，其中，青岡模擬的AUC值為0.995，滇青岡模擬的AUC值為0.986，表明預(yù)測(cè)結(jié)果具有較高的可靠性，可得到兩個(gè)物種準(zhǔn)確的適生區(qū)分布范圍（圖1）。青岡分布于中國(guó)亞熱帶地區(qū)及熱帶北緣，以及日本南部和印度北部，其中最適生區(qū)是中國(guó)貴州東部以東、南嶺以北和秦嶺以南，基本涵蓋了中國(guó)亞熱帶東部地區(qū)（圖1：A）;滇青岡的分布區(qū)則集中在中國(guó)西南部云南、四川、貴州、重慶、西藏等地，以及喜馬拉雅山地區(qū)，最適生區(qū)則以中國(guó)云南中北部、四川南部和貴州西北部為主（圖1：B）。從兩個(gè)樹(shù)種的最適生區(qū)（>0.75分布概率）來(lái)看，二者幾乎在中國(guó)云南東部和貴州西部形成明顯的替代分布（圖1：C）。

2.2 生態(tài)位差異分析

MaxEnt模擬的刀切法（Jackknife）給出19個(gè)氣候變量對(duì)青岡和滇青岡的相對(duì)貢獻(xiàn)率（表2）。

青岡的分布主要受最干月降水量（bio14）、年降水量（bio12）和最暖季降水量（bio18）的影響，三者的貢獻(xiàn)率分別為38.0%、22.1%和15.8%，累計(jì)貢獻(xiàn)率為75.9%，這說(shuō)明了降水及降水的季節(jié)分配是影響青岡分布的主導(dǎo)變量;滇青岡的分布主要受氣溫年變化范圍（bio7）、最暖季降水量（bio18）和溫度季節(jié)性變動(dòng)系數(shù)（bio4）的影響，三者的貢獻(xiàn)率分別為25.3%、25.1%和18.0%，累計(jì)貢獻(xiàn)率為68.4%，這說(shuō)明滇青岡分布受溫度和降水量季節(jié)性變化的限制。

通過(guò)判別函數(shù)分析（DFA）可以看出（表2），造成兩物種間分布差異的主要?dú)夂蜃兞渴堑葴匦裕╞io3）、晝夜溫差月均值（bio2）和降水量季節(jié)性變化（bio15），其差異貢獻(xiàn)率的絕對(duì)值分別為0.81、0.73和0.71;溫度季節(jié)性變動(dòng)系數(shù)（bio4）、最熱月最高氣溫（bio5）、最暖季平均氣溫（bio10）、最干月降水量（bio14）、最干季降水量（bio17）和最冷季降水量（bio19）對(duì)兩物種分布差異貢獻(xiàn)率的絕對(duì)值均在0.6～0.7之間。可見(jiàn)，反映溫度和降水的季節(jié)性差異是造成兩物種分布差異的主要?dú)夂蜃兞俊?/p>

此外，通過(guò)方差分析發(fā)現(xiàn)，晝夜溫差月均值（bio2）、等溫性（bio3）、降水量季節(jié)性變化（bio15）的種間差異大于種內(nèi)差異，最暖季降水量（bio18）種間差異與種內(nèi)差異的方差值相等，而其余變量的種間差異均大于種內(nèi)差異（表3），這說(shuō)明青岡與滇青岡生態(tài)位分化的原因是兩物種對(duì)溫度變化幅度以及降水季節(jié)性分配的適應(yīng)存在差異性。

核密度分析直觀反映兩個(gè)物種對(duì)19個(gè)氣候變量的生態(tài)位幅度變化（圖2）。兩物種在晝夜溫差月均值（bio2）、等溫性（bio3）、溫度季節(jié)性變動(dòng)系數(shù)（bio4）三個(gè)溫度變量表現(xiàn)出明顯的分異，這說(shuō)明溫度年內(nèi)變化及日內(nèi)變化差異是造成二者分布差異的主要溫度制約因素; 而其余溫度變量變化幅度較為相近，尤其在最冷月最低氣溫（bio6）上較為相似，說(shuō)明兩樹(shù)種對(duì)極端低溫的適應(yīng)是相似的（圖2）。另外，兩物種對(duì)降水變量上的適應(yīng)差異較大，尤其是在降水量季節(jié)性變化（bio15）上分異最為突出，且滇青岡在每個(gè)降水變量上表現(xiàn)出較高的波動(dòng)，說(shuō)明滇青岡對(duì)降水季節(jié)性變化產(chǎn)生了適應(yīng)。綜上所述，通過(guò)兩物種在不同氣候變量生態(tài)位幅度變化差異可以看出，表征溫度和降水年內(nèi)或日內(nèi)分配差異的變量是造成二者分布差異的主導(dǎo)氣候變量。

3 討論與結(jié)論

自Phillips等（2006）首次引入MaxEnt模型預(yù)測(cè)物種潛在分布之后，讓生態(tài)位量化研究成為可能（Kozak et al.， 2008），本研究青岡和滇青岡生境模擬的AUC值分別為0.995和0.986，均高于0.90，劉艷梅等（2018）認(rèn)為，AUC高于0.9說(shuō)明兩個(gè)物種的生境預(yù)測(cè)結(jié)果具有較高的可信度。本研究預(yù)測(cè)的青岡與滇青岡適宜生境中，在中國(guó)的分布范圍與《中國(guó)植物志》的分布區(qū)位置基本一致（陳煥鏞和黃成就，1998），國(guó)外分布向南延伸至印度、緬甸北部，向東、東北散布到朝鮮半島、日本，與周浙昆（1992）繪制的青岡屬物種的分布區(qū)相符，進(jìn)一步印證了本文青岡與滇青岡生境模擬的準(zhǔn)確性。此外，從青岡和滇青岡的最適分布區(qū)（>0.75）可以看出，青岡主要分布于20°～30° N之間的亞熱帶及北熱帶地區(qū)，滇青岡集中分布在中國(guó)滇中地區(qū)以及四川南部，二者以中國(guó)云南東部、貴州西部為界形成明顯分異，這也為分析兩個(gè)物種的生態(tài)位差異，以及基于這兩個(gè)代表性優(yōu)勢(shì)樹(shù)種的濕潤(rùn)和半濕潤(rùn)常綠闊葉林的替代分布規(guī)律分析提供了準(zhǔn)確的數(shù)據(jù)保障。

本研究通過(guò)影響青岡與滇青岡分布的19個(gè)氣候因子的方差分析（ANOVA）、判別分析（DFA）和核密度分析發(fā)現(xiàn)，等溫性（bio3）和降水量季節(jié)性變化（bio15）是兩物種之間的差異貢獻(xiàn)率最高的溫度變量和降水變量，其貢獻(xiàn)率的絕對(duì)值分別為0.81和0.71，晝夜溫差月均值（bio2）的貢獻(xiàn)率的絕對(duì)值也達(dá)到了0.73，且三個(gè)氣候變量的種間差異顯著大于種內(nèi)差異;核密度分析表明青岡和滇青岡在晝夜溫差月均值（bio2）、等溫性（bio3）、溫度季節(jié)性變動(dòng)系數(shù)（bio4）、降水量季節(jié)性變化（bio15）等氣候變量上呈現(xiàn)出顯著的生態(tài)位分化。氣候及因?yàn)楹０尾町悗?lái)的氣候因子間接變化會(huì)影響植被分布（Ge & Xie， 2017;Elsen et al.， 2018），F(xiàn)ang & Yoda（1991）研究發(fā)現(xiàn)溫度和降水因素限制了常綠闊葉林的分布界限;倪健和宋永昌（1997）基于常綠闊葉林的主要優(yōu)勢(shì)種與氣候指數(shù)的關(guān)系將青岡和滇青岡分別歸為低中溫濕潤(rùn)型和低溫濕潤(rùn)型，這些研究說(shuō)明了溫度因子是造成二者分異的主導(dǎo)因子，但是該研究?jī)H采用反映溫度和降水總體特征的Kira指數(shù)等氣候指數(shù)來(lái)分析。但反映整體年際氣候變化的指數(shù)不能很好解釋物種及植被分布界限的影響（Chiu et al.， 2014），不能解釋溫度和降水的季節(jié)分配差異對(duì)物種及植被分布的影響，而氣候的季節(jié)性變化甚至是日變化特征都會(huì)對(duì)植被的分布產(chǎn)生影響（Ernakovich et al.， 2014;Chan et al.， 2016;Allen et al.， 2017;Ge et al.， 2019），甚至起到相反的作用（Chan et al.， 2016）?？梢?jiàn)，晝夜溫差月均值（bio2）、等溫性（bio3）、降水量季節(jié)性變化（bio15）等表征溫度和氣候季節(jié)性差異的氣候變量是造成二者生態(tài)位分化的主導(dǎo)變量，進(jìn)而限制了青岡與滇青岡的地理分布范圍，使得兩個(gè)物種的集中分布區(qū)在中國(guó)東西部呈現(xiàn)明顯的替代分布。

優(yōu)勢(shì)種是對(duì)群落生境及其功能發(fā)揮具有重要控制作用的物種（宋永昌，2016），基于優(yōu)勢(shì)樹(shù)種分布可以揭示植被分布規(guī)律及影響因素（倪健和宋永昌，1997）。滇青岡作為西部半濕潤(rùn)常綠闊葉林區(qū)的代表性物種，其集中分布的云貴高原具有海拔高和降水季節(jié)分配不均的特點(diǎn)，在每年11月至次年4月份存在一個(gè)相對(duì)干旱的旱季，植被大多具有耐旱特性（吳征鎰和朱彥丞，1987），而青岡集中分布的亞熱帶東部地區(qū)則是濕潤(rùn)型氣候，海拔較低、熱量充足、生境濕潤(rùn)（倪健和宋永昌，1997）。本研究發(fā)現(xiàn)溫度和氣候季節(jié)性差異的氣候變量是造成青岡和滇青岡兩個(gè)物種生態(tài)位分化的主導(dǎo)變量，也就是說(shuō)中國(guó)亞熱帶區(qū)域自東向西等溫性逐漸降低、溫度季節(jié)性變異增大、降水量相對(duì)減少，尤其是冬春兩季降水顯著減少，使得常綠闊葉林植被亞型呈現(xiàn)由喜溫喜濕逐步變?yōu)橄矞啬秃底兓?，從而使得濕?rùn)常綠闊葉林到半濕潤(rùn)常綠闊葉林沿經(jīng)度和海拔梯度的替代分布變化。由此可見(jiàn)，反映溫度和降水的年際變化差異會(huì)對(duì)常綠闊葉林替代分布產(chǎn)生一定影響，但溫度和降水的季節(jié)分配差異才是青岡與滇青岡生態(tài)位分化的主要?dú)夂蜃兞浚彩窃斐蓾駶?rùn)常綠闊葉林與半濕潤(rùn)常綠闊葉林替代分布的重要?dú)夂蛞蛩亍?/p>

此外，中國(guó)亞熱帶常綠闊葉林區(qū)域跨度大，水熱條件差異大（蔡永立和達(dá)良俊，2002），常綠闊葉林中普遍存在替代分布規(guī)律，以青岡與滇青岡這對(duì)優(yōu)勢(shì)樹(shù)種表征濕潤(rùn)常綠闊葉林和半濕潤(rùn)常綠闊葉林的替代分布具有一定的代表性，有助于學(xué)者在大尺度上清晰認(rèn)識(shí)濕潤(rùn)與半濕潤(rùn)常綠闊葉林植被分布特點(diǎn)與氣候的相互關(guān)系，并明確造成常綠闊葉林替代分布的氣候因素，但不同植被分布往往還受到地形、土壤條件、物種間關(guān)系等其他因素的影響（丁圣彥和宋永昌，2004;Wei et al.， 2020;Long et al.， 2021），考慮多種環(huán)境因子及森林結(jié)構(gòu)差異，從多種角度分析解釋其替代規(guī)律，揭示中國(guó)常綠闊葉林的分布變化規(guī)律有待進(jìn)一步研究。

參考文獻(xiàn)：

ALLEN K， DUPUY JM， GEI M， et al.， 2017. Will seasonally dry tropical forests be sensitive or resistant to future changes in rainfall regimes？ [J]. Environ Res Lett， 12： 023001.

BOISVERT-MARSH L， PRI C，? BLOIS SD， et al.， 2014. Shifting with climate？ Evidence for recent changes in tree species distribution at high latitudes[J]. Ecosphere， 5（7）： art83.

BUITENWERF R， HIGGINS SI， 2016. Convergence among global biogeographical realms in the physiological niche of? evergreen and deciduous vegetation[J]. Global Ecol Biogeogr， 25： 704-715.

CAI YL， DA LJ， 2002. Ecological leaf anatomy of six evergreen species of Fagaceae in the eastern subtropical area of China[J]. Chin J Appl Environ Biol， 8（5）： 460-466. [蔡永立， 達(dá)良俊， 2002. 亞熱帶東部殼斗科6種常綠植物葉的生態(tài)解剖[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào)， 8（5）： 460-466.]

CAO XF， QIAN GL， HU BS， et al.， 2010. Prediction of potential suitable distribution area of Flaveria bidentis in China based on niche models[J]. Chin J Appl Ecol， 21（12）： 3063-3069. [曹向鋒， 錢國(guó)良， 胡白石， 等， 2010. 采用生態(tài)位模型預(yù)測(cè)黃頂菊在中國(guó)的潛在適生區(qū)[J]. 應(yīng)用生態(tài)學(xué)報(bào)， 21（12）： 3063-3069.]

CHAN WP， CJEN IC， COLWELL RK， et al.， 2016. Seasonal and daily climate variation have opposite effects on species elevational range size[J]. Science， 351（6280）： 1437-1439.

CHEN HY， HUANG CJ， 1978. Flora Reipublicae Popularis Sinicae （Vol. 22）[M]. Beijing： Science Press： 315-317. [陳煥鏞， 黃成就， 1978. 中國(guó)植物志（第22卷）[M]. 北京： 科學(xué)出版社： 315-317]

CHIU CA， CHIOU CR， LIN JR， et al.，? 2014. Coldness index does not indicate the upper limit of evergreen broad-leaved forest on a subtropical island[J]. J For Res， 19： 115-124.

DING SY， SONG YC， 2004. Research advances in vegetation dynamic of evergreen broad-leaved forest[J]. Acta Ecol Sin， （8）： 1769-1779. [丁圣彥， 宋永昌， 2004. 常綠闊葉林植被動(dòng)態(tài)研究進(jìn)展[J]. 生態(tài)學(xué)報(bào)， （8）： 1769-1779.]

DUNMEI L， JIANGSHAN L， HELENCE C， et al.， 2012. Topographic variation in aboveground biomass in a subtropical evergreen broad-Leaved forest in China[J]. PLoS ONE， 7（10）： e48244.

ELITH J， LEATHWICK JR， 2009. Species distribution models： ecological explanation and prediction across space and time[J]. Ann Rev Ecol Evol Syst， 40： 677-697.

ELSEN PR， MONAHAN WB， MERENLENDER AM，2018. Global patterns of protection of elevational gradients in mountain ranges[J]. Proc Natl Acad Sci USA， 115： 6004-6009.

ERNAKOVICH JG， HOPPING KA， BERDANIER AB， et al.， 2014. Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change[J]. Global Chang Biol， 20（10）： 3256- 3269.

FANG JY， SONG YC， LIU HY， et al.， 2002. Vegetation-climate relationship and its application in the division of vegetation zone in China[J]. J Integr Plant Biol， 44（9）： 1105-1122.

FANG JY， YODA K， 1991. Climate and vegetation in China V. Effect of climatic factors on the upper limit of distribution of evergreen broadleaf forest[J]. Ecol Res， 6： 113-125.

GE JL， BERG B， XIE ZQ， 2019. Climatic seasonality is linked to the occurrence of the mixed evergreen and deciduous broad-leaved forests in China[J]. Ecosphere， 10（9）： e02862.

GE JL， XIE ZQ， 2017. Geographical and climatic gradients of evergreen versus deciduous broadleaved tree species in subtropical China： implications for the definition of the mixed forest[J]. Ecol? Evol， 7： 3636-3644.

GRNIG M， BEERLI N，BALLESTROS-MEJIA L， et al.， 2017. How climatic variability is linked to the spatial distribution of range sizes： seasonality versus climate change velocity in sphingid moths[J]. J Biogeogr， 44： 2441-2450.

HIJMANS RJ， CAMERON SE， PARRA JL， et al.， 2005. Very high resolution interpolated climate surfaces for global land areas[J]. Int J Climatol， 25： 1965-1978.

KOZAK KH， GRAHAM CH， WIENS JJ， 2008. Integrating GIS-based environmental data into evolutionary biology[J]. Trends Ecol Evol， 23： 141-148.

LI CH， 1997. The distribution of evergreen broad-leaves forests in East Asia[J]. Nat Resouc， （2）： 36-45. [李昌華， 1997. 亞洲東部常綠闊葉林的分布[J]. 自然資源， （2）： 36-45.]

LI J， CAO HL， LIAN JY， et al.， 2006. Study on classification and geographical distribution of Castanopsis hystrix forest in China[J]. Guihaia， 26（1）： 22-27. [李靜， 曹洪麟， 練琚蕍， 等， 2006. 中國(guó)刺栲林的分類與分布初探[J]. 廣西植物， 26（1）： 22-27.]

LI Y， 2019. Species composition and geographical differentiation of the Fagaceae community in subtropical evergreen broad-leaved forest in China[D]. Hefei： Anhui Agricultural Univeristy. [李緩， 2019. 中國(guó)亞熱帶常綠闊葉林殼斗科群落的物種組成及其地理分異[D]. 合肥： 安徽農(nóng)業(yè)大學(xué).]

LIU BB， LOU LH， LIU GY， 2013. A floristic analysis of Fagaceae in Zhejiang Province， China[J]. J Zhejiang A & F Univ， 30（5）： 698-705. [劉彬彬， 樓爐煥， 劉廣寧， 2013. 浙江省殼斗科植物區(qū)系特征分析[J]. 浙江農(nóng)林大學(xué)學(xué)報(bào)， 30（5）： 698-705.]

LIU YM， ZHOU SD， XIE DF， et al.， 2018. Potential distribution of Fritillaria unibracteata predicted by the MaxEnt model[J]. Guihaia， 38（3）： 352-360. [劉艷梅， 周頌東， 謝登峰， 等， 2018. 基于最大熵模型（MaxEnt）預(yù)測(cè)暗紫貝母的潛在分布[J]. 廣西植物， 38（3）： 352-360.]

LONG JS， TANG MP， CHEN GS， 2021. Influence of strata-specific forest structural features on the regeneration of the evergreen broad-leaved forest in Tianmu Mountain[J]. PLoS ONE， 16（2）： e0247339.

LUO Y， ZHOU ZK， 2001.Phytogeograph of Quercus subg. Cyclobalanopsis[J]. Plant Divers， 23（1）： 1-16. [羅艷， 周浙昆， 2001. 青岡亞屬植物的地理分布[J]. 植物分類與資源學(xué)報(bào)， 23（1）： 1-16.]

MIGUEL B， ARAUJO RG， PEARSON WT， et al.， 2010. Validation of species-climate impact models under climate change[J]. Global Change Biol， 11（9）： 1504-1513.

NI J， SONG YC， 1997. Relationships between geographical distribution of Cyclobalanopsis glauca and climate in China[J]. Acta Bot Sin， 39（5）： 451-460. [倪健， 宋永昌， 1997. 中國(guó)青岡的地理分布與氣候的關(guān)系[J]. 植物學(xué)報(bào)， 39（5）： 451-460.]

PEARSON RG， RAXWORTHY CJ， NAKAMURA M， et al.， 2007. Predicting species distributions from small numbers of occurrence records： a test case using cryptic geckos in Madagascar[J]. J Biogeogr， 34（1）： 102-117.

PHILLIPS SJ， 2008. Transfer ability， sample selection bias and background data in presence-only modelling： A response to Peterson [J]. Ecography， 31： 272-278.

PHILLIPS SJ， ANDERSON RP， SCHAPIRE RE， 2006. Maximum entropy modeling of species geographic distributions[J]. Ecol Model， 190： 231-259.

SONG YC， 2016. Vegetation ecology [M]. 2nd ed. Beijing： Higher Education Press： 353-360. [宋永昌， 2016. 植被生態(tài)學(xué) [M]. 2版. 北京： 高等教育出版社： 353-360.]

SONG YC， CHEN XY， WANG XH， 2005. Study on evergreen broad-leaved forests of China： A retrospect and prospect[J]. J East Chin Norm Univ（Nat Sci Ed）， （1）： 1-8. [宋永昌， 陳小勇， 王希華， 2005. 中國(guó)常綠闊葉林研究的回顧與展望[J]. 華東師范大學(xué)學(xué)報(bào)（自然科學(xué)版）， （1）： 1-8.]

WANG TJ， 2004. Studies on kinships， geographical replacement law and reasons of the six sand plants of Artemisia in northern China[D]. Huhhot： Inner Mongolia University. [王鐵娟， 2004. 中國(guó)北方六種沙蒿親緣關(guān)系、替代分布規(guī)律與成因的研究[D]. 呼和浩特： 內(nèi)蒙古大學(xué).]

WANG XH， KENT M， FANG XF， 2007. Evergreen broad-leaved forest in Eastern China： Its ecology and conservation and the importance of resprouting in forest restoration[J]. For Ecol Manag， 245（1-3）： 76-87.

WEI SZ， TIE LH， LIAO J， et al.， 2020. Nitrogen and phosphorus co-addition stimulates soil respiration in a subtropical evergreen broad-leaved forest[J]. Plant Soil， 450（10）： 171-182.

WU ZY， 1980. China vegetation[M]. Beijing： Science Press： 306-355. [吳征鎰， 1980. 中國(guó)植被[M]. 北京： 科學(xué)出版社： 306-355.]

WU ZY， ZHU YC， 1987. Yunnan vegetation[M]. Beijing： Science Press： 193-230. [吳征鎰， 朱彥丞， 1987. 云南植被[M]. 北京： 科學(xué)出版社： 193-230.]

YU MQ， 2015. The xylem anatomy and germination of geographical substitute Garagana species[D]. Lanzhou： Lanzhou University. [于明茜， 2015. 錦雞兒屬植物典型地理替代分布種解剖結(jié)構(gòu)特征和種子萌發(fā)特性研究[D]. 蘭州： 蘭州大學(xué).]

ZENG JM， 2018. The classification system of natural forests and its geographic distribution in Yunnan[J].J SW For Univ（Nat Sci Ed）， 38（6）： 1-18. [曾覺(jué)民， 2018. 云南自然森林分類系統(tǒng)及地理分布研究[J]. 西南林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)）， 38（6）： 1-18.]

ZHENG WY， CAO KF， 2020. Impact of future climate change on potential geographical distribution of four Litsea species in China [J]. Guihaia， 40（11）： 1584-1594. [鄭維艷， 曹坤芳， 2020. 未來(lái)氣候變化對(duì)四種木姜子地理分布的影響[J]. 廣西植物， 40（11）： 1584-1594.]

ZHOU ZK， 1992. Origin phylogeny and dispersal of Quercus from China[J]. Acta Bot Yunnan， 14（3）： 227-236. [周浙昆， 1992. 中國(guó)櫟屬的起源演化及其擴(kuò)散[J]. 云南植物研究， 14（3）： 227-236.]

ZHU H， 2019. Studies on the evergreen broad-leaved forests of Yunnan， Southwestern China[J]. Bot Rev， 85（1）： 1-20.

ZHUANG HF， ZHANG YB， WANG W， et al.， 2018. Optimized hot spot analysis for probability of species distribution under different spatial scales based on MaxEnt model：Manglietia insignis case[J]. Biodivers Sci， 26（9）： 931-940. [莊鴻飛， 張殷波， 王偉， 等， 2018. 基于最大熵模型的不同尺度物種分布概率優(yōu)化熱點(diǎn)分析： 以紅色木蓮為例[J]. 生物多樣性， 26（9）： 931-940.]

ZOU S， ZHANG QM， ZHOU GY， et al.， 2020. Functional composition changes of a subtropical monsoon evergreen broad-leaved forest under environmental change[J]. Forests， 11（2）： 191.

（責(zé)任編輯 周翠鳴）

收稿日期：2021-03-21

基金項(xiàng)目：國(guó)家自然科學(xué)基金（31760206， 31660202）;云南省唐守正院士工作站（ 2018IC066）;云南省萬(wàn)人計(jì)劃青年拔尖人才專項(xiàng)（YNWR-QNBJ-2018-184）[Supported by National Natural Science Foundation of China（31760206， 31660202）; Tang Shouzheng Academician Workstation of Yunnan Province （2018IC066）; Ten-Thousand Talents Program of Yunnan Province （YNWR-QNBJ-2018-184）]。

第一作者： 劉穎 （1997-），碩士研究生，研究方向?yàn)樯纸?jīng)理學(xué)，（E-mail） 1017923405@qq.com。

通信作者：歐光龍，博士，教授，博士研究生導(dǎo)師，研究方向?yàn)樯稚鷳B(tài)學(xué)，（E-mail） olg2007621@126.com。