山西呂梁山草本群落物種多樣性的海拔梯度格局及與環(huán)境因子的關(guān)系

張世雄 秦瑞敏 楊曉艷 溫靜 徐滿厚

摘?要：山地草本群落物種多樣性的海拔梯度格局及其與環(huán)境因子的關(guān)系一直是生態(tài)學(xué)和地理學(xué)研究的熱點(diǎn)問題。為探究呂梁山草本群落不同尺度物種多樣性（α、β、γ多樣性）的海拔梯度格局及其與環(huán)境因子的關(guān)系，該研究采用樣方法在呂梁山從北向南選取管涔山、關(guān)帝山、五鹿山，并對(duì)每個(gè)山地劃分出高、中、低3個(gè)海拔梯度，對(duì)草本群落物種多樣性和環(huán)境因子進(jìn)行調(diào)查分析。結(jié)果表明：（1）α多樣性隨海拔升高呈先升后減的單峰變化格局，峰值出現(xiàn)在中海拔;β多樣性中Cody指數(shù)在1 900～2 000 m間出現(xiàn)峰值，而Bray-Curtis 指數(shù)出現(xiàn)最小值，表明1 900～2 000 m海拔帶是草本群落物種更新速率和組成變化較快的過渡地帶;γ多樣性隨海拔升高也呈先升后減的單峰變化格局（R2=0.406，P<0.01）。（2）海拔對(duì)溫度、濕度變化的影響極其顯著（P<0.01）;海拔、溫度、濕度等環(huán)境因子共同影響物種的α多樣性，其中海拔是主要影響因素（P<0.05），溫度、濕度影響不顯著（P>0.05）。綜上所述，呂梁山草本群落多樣性在海拔梯度上多表現(xiàn)為 “中間膨脹”變化格局，海拔變化是影響多樣性變化的主要因素。

關(guān)鍵詞：呂梁山， 草本群落， 物種多樣性， 海拔梯度， 環(huán)境因子

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

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

文章編號(hào)：1000-3142（2020）12-1860-09

Abstract：Patterns of species diversity in mountain herbaceous communities at different altitude gradients and their relationships with environmental factors was an enduring hot topic in ecological and geographical researches.This topic was explored in the context of the mountains of Guancen， Guandi， and Wulu that extend from north to south across the Lüliang Mountain range in Shanxi Province in China.Each mountain was divided into three altitude gradients：high， middle and low.During the vegetation growing season in 2018， an investigation was conducted on the growth characteristics indicators of plants and hydrothermal factors located at these altitude gradients and species diversity index values （α-， β-， and γ-diversity） were calculated to explore the patterns of species diversity at different altitude gradients?and their relationships with environmental factors.The results were as follows：（1） Unimodal patterns of variation for α-diversity （calculated using the Simpson， Shannon-Wiener and Pielou indexes） that corresponded to increasing altitudes， with the highest species diversity occurring at the mid-altitude gradient; The results calculated for β-diversity using different indexes revealed the highest value for the Cody index and the lowest value for the Bray-Curtis index at altitudes between 1 900 and 2 000 m， indicating that areas located between 1 900 and 2 000 m form a transition zone in which the herbaceous community undergoes a rapid process of species renewal and changes in its composition; The results for γ-diversity indicated a pattern of unimodal variation in relation to altitude.（2） Changes in altitude gradient had highly significant impacts on changes in temperature and humidity （P<0.01）; Various environmental factors， notably altitude， humidity， and temperature had a combined effects on changes in the values of the α-diversity indexes， with altitude being the main influencing factor （P<0.05）， and changes in temperature and humidity had no significant effects on species diversity （P>0.05）.Therefore， species diversity within the herbaceous community in the Lüliang Mountains predominantly showed a “mid-altitude bulge” pattern， with altitude being the main factor affecting changes in species diversity.The finding of this study can contribute to refining and expanding studies of herbaceous communities in this region and provide inputs for advancing the systematic conservation and management of herbaceous resources.

1.5 數(shù)據(jù)統(tǒng)計(jì)分析

1.5.1 α物種多樣性?α多樣性為生境內(nèi)多樣性，表征局部均質(zhì)生境下的物種數(shù)量。采用相對(duì)高度、相對(duì)多度、相對(duì)蓋度和相對(duì)頻度4個(gè)指標(biāo)計(jì)算物種重要值，進(jìn)而得到α多樣性指數(shù)，主要包括 Simpson指數(shù)、 Shannon-Wiener指數(shù)和 Pielou均勻度指數(shù)，其計(jì)算公式如下。

式中：IV為重要值;rh為相對(duì)高度;ra為相對(duì)多度;rc為相對(duì)蓋度;rf為相對(duì)頻度;Pi為相對(duì)重要值;H′ 為Simpson指數(shù);H為Shannon-Wiener指數(shù);E為Pielou均勻度指數(shù);i 為樣地內(nèi)的植物物種i;S 為樣地內(nèi)所有物種總和。

1.5.2 β多樣性?β多樣性為生境間多樣性，表示環(huán)境梯度下不同生境間群落物種組成的相異性或物種沿環(huán)境梯度的更替速率，通過群落β多樣性的測(cè)度可以比較不同地段的生境多樣性。將不同山地的地理梯度帶劃分：海拔以100 m為間隔，從低到高將實(shí)驗(yàn)樣地劃分成9個(gè)海拔梯度帶（表1）。選擇基于二元屬性數(shù)據(jù)測(cè)度的Cody指數(shù)和基于數(shù)量數(shù)據(jù)的Bray-Curtis指數(shù)進(jìn)行β多樣性的海拔梯度特征分析，計(jì)算公式如下。

式中：βC為Cody指數(shù)，反映物種沿環(huán)境梯度的替代速率;βB-C為Bray-Curtis 指數(shù)，體現(xiàn)群落物種組成的變化;a和b分別為樣地A和B群落各自的物種數(shù)，c為兩個(gè)樣地群落共有的物種數(shù);cIV為樣地A和B共有種中重要值較小者之和，即cIV =∑min（cIVa， cIVb）。

1.5.3 γ多樣性?γ多樣性表示一個(gè)地理空間內(nèi)各生境中的物種豐富度，多指區(qū)域或更大空間尺度下的物種數(shù)量。以呂梁山的9個(gè)海拔梯度帶為單位，作為呂梁山草本群落的區(qū)域尺度，以每個(gè)區(qū)域內(nèi)的總物種數(shù)目（總物種豐富度S）為分析指標(biāo)，分析呂梁山草本群落 γ 多樣性的空間分布格局。

1.5.4 數(shù)據(jù)分析?運(yùn)用IBM SPSS Statistics 24.0和Origin 2015軟件，對(duì)不同數(shù)據(jù)進(jìn)行趨勢(shì)分析、相關(guān)分析、多因素方差分析以及回歸分析等處理和圖形繪制。

2?結(jié)果與分析

2.1 物種多樣性的海拔梯度變化

2.1.1 α多樣性?在垂直空間上，Simpson指數(shù)（H′）、Shannon-Wiener指數(shù)（H）和Pielou均勻度指數(shù)（E）的整體變化趨勢(shì)一致，隨海拔呈二次多項(xiàng)式函數(shù)變化，即中間高兩邊低的單峰變化格局（圖2），其中H′的擬合度最高（R2 =0.611，P<0.01），其次是H（R2 =0.575，P<0.01）、E（R2 =0.122，P>0.05）。E在中海拔較為接近，表明中海拔內(nèi)部環(huán)境較為均勻且具有較高的物種多樣性，與H′、H的變化表征相一致。從多因素方差分析看（表2），在不同山地的不同海拔影響下，H均達(dá)到極顯著水平（P<0.01），H′達(dá)到顯著水平（P<0.05），E均表現(xiàn)為不顯著差異（P>0.05）?？傮w來看，α多樣性指數(shù)隨海拔梯度變化一致，都隨海拔的升高先升高后降低，即在中海拔最高，且不同海拔和山地對(duì)H′、H影響顯著，對(duì)E影響不顯著。

2.1.2 β多樣性?隨著海拔的升高，呂梁山草本群落間Cody指數(shù)（βC）呈先急劇升高后平穩(wěn)變化的趨勢(shì)，最小值在1 400～1 500 m之間，在1 900～2 000 m和2 200～2 300 m海拔帶間達(dá)到最大值，且中海拔帶呈平穩(wěn)變化趨勢(shì);Bray-Curtis 指數(shù)（βB-C）隨著海拔的升高呈先急劇下降后平緩變化的趨勢(shì)，與βC指數(shù)變化相反（圖3），最大值在1 400～1 500 m之間，最小值出現(xiàn)在2 200～2 300 m。整體上，β多樣性在1 400～2 000 m海拔間出現(xiàn)急劇變化，因此，1 900～2 000 m海拔帶是呂梁山草本群落物種組成變化的過渡地帶;在2 000～2 300 m海拔之間差異較小，即中海拔表現(xiàn)為相對(duì)平緩的變化，表明草本群落從低海拔到高海拔的物種更新速率加快，中海拔間群落物種變化相對(duì)較小。

2.1.3 γ多樣性?草本群落的γ多樣性隨著海拔的升高呈二次多項(xiàng)式函數(shù)變化（R2=0.406，P<0.01），即先升高后降低的單峰分布格局（圖4）?？偽锓N豐富度在海拔梯度上的變化范圍為1～14，平均值為9，峰值在2 000～2 100 m的中海拔帶，表明在呂梁山草本群落的物種在垂直空間上分布不均，中海拔帶明顯高于低、高海拔帶。

2.2 物種多樣性與環(huán)境因子的關(guān)系

呂梁山草本群落的α多樣性指數(shù)與空氣溫度均呈負(fù)相關(guān)關(guān)系，即隨著空氣溫度的升高物種多樣性呈減少趨勢(shì);H與土壤溫度呈正相關(guān)關(guān)系，H′、E與土壤溫度呈負(fù)相關(guān)關(guān)系;H′、H、E與空氣和土壤濕度均呈正相關(guān)關(guān)系，即隨著濕度的升高物種多樣性呈增加趨勢(shì)，但水熱因子與α多樣性的關(guān)系均未達(dá)到顯著水平（P>0.05）（表3）。對(duì)α多樣性指數(shù)與水熱因子進(jìn)行回歸分析（圖5），H′和H與空氣溫度、H與空氣濕度均符合二次函數(shù)關(guān)系（P<0.05）。隨著空氣溫度升高，H′和H呈先下降后上升的變化; 隨著空氣濕度的升高，H呈先下降后上升的變化;除圖5中的三種回歸分析外，α多樣性指數(shù)與溫濕度的回歸分析均未達(dá)到顯著水平（P>0.05）。整體而言，溫度、濕度的變化對(duì)α多樣性的變化具有一定影響，但并不顯著。

3?討論

3.1 物種多樣性的空間分布

物種多樣性沿地理梯度的變化包括熱量、水分、地形等多種環(huán)境因子的綜合影響，從而導(dǎo)致物種多樣性的空間分布格局有很大差異性（陳斌等，2016）。沿海拔梯度的分布格局一般有5種形式（唐志堯和方精云，2004），分別是沿海拔梯度先降后升、單峰曲線、單調(diào)升高、單調(diào)下降和無明顯格局（李晉鵬等，2008）。本研究中α多樣性和γ多樣性隨著海拔的升高呈中間高兩邊低的單峰曲線變化趨勢(shì)，在中海拔（2 000～2 300 m）地區(qū)達(dá)到最大值，這與大多數(shù)的研究相一致。在大多數(shù)情況下，低海拔區(qū)域溫度較高而水分相對(duì)缺少，加之受人為干擾的影響大而產(chǎn)生負(fù)面影響（馮建孟等，2006）;而在高海拔地區(qū)主要受熱量限制，太陽輻射強(qiáng)、晝夜溫差大等環(huán)境條件嚴(yán)重影響植物生長發(fā)育，因此物種多樣性也較低。相對(duì)于低海拔和高海拔地區(qū)，中海拔地區(qū)的溫度和降水量適宜，土壤深厚發(fā)育好，保水性較好，其水熱條件更適宜植被生長，加之人為干擾相對(duì)較少，從而形成物種多樣性的峰值區(qū)域（曲波等，2012）。

β多樣性能夠反映隨著環(huán)境梯度變化植物群落間的差異性，或是快速變化的環(huán)境梯度中的物種組成的更替（張仕豪等，2019）。本研究中Cody指數(shù)和Bray-Curtis 指數(shù)在1 400～2 000 m 海拔間出現(xiàn)急劇變化，且在1 900～2 000 m海拔帶間出現(xiàn)極值點(diǎn)，這表明從低海拔（1 400 m）到中海拔（2 000 m）草本群落間物種更新速率加快，物種間組成差異大，且1 900～2 000 m海拔帶是明顯的物種組成變化的過渡地帶。有研究表明山地植被群落間的過渡帶具有高的物種多樣性和更替速率（盧訓(xùn)令等，2010;鄧清月等，2019），這一結(jié)果也與草本群落的α多樣性和γ多樣性在中海拔帶出現(xiàn)峰值的結(jié)果相吻合，這一結(jié)論有利地驗(yàn)證了“中間膨脹效應(yīng)”（Willis & Whittaker，2002）。

3.2 物種多樣性與環(huán)境因子的關(guān)系

海拔的變化影響著山地生態(tài)系統(tǒng)的地理結(jié)構(gòu)和水熱過程，從而影響物種多樣性的空間梯度變化（胡光成等，2009）。劉洋等（2007）在物種豐富度的研究中發(fā)現(xiàn)物種多樣性的變化受到熱量和水分的綜合調(diào)控，不同的因子對(duì)物種豐富度海拔梯度格局的解釋程度不同。OBrien & Whittaker（2000）在南非關(guān)于物種豐富度格局研究中建立的“水分-能量動(dòng)態(tài)”模型認(rèn)為水分和熱量間的關(guān)系控制了植物的生理活動(dòng)，進(jìn)而影響物種多樣性的變化格局。

在山地系統(tǒng)中有大量的環(huán)境因子用于解釋物種多樣性的分布格局，許多表征熱量、水分的氣候因子被認(rèn)為與不同地區(qū)、不同尺度的多樣性格局密切相關(guān)（Currie et al.，2004）。本研究中將海拔作為環(huán)境因子分析了其對(duì)α多樣性的影響，結(jié)果表明海拔對(duì)多樣性變化有顯著的正相關(guān)影響。生境異質(zhì)性高的地區(qū)能夠提供更多的生態(tài)位，氣候隨著海拔的升高而變化，垂直立體氣候特征越明顯，生境類型越多樣化，從而滿足不同物種對(duì)生境的不同需求，使得物種多樣性呈現(xiàn)增加趨勢(shì)（葉鵬程等，2019）。

溫度與α多樣性呈現(xiàn)不顯著的負(fù)相關(guān)，即隨著溫度升高物種多樣性呈降低趨勢(shì);濕度對(duì)物種多樣性的變化產(chǎn)生了正向的影響但不顯著，這與前人的多數(shù)研究結(jié)果相反（Zhao & Fang，2006），這可能是由于研究區(qū)域尺度不同所導(dǎo)致。本研究區(qū)域呂梁山的整體跨度較小，跨緯度僅有

2°～3°，其所反映的水熱因子變化不顯著，從而對(duì)物種多樣性的影響不顯著。因此，為了探明呂梁山環(huán)境因子與植被物種多樣性的關(guān)系，在后續(xù)的研究中應(yīng)把握好研究尺度和因子的選擇問題。

4?結(jié)論

呂梁山草本群落的α多樣性在海拔梯度上呈現(xiàn)單峰變化模式，峰值在中海拔地區(qū);β多樣性指數(shù)的變化表明1 900～2 000 m海拔帶是草本群落物種組成變化的過渡地帶;γ多樣性在海拔梯度上表現(xiàn)為單峰變化模式，草本群落物種多樣性的變化符合“中間膨脹效應(yīng)”。海拔是影響呂梁山草本植物群落物種多樣性的主要因素，溫度、濕度的變化對(duì)物種多樣性的影響不顯著。

參考文獻(xiàn)：

CHEN B， LI HD，CAO XZ， et al.， 2016.Dynamic changes in vegetation coverage in the Yarlung Zangbo River basin based on SPOT-VGT NDVI[J].J Mount Sci， 34 （2）：249-256.[陳斌， 李海東， 曹學(xué)章， 等， 2016.基于SPOT-VGT NDVI的雅魯藏布江流域植被動(dòng)態(tài)變化[J].山地學(xué)報(bào)， 34（2）：249-256.]

CURRIE DJ， MITTELBACH GG， CORNELL HV， et al.， 2004.Predictions and tests of climate-based hypotheses of broad-scale variation in taxonomic richness[J].Ecol Lett， 7（12）： 1121-1134.

DENG QY， ZHANG XL， NIU JJ， et al.， 2019.Species diversity of plant communities along an altitude gradient in Yinmachi Mountain， Northwestern Shanxi， China[J].Ecol Environ， 28（5）：865-872.[鄧清月， 張曉龍， ?？〗?， 等， 2019.晉西北飲馬池山植物群落物種多樣性沿海拔梯度的變化[J].生態(tài)環(huán)境學(xué)報(bào)， 28（5）：865-872.]

FENG JP， WANG XP， XU CD， et al.， 2006.Altitudinal patterns of plant species diversity and community structure on Yulong Mountains， Yunnan， China[J].J Mount Sci， 24（1）：110-116.[馮建孟， 王襄平， 徐成東， 等， 2006.玉龍雪山植物物種多樣性和群落結(jié)構(gòu)沿海拔梯度的分布格局[J].山地學(xué)報(bào)， 24（1）：110-116.]

FRASER LH， PITHER J， JENTSCH A， et al.， 2015.Worldwide evidence of a unimodal relationship between productivity and plant species richness[J].Science， 349（6245）：302-305.

GASTON KJ， DAVIES RG， ORME CDL， et al.， 2007.Spatial turnover in the global avifauna[J].Proc Roy Soc B-Biol Sci， 274（1618）：1567-1574.

HE JS， CHEN WL， 1997.A review of gradient changes in species diversity of land plant communities[J].Acta Ecol Sin， 17（1）：91-99.[賀金生， 陳偉烈， 1997.陸地植物群落物種多樣性的梯度變化特征[J].生態(tài)學(xué)報(bào)， 17（1）：91-99.]

HE YH， YAN M， ZHANG QD， et al.， 2013.Altitudinal pattern of plant species diversity in the Wulu Mountain Nature Reserve， Shanxi， China[J].Acta Ecol Sin， 33（8）：2452-2462.[何艷華， 閆明， 張欽弟， 等， 2013.五鹿山國家級(jí)自然保護(hù)區(qū)物種多樣性海拔格局[J].生態(tài)學(xué)報(bào)， 33（8）：2452-2462.]

HU GC， JIN XM， WAN L， et al.， 2009.The relationship between vegetation growth and the combination of water and heat in the Qilian Mountain[J].J Arid Land Resour Environ， 23（2）：17-20.[胡光成， 金曉媚， 萬力， 等， 2009.祁連山區(qū)植被生長與水熱組合關(guān)系研究[J].干旱區(qū)資源與環(huán)境， 23（2）：17-20.]

JIN SH， LIU T， PANG XP， et al.， 2017.Effects of plateau pika （Ochotona crzoniae） disturbances on plant species diversity and aboveground plant biomass in a Kobresia pygmaea meadow in the Qinghai Lake Region[J].Acta Pratac Sin， 26（5）：29-39.[金少紅， 劉彤， 龐曉攀， 等， 2017.高原鼠兔干擾對(duì)青海湖流域高山嵩草草甸植物多樣性及地上生物量的影響[J].草業(yè)學(xué)報(bào)， 26（5）：29-39.]

KESSLER M， 2000.Elevational gradients in species richness and endemism of selected plant groups in the central Bolivian Andes[J].Plant Ecol， 149（2）：181-193.

LI JP， SHANGGUAN TL， GUO DG， et al.， 2008.Study on relationship between species diversity of plant community and environment in south Lüliang Mountains， Shanxi[J].J Mount Sci， 26（5）：612-619.[李晉鵬， 上官鐵梁， 郭東罡， 等， 2008.山西呂梁山南段植物群落物種多樣性與環(huán)境的關(guān)系[J].山地學(xué)報(bào)， 26（5）：612-619.]

LIU Y， ZHANG J， YANG WQ， 2009.Responses of alpine biodiversity to climate change[J].Biodivers Sci， 17（1）：88-96.[劉洋， 張健， 楊萬勤， 2009.高山生物多樣性對(duì)氣候變化響應(yīng)的研究進(jìn)展[J].生物多樣性， 17（1）：88-96.]

LIU Y， ZHANG YP， HE DM， et al.， 2007.Elevational patterns of species richness and climate[J].Chin Sci Bull， 52（S2）：43-50.[劉洋， 張一平， 何大明， 等， 2007.縱向嶺谷區(qū)山地植物物種豐富度垂直分布格局及氣候解釋[J].科學(xué)通報(bào)， 52（S2）：43-50.]

LU XL， HU N， DING SY， et al.， 2010.The pattern of plant species diversity of Funiu Mountain Nature Reserve[J].Acta Ecol Sin， 30（21）：5790-5798.[盧訓(xùn)令， 胡楠， 丁圣彥， 等， 2010.伏牛山自然保護(hù)區(qū)物種多樣性分布格局[J].生態(tài)學(xué)報(bào)， 30（21）：5790-5798.]

LUO LM， MIAO YJ， WU JS， et al.， 2014.Variation in the biodiversity of montane shrub grassland communities along an altitudinal gradient in a Lhasa River basin valley[J].Acta Prat Sin， 23（6）：320-326.[羅黎鳴， 苗彥軍， 武建雙， 等， 2014.拉薩河谷山地灌叢草地物種多樣性隨海拔升高的變化特征[J].草業(yè)學(xué)報(bào)， 23（6）：320-326.]

OBRIEN EM， WHITTAKER FRJ， 2000.Climatic gradients in woody plant （tree and shrub） diversity：water-energy dynamics， residual variation， and topography[J].Oikos， 89（3）： 588-600.

QU B， MIAO YM， ZHANG QD， et al.， 2012.Plant diversity and its elevational gradient patterns in Wulu Mountain， Shanxi， China[J].Plant Divers Resour， 34（4）：376-382.[曲波， 苗艷明， 張欽弟， 等， 2012.山西五鹿山植物物種多樣性及其海拔梯度格局[J].植物分類與資源學(xué)報(bào)， 34（4）：376-382.]

REY BENAYAS JM， SCHEINER SM， 2002.Plant diversity， biogeography and environment in Iberia：patterns and possible causal factors[J].J Veg Sci， 13（2）：245-258.

SHELDON KS， YANG S， TEWKSBURY JJ， 2011.Climate change and community disassembly：impacts of warming on tropical and temperate montane community structure[J].Ecol Lett， 14（12）：1191-1200.

SU C， ZHANG XY， MA WH， et al.， 2018.Altitudinal pattern and environmental interpretation of species diversity of scrub communty in the Helan Mountains， China[J].J Mount Sci， 36（5）：699-708.[蘇闖， 張芯毓， 馬文紅， 等， 2018.賀蘭山灌叢群落物種多樣性海拔格局及環(huán)境解釋[J].山地學(xué)報(bào)， 36（5）：699-708.]

TANG ZY， FANG JY， CHI XL， et al.， 2012.Patterns of plant beta diversity along elevational and latitudinal gradients in mountain forests of China[J].Ecography， 35（12）：1083-1091.

TANG ZY， FANG JY， 2004.A review on the elevational patterns of plant species diversity[J].Biodivers Sci， 12（1）：20-28.[唐志堯， 方精云， 2004.植物物種多樣性的垂直分布格局[J].生物多樣性， 12（1）：20-28.]

WANG GH， ZHOU GS， YANG LM， et al.， 2002.Distribution， species diversity and life-form spectra of plant communities along an altitudinal gradient in the northern slopes of Qilianshan Mountains， Gansu， China[J].Plant Ecol， 165（2）： 169-181.

WHITTAKER RJ， WILLIS KJ， FIELD R， 2001.Scale and species richness：Towards a general， hierarchical theory of species diversity[J].J Biogeogr， 28（4）：453-470.

WILLIS KJ， WHITTAKER RJ， 2002.Species diversity scale matters[J].Science， 295（5558）：1245-1248.

XU MH， JIA YY， ZHANG XY， et al.， 2018.Evaluation of agricultural economic development and its spatial-temporal distribution in contiguous dire poverty-stricken areas of Lüliang District in Shanxi[J].J Agric， 8（5）：73-80.[徐滿厚， 賈燕燕， 張瀟月， 等， 2018.山西呂梁連片特困區(qū)農(nóng)業(yè)經(jīng)濟(jì)發(fā)展評(píng)價(jià)及其時(shí)空分布特征[J].農(nóng)學(xué)學(xué)報(bào)， 8（5）：73-80.]

YANG XY， ZHANG SX， WEN J， et al.， 2018.Spatial pattern of herbaceous plant species diversity and its changes due to simulated warming in the forest community of the Lüliang Mountains[J].Acta Ecol Sin， 38（18）：6642-6654.[楊曉艷， 張世雄， 溫靜， 等， 2018.呂梁山森林群落草本層植物物種多樣性的空間格局及其對(duì)模擬增溫的響應(yīng)

[J].生態(tài)學(xué)報(bào)， 38（18）：6642-6654.]

YE PC， CHEN H， WU JY， et al.， 2019.Distribution pattern and correlation with main environmental factors of higher plant diversity in Northwest Yunnan[J].J Ecol Rur Environ， 35（11）：1-6.[葉鵬程， 陳慧， 武建勇， 等， 2019.滇西北高等植物多樣性分布格局及其與主要環(huán)境因子的相關(guān)性分析[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào)， 35（11）：1-6.]

ZHANG SH， XIONG KN， ZHANG Y， et al.， 2019.Response of plant community species diversity and leaf traits of dominant species to environmental factors in different grades of rocky desertification areas[J].Guihaia， 39（8）：1069-1080.[張仕豪， 熊康寧， 張俞， 等， 2019.不同等級(jí)石漠化地區(qū)植物群落物種多樣性及優(yōu)勢(shì)種葉片性狀對(duì)環(huán)境因子的響應(yīng)[J].廣西植物， 39（8）：1069-1080.]

ZHAO SQ， FANG JY， 2006.Patterns of species richness for vascular plants in Chinas nature reserves[J].Divers Distribut， 12（4）：364-372.

ZHU L， XU GQ， LI Y， et al.， 2017.Relationships among plant species diversity， biomass， and the groundwater table in the Hailiutu River basin[J].Acta Ecol Sin， 37（6）：1912-1921.[朱麗， 徐貴青， 李彥， 等， 2017.物種多樣性及生物量與地下水位的關(guān)系——以海流兔河流域?yàn)槔齕J].生態(tài)學(xué)報(bào)， 37（6）：1912-1921.]

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