摘要：為探究黃土高原植物群落生物多樣性與生態(tài)系統(tǒng)功能的關(guān)系，本研究以黃土丘陵區(qū)不同環(huán)境條件下穩(wěn)定的自然植物群落為對象，采用3個物種多樣性指數(shù)（Shannon-Wiener指數(shù)、Simpson優(yōu)勢度指數(shù)、Pielou均勻度指數(shù)）和4個功能多樣性指數(shù)（FRic功能豐富度、FDiv功能趨異指數(shù)、FEve功能均勻度和FDis功能離散度），選取地上生物量（Aboveground biomass，AGB）、土壤全氮（Soil total nitrogen，STN）、土壤有機碳（Soil organic carbon，SOC）和土壤全磷（Soil total phosphorus，STP）作為生態(tài)系統(tǒng)功能指標(biāo)，運用冗余分析和全因子回歸的方法對影響生態(tài)系統(tǒng)功能的因素進行分析。結(jié)果表明：Shannon指數(shù)、STP隨降雨量和氣溫遞增呈現(xiàn)先遞增后遞減的趨勢，F(xiàn)Ric呈遞減趨勢，F(xiàn)Div，AGB，SOC，STN呈遞增趨勢（Plt;0.05）；年均降雨量Pa對生態(tài)系統(tǒng)功能的貢獻值達到21.5%，功能多樣性指數(shù)（FDiv，F(xiàn)Dis，F(xiàn)Eve）對生態(tài)系統(tǒng)功能的貢獻值要高于物種多樣性（Shannon指數(shù)），且顯著影響AGB和STP（Plt;0.05）。綜上所述，在黃土丘陵區(qū)降雨量主要影響著生態(tài)系統(tǒng)功能，功能多樣性對生態(tài)系統(tǒng)功能的影響程度比物種多樣性更大。本研究可為黃土高原生物多樣性和生態(tài)系統(tǒng)功能的恢復(fù)提供理論依據(jù)。

關(guān)鍵詞：功能多樣性；物種多樣性；生態(tài)系統(tǒng)功能；植被恢復(fù)；冗余分析

中圖分類號：Q948.15;S812.8文獻標(biāo)識碼：A文章編號：1007-0435（2023）05-1490-11

Analysis on the Relationship between Biodiversity and

Ecosystem Function in Loess Hilly Region

WEI Jia-qi ZHENG Cheng CUI Meng-ying SHEN Qiao-tian JIN Xiang-qian WEN Zhong-ming

（1.College of Grassland Agriculture， Northwest A ＆ F University， Yangling， Shaanxi Province 712100，China; 2.Research Center

of Soil and Water Conservation and Ecological Environment， Chinese Academy of Sciences and Ministry of Education， Yangling，

Shaanxi Province 712100， China）

Abstract：Under different environmental conditions in Loess Hilly Region，we analyzed three species diversity indexes （Shannon-Wiener index，Simpson index，Pielou index），four functional diversity indexes （functional richness，functional divergence，functional evenness，functional dispersion） in the stable natural plant communities，and four ecosystem function indicators （aboveground biomass，soil total nitrogen，soil organic carbon，soil total phosphorus），to examine the relationship between the biodiversity of plant communities and ecosystem function in Loess Hilly Region. The redundancy analysis and full factor regression were used to identify the factors which affect ecosystem functions. The results showed that the Shannon index and STP had a trend of increasing first and then decreasing with the increasing of precipitation and temperature（Plt;0.05）. Under different environmental conditions，F(xiàn)Ric showed a decreasing trend，while FDiv，AGB，SOC and STN an increasing trend（Plt;0.05）. According to redundancy analysis，the contribution value of average annual precipitation to ecosystem function reached to 21.5%. The contribution value of functional diversity index to ecosystem function was higher than that of species diversity，and significantly affected the aboveground biomass and soil total phosphorus. In summary，precipitation mainly affected the ecosystem function in Loess Hilly Region，and functional diversity had a greater impact on ecosystem functions than species diversity. This study could provide theoretical basis for the restoration of biodiversity and ecosystem functions in the Loess Plateau.

Key words：Functional diversity;Species diversity;Ecosystem functions;Revegetation;Redundancy analysis

黃土高原是我國水土流失和土地退化程度最為嚴重的區(qū)域［1］，退耕還林還草工程的實施顯著增加了該區(qū)植被的覆蓋面積，但植被恢復(fù)選擇的物種較為單一，雖然在數(shù)值上提升了植被覆蓋度，但過度種植單一物種使得當(dāng)?shù)匚锓N豐富度降低，群落結(jié)構(gòu)簡單且生產(chǎn)力低［2］，甚至造成土壤干燥化和水資源短缺等問題［3］，威脅到生態(tài)系統(tǒng)演化進程和生態(tài)系統(tǒng)服務(wù)［4］，特別是在干旱半干旱地區(qū)，更加難以形成良好的生態(tài)系統(tǒng)功能［5］，也嚴重阻礙植被體系生態(tài)功能的發(fā)揮［6］。

生態(tài)系統(tǒng)功能即生態(tài)系統(tǒng)的過程或性質(zhì)，其中生態(tài)系統(tǒng)過程指構(gòu)成生態(tài)系統(tǒng)的生物以及非生物因素為了達到一定的結(jié)果（如養(yǎng)分循環(huán)、初級生產(chǎn)、分解作用等）而發(fā)生的復(fù)雜相互作用［7］。生物多樣性與生態(tài)系統(tǒng)功能的關(guān)系一直是生態(tài)學(xué)研究的熱點內(nèi)容之一［8-9］，在變化的景觀中兩者之間的關(guān)系十分復(fù)雜［10-11］，其影響因素包含土地利用類型、空間尺度、環(huán)境類型和研究選擇的生態(tài)系統(tǒng)功能等［12-14］。物種多樣性和功能多樣性作為生物多樣性中的兩個重要組分，物種多樣性體現(xiàn)了區(qū)域生物資源中物種種類的豐富程度［15］，功能多樣性通過植物功能性狀的不同組合以準確反映生態(tài)系統(tǒng)中植物功能之間的差異程度和時空分布的變動范圍，反映群落中物種對所在地區(qū)環(huán)境復(fù)雜變化的主動適應(yīng)［16］。隨著人類活動的加劇和氣候的變化，生物多樣性的喪失會降低生態(tài)系統(tǒng)的功能和服務(wù)［17］，所以探究生物多樣性和生態(tài)系統(tǒng)功能之間的關(guān)系則極為重要［18］。

面對環(huán)境變化，生物多樣性是維持［19-20］和緩沖生態(tài)系統(tǒng)功能的關(guān)鍵因素［21-22］，然而環(huán)境變化也會改變物種間的相互作用，導(dǎo)致生物多樣性在不同空間尺度上發(fā)生變化［23］，因此變化的環(huán)境可以通過改變生態(tài)系統(tǒng)過程的速率直接影響生態(tài)系統(tǒng)的功能［24］，也可以通過改變生物多樣性間接影響生態(tài)系統(tǒng)功能［25］。國內(nèi)學(xué)者黃小波等［26］在云南松天然次生林開展的研究中發(fā)現(xiàn)，當(dāng)群落達到中等功能水平時，物種豐富度是對多功能性影響最大的因素；而根據(jù)熊定鵬等［27］在西藏北部羌塘高原草地的研究可知，生態(tài)系統(tǒng)多功能性與種群豐富度間存在正相關(guān)關(guān)系。因此物種多樣性確實在一定程度上影響著生態(tài)系統(tǒng)功能，但由于物種多樣性在計算時認為所有物種對群落具有同等貢獻，因此會導(dǎo)致其缺失各個物種功能性狀信息［28］。

生物多樣性對生態(tài)系統(tǒng)的影響不僅在于物種種類和物種數(shù)量，更與物種所有的功能特性有關(guān)，植物功能性狀與環(huán)境變化聯(lián)系緊密，植物群落中物種個體的功能性狀在很大程度上影響著整個群落的生態(tài)發(fā)展過程［29］，這就意味著功能多樣性根據(jù)植物性狀的變化能更有效地說明其對生態(tài)系統(tǒng)功能變動的影響，也體現(xiàn)功能多樣性能夠更深入地解釋此類生態(tài)學(xué)問題［30］。Huang等人［31］在對云南松天然次生林的多樣性研究中發(fā)現(xiàn)，功能多樣性也與生態(tài)系統(tǒng)多功能存在著正相關(guān)關(guān)系，并且功能多樣性比物種多樣性能夠更好地用于解釋整個生態(tài)系統(tǒng)功能的變化。國外也有學(xué)者對258項已發(fā)表的關(guān)于自然群落生物多樣性與生態(tài)系統(tǒng)功能的研究分析發(fā)現(xiàn)，大多數(shù)生態(tài)系統(tǒng)功能的研究側(cè)重于物種多樣性的影響，然而功能多樣性的指標(biāo)通常更能預(yù)測生態(tài)系統(tǒng)的功能［18］。因此從功能多樣性的角度去研究生態(tài)系統(tǒng)功能的變化逐漸成為新的研究方向。

針對我國黃土高原地區(qū)的生態(tài)問題一直有不少的相關(guān)研究，如植被恢復(fù)的承載力［32］、人工植被恢復(fù)對物種多樣性的影響［33］、基于功能性狀的植物對環(huán)境變化的響應(yīng)與適應(yīng)機制［34］等，對于黃土高原生物多樣性和生態(tài)系統(tǒng)功能關(guān)系的研究卻多以人工林為主［35］，在研究退耕還林工程成效的同時，我們?nèi)孕桕P(guān)注黃土高原自然穩(wěn)定植物群落的情況。為此，本研究選取黃土丘陵區(qū)結(jié)構(gòu)合理、生態(tài)功能良好的穩(wěn)定自然植物群落為研究對象，對不同環(huán)境條件下的植物群落生物多樣性和生態(tài)系統(tǒng)功能的變化及其影響關(guān)系開展研究，旨在探索黃土丘陵區(qū)自然群落物種多樣性、功能多樣性與生態(tài)系統(tǒng)功能之間的相互關(guān)系，為黃土高原生物多樣性和生態(tài)系統(tǒng)功能的恢復(fù)提供依據(jù)，對該區(qū)植被恢復(fù)與環(huán)境改善具有更重要的理論與實踐意義。

1材料與方法

1.1自然概況

本研究以典型黃土丘陵區(qū)延河流域為研究區(qū)（36°23′～37°17′N，108°45′～110°28′E），該區(qū)氣候?qū)俅箨懶约撅L(fēng)氣候，年均降水量500 mm，年均氣溫9℃。延河的干流全長286.9 km，流域面積7 687 km²，平均坡度4.4°。從東南向西北，降雨、溫度呈現(xiàn)出明顯的梯度變化特征，植被分布具有明顯的地帶性規(guī)律，依次劃分為森林區(qū)（以蒙古櫟（Quercus mongolica）、雞爪槭（Acer palmatum）、狼牙刺（Sophora viciifolia）為主）、森林草原區(qū)（以荊條（Vitex negundo）、野棗樹（Zizyphus jujuba）、白羊草（Bothriochloa ischaemun）為主）、典型草原區(qū)（以鐵桿蒿（Artemisia sacrorum）、達烏里胡枝子（Lespedeza daurica）、百里香（Thymus mongolicus）、長芒草（Stipa bungeana）為主）。流域土壤類型主要為黃綿土，具有土質(zhì)疏松，抗侵蝕能力差等特點。

1.2取樣方法

為有效保證樣品的采集及其在空間分布上具有代表性，在黃土高原丘陵區(qū)根據(jù)水熱條件形成的三個植被帶上分別選取了完整斷面（陰坡、陽坡）進行野外調(diào)查，均為人為干擾較輕、自然植被發(fā)育較好、具有代表性的樣地，采樣時間為2021年7—8月（生長旺盛季）。每個斷面包含陰、陽坡，在不同坡位（上坡、中坡、下坡）分別設(shè)置采樣點，每個調(diào)查樣點內(nèi)設(shè)置獨立的3個1 m×1 m的草本樣方、1個5 m×5 m的灌木樣方和1個10 m×10 m的喬木樣方，分別進行了植物群落調(diào)查，總計108個樣方。草本樣方調(diào)查物種名、物種數(shù)、株數(shù)、高度、蓋度等，將草本樣方中全部地上植物部分剪下來，裝入樣品袋并做好標(biāo)記，帶回實驗室稱量鮮重（g），用以計算生物量；灌木樣方調(diào)查植株物種名、物種數(shù)、高度、冠幅等；喬木樣方調(diào)查樹木物種名、物種數(shù)、樹高、胸徑、地徑等。使用GPS實際測量每個樣點的經(jīng)緯度、海拔，使用坡度儀測量坡度。土壤采樣樣點布置與植被樣點相同，用土鉆法進行取樣，剖面按垂直高度0～10 cm，10～20 cm，20～40 cm進行分層，每層按照“S”型進行3次重復(fù)取樣后充分混合，裝入塑封袋后編號帶回室內(nèi)，將其過0.149 mm篩，去除植物根系和石塊等，再使其自然風(fēng)干，用以土壤養(yǎng)分的測定。

1.3植物功能性狀的計測

本研究選取了6種植物功能性狀，分別為：植物高度（Plant height，H）、比葉面積（Specific leaf area，SLA）、葉組織密度（Leaf tissue density，LTD）、葉氮含量（Leaf nitrogen content，LN）、葉碳含量（Leaf carbon content，LC）和葉磷含量（Leaf phosphorus content，LP）。樣方內(nèi)出現(xiàn)的每個物種隨機選擇3—5株測量其株高（H）；采集每個樣方內(nèi)每種完全伸展且沒有病蟲害的物種葉片10～20片，用掃描儀獲取植物的平展葉片圖像，使用分析軟件Image-J計算葉面積；葉厚度選用精度為0.01 mm的電子游標(biāo)卡尺；葉片在105℃下殺青15 min，85℃下烘干48～72 h后稱重。比葉面積=葉面積/葉干重，葉組織密度=葉片干重/葉片體積，葉片體積=葉面積×葉厚度。葉氮含量（LN）、葉碳含量（LC）和葉磷含量（LP）分別用凱氏定氮法、重鉻酸鉀外加熱法和鉬銻抗比色法測定，單位為g·kg^-1。植物群落性狀值描述性統(tǒng)計見表2。

1.4環(huán)境氣象數(shù)據(jù)的獲取

環(huán)境氣象數(shù)據(jù)是用Arcgis 10.2軟件根據(jù)地形圖、數(shù)字高程圖像、延河流域及周邊57個氣象站的站點觀測數(shù)據(jù)進行匯總（2021年），通過空間插值的方法得到延河流域的環(huán)境因子專題圖像［36］，再結(jié)合野外記錄的各樣點的經(jīng)緯度信息從現(xiàn)有的環(huán)境因子專題圖中提取得到，包括年均降雨量（Average annual precipitation，Pa）、7—9月均降雨量（Average annual precipitation from July to September，P789）、年均溫（Annual mean temperature，Ta）、4—10月平均溫度（Average temperature from April to October，T410）。

1.5生物多樣性指數(shù)的計算

1.6生態(tài)系統(tǒng)功能指標(biāo)的選取和測定結(jié)合以往相關(guān)研究，選取出以下4種生態(tài)系統(tǒng)功能指標(biāo)來量化生態(tài)系統(tǒng)功能［37-38］。地上生物量（Aboveground biomass，AGB）用于表征植物生長功能，也是研究生態(tài)系統(tǒng)功能的關(guān)鍵；土壤有機碳（Soil organic carbon，SOC）含量用于表征土壤有機碳蓄積功能；土壤全氮（Soil total nitrogen，STN）含量、土壤全磷（Soil total phosphorus，STP）含量用于表征土壤養(yǎng)分循環(huán)功能。

將在野外采集的草本植物樣品放入75℃的烘箱中烘24 h，再稱量其干重；灌木的地上生物量根據(jù)其地徑（Dt）、株高（H）、冠幅（C）、冠幅面積（A=π（C/2）²）、冠幅體積（V=A×H），參考羅永開等［39］關(guān)于灌木生物量模型總結(jié)的最佳方程和參數(shù)計算灌木的地上生物量，例如計算黃刺玫（Rosa xanthina）地上生物量的最佳擬合方程為：W=aV^b，其中參數(shù)a為0.003，參數(shù)b為0.88；喬木的地上生物量參考魏勝利［40］關(guān)于喬木生物量模型的研究，通過樹高（H）、胸徑（D）與生物量的模型方程計算喬木生物量，例如蒙古櫟（Quercus mongolica）的最佳擬合方程為：W=a×D^b×H^c，其中參數(shù)a為0.01，參數(shù)b為1.5，參數(shù)c為1.71。土壤有機碳采用重鉻酸鉀外加熱法測定，土壤全氮采用凱氏定氮法測定，土壤全磷使用鉬銻抗比色法測定。

1.7數(shù)據(jù)分析與處理

采用Excel 2021對所有數(shù)據(jù)進行計算，利用R 4.1.3中的vegan、FD軟件包計算群落的物種多樣性和功能多樣性，并用SPSS 26.0進行如下統(tǒng)計分析：（1）不同樣地物種多樣性、功能多樣性、生態(tài)系統(tǒng)功能的單因素方差分析（One-Way ANOVA）并用LSD后置檢驗法檢驗；（2）物種多樣性、功能多樣性與生態(tài)系統(tǒng)功能的逐步回歸分析；環(huán)境因子、生物多樣性對生態(tài)系統(tǒng)功能的影響采用RDA冗余分析（Canoco 5）。繪圖軟件采用Origin 2022b。

2結(jié)果與分析

2.1植物群落物種多樣性隨環(huán)境變化的響應(yīng)

隨著樣地環(huán)境的變化，Shannon指數(shù)、Simpson指數(shù)和Pielou指數(shù)呈現(xiàn)有不同程度的變化。其中，Shannon指數(shù)呈現(xiàn)先增加后減少的趨勢，EN3顯著高于EN2和EN5（Plt;0.05），EN1與其他四個樣地?zé)o顯著差異；EN5的Simpson指數(shù)顯著低于EN1～EN4 （Plt;0.05），但EN1～EN4之間無顯著性差異；Pielou指數(shù)在EN4出現(xiàn)最小值（0.72），顯著低于EN2，EN3和EN5（Plt;0.05）（圖1）。

2.2植物群落功能多樣性隨環(huán)境變化的響應(yīng)

由圖2可知，不同環(huán)境條件下植物群落功能多樣性的變化。功能豐富度（FRic） 隨著樣地環(huán)境的變化呈遞減趨勢，EN1顯著高于其他樣地（Plt;0.05）；而功能趨異指數(shù)（FDiv）呈遞增趨勢，EN4，EN5顯著高于EN1，EN2（Plt;0.05）；功能均勻度和功能離散度的變化總體較為平緩，不同樣地間差異性相對不顯著。

2.3生態(tài)系統(tǒng)功能隨環(huán)境變化的響應(yīng)

隨著不同樣地環(huán)境條件的變化，本研究選取的生態(tài)系統(tǒng)功能指標(biāo)呈現(xiàn)出不同的響應(yīng)（圖3）。EN1的地上生物量（AGB）顯著低于其他四個樣地（Plt;0.05）；土壤全氮（STN）和全碳（SOC）呈現(xiàn)相同的變化趨勢，即EN5的STN和SOC最高且顯著高于其他四個樣地（Plt;0.05），EN3顯著高于EN1，EN2（Plt;0.05），EN1，EN2和EN4之間無顯著差異；EN2，EN4和EN5的土壤全磷（STP）間無顯著差異，但均顯著高于EN1，EN3（Plt;0.05），而EN1和EN3的STP不存在顯著差異性。

2.4環(huán)境因子、生物多樣性與生態(tài)系統(tǒng)功能的關(guān)系

以生態(tài)系統(tǒng)功能指標(biāo)作為響應(yīng)變量，環(huán)境因子和生物多樣性作為解釋變量進行冗余分析（RDA）。由圖4可以表明，第1排序軸的解釋量為20.41%，第2排序軸為16.16%。分析箭頭長短、夾角大小及冗余分析的相關(guān)系數(shù)可以得出，環(huán)境因子和生物多樣性對生態(tài)系統(tǒng)功能的影響大小為：Pagt;FDivgt;Shannongt;FDisgt;FDvegt;Tagt;T410（表5排除貢獻值lt;1%的解釋變量）。年均降雨量Pa與AGB，STN，SOC，STP均呈顯著正相關(guān)關(guān)系（Plt;0.05），解釋量為10.0%，貢獻值達到21.5%（表5）；Shannon指數(shù)與SOC，STN呈正相關(guān)關(guān)系，與AGB，STP呈負相關(guān)關(guān)系；而功能多樣性指數(shù)中FDiv，F(xiàn)Dve及年均溫Ta與AGB，STP呈正相關(guān)關(guān)系，與SOC，STN呈負相關(guān)關(guān)系；FDis與STN，SOC，STP呈正相關(guān)關(guān)系，與AGB呈負相關(guān)關(guān)系。

全因子回歸分析表明（表6），AGB，STP與本研究的多樣性指數(shù)達到顯著水平（Plt;0.05）。根據(jù)表7中的全因子回歸分析系數(shù)得到生態(tài)系統(tǒng)功能與多樣性指數(shù)的相關(guān)方程為Y₁=-13 727.064+8 502.281X₅+28 083.649X₆-6 099.455X_７和Y₄=0.360+0.181X₆。由方程Y₁=-13 727.064+8 502.281X₅+28 083.649X₆-6 099.455X₇可知（表7），隨著功能均勻度和功能趨異指數(shù)的增加，AGB呈增加趨勢，而隨功能離散度的降低而減少；從方程Y₄=0.360+0.181X₆可知（表7），STP隨著功能趨異指數(shù)升高而顯著增加。

3討論

3.1植物群落物種多樣性對環(huán)境的響應(yīng)變化

隨著環(huán)境條件的變化，Shannon指數(shù)呈現(xiàn)先增加后減少的趨勢（圖1），這是由于EN1～EN5的降雨量在逐漸增加，EN1樣地主要為草本植物，而EN2～EN5存在灌木和喬木，特別是EN3和EN4的灌木植物種類更為豐富。因此降雨量和溫度的增加，使得物種的組成更加豐富，Shannon指數(shù)在EN1～EN3出現(xiàn)增加趨勢；在干旱半干旱地區(qū)，物種多樣性差異不大的樣地間，物種種類卻在發(fā)生變化，灌木和喬木會使草本植物能夠獲取的資源變少，因此在EN3和EN4這兩個環(huán)境中，中華卷柏（Selaginella sinensis）等耐旱植物的數(shù)量逐漸增多，這說明由于水資源的限制，耐旱植物在植物群落分布中的數(shù)量會逐漸增多，從而占據(jù)優(yōu)勢地位［41］。EN5物種多樣性的降低是由于該樣地主要優(yōu)勢種是喬木，灌木和草本植物的種類較少，因此物種多樣性呈現(xiàn)出減少的趨勢。Pielou均勻度指數(shù)沒有呈現(xiàn)明顯變化，這是因為本研究中選取的是無人為干擾的穩(wěn)定自然草地，因此生態(tài)系統(tǒng)大多處于穩(wěn)定狀態(tài)，能夠在不同環(huán)境條件下穩(wěn)定存在，也說明在大尺度的環(huán)境變化下，群落物種組成發(fā)生的改變即使會造成群落生態(tài)位發(fā)生較大變動，也可能繼續(xù)有新的物種出現(xiàn)來重新維持生態(tài)系統(tǒng)結(jié)構(gòu)的平衡［42］。

3.2植物群落功能多樣性對環(huán)境的響應(yīng)變化

植物性狀對于環(huán)境變化的反饋體現(xiàn)了物種對環(huán)境的適應(yīng)性，因此功能多樣性直接反映了群落在整個生態(tài)系統(tǒng)中主要的功能分布范圍［43］。本研究發(fā)現(xiàn)，EN1植物群落的功能豐富度顯著高于其他4個樣地（圖2），功能豐富度即是指物種在群落中所占據(jù)的功能空間的大小，體現(xiàn)著生態(tài)空間資源的利用程度［44］，這說明EN1群落中草本植物能夠通過功能性狀的改變對生態(tài)系統(tǒng)中的資源進行更加充分的利用，因此在EN1群落中生態(tài)位空間資源的利用程度明顯較高。另外，功能趨異指數(shù)主要體現(xiàn)了群落內(nèi)各物種的生態(tài)位分化情況和空間資源競爭程度［44］，從圖2功能趨異指數(shù)的變化可以發(fā)現(xiàn)，EN1和EN2群落內(nèi)植物生態(tài)位互補程度都較低，之間的競爭較強，而EN3～EN5群落中喬灌草的空間分布使得其生態(tài)位互補程度更強，競爭就會較弱。本研究中功能均勻度和功能離散度在不同環(huán)境梯度間差異不顯著（圖2），雖然降雨量和溫度存在差異，但本研究中的樣地均為穩(wěn)定的自然植物群落，生態(tài)系統(tǒng)處于穩(wěn)定狀態(tài)，因此物種功能性狀分布均勻，植物對外界干擾具有防御能力［30］。

3.3生態(tài)系統(tǒng)功能對環(huán)境變化的響應(yīng)

不同樣地環(huán)境的變化影響著生態(tài)系統(tǒng)功能，由圖3可知，AGB，SOC，STN均呈現(xiàn)隨樣地的氣溫和降水量的增加而增大的趨勢。生物量代表了生態(tài)系統(tǒng)功能的最基本特征，EN1由于地理位置溫度和水分的限制，其植被類型主要為草本植物，因此地上生物量顯著低于其他四個樣地；而整體來看，EN2～EN5均存在有木本層和草本層，但環(huán)境條件的不同導(dǎo)致其存在不同的群落垂直結(jié)構(gòu)，在水分含量較低的生態(tài)系統(tǒng)中，草本層發(fā)揮重要的作用，而在環(huán)境資源充足的樣地（EN5），木本層則是群落生態(tài)系統(tǒng)功能變化的主要因素。這與實際調(diào)查中黃土丘陵區(qū)植被分布情況相符合：黃土丘陵區(qū)延河流域從東南向西北，降雨、溫度呈現(xiàn)出明顯的梯度變化特征，植被分布具有明顯的地帶性規(guī)律，依次為森林區(qū)、森林草原區(qū)和草原區(qū)。Tilman［45］研究也發(fā)現(xiàn)群落中的生物及其功能特性通過資源利用的互補性對生態(tài)系統(tǒng)功能產(chǎn)生影響，即高多樣性的群落中整體可利用資源的比例會增加，生態(tài)系統(tǒng)功能也會增加［4，46］。因此群落的組成是生態(tài)系統(tǒng)功能的重要影響因素，草本層和木本層對環(huán)境適應(yīng)策略的不同會導(dǎo)致其生態(tài)系統(tǒng)功能的差異。

3.4環(huán)境因子、生物多樣性對生態(tài)系統(tǒng)功能的影響

冗余分析結(jié)果顯示，年均降雨量Pa對生態(tài)系統(tǒng)功能的解釋量最高為10.0%，貢獻值達到21.5%，且Pa與AGB，STN，SOC，STP均呈顯著正相關(guān)關(guān)系（Plt;0.05），說明本研究中降雨量是黃土丘陵區(qū)生態(tài)系統(tǒng)功能主要的氣候驅(qū)動因素，這也符合水分條件是黃土高原生態(tài)系統(tǒng)功能的限制因子。有研究發(fā)現(xiàn)，在生產(chǎn)力較低的地區(qū)（如更寒冷、更干燥等），積極的生物多樣性效應(yīng)會變得更弱［47-48］；也有研究表明，雖然生物多樣性在各種情況下對許多類型的生態(tài)系統(tǒng)功能產(chǎn)生積極影響［49-51］，但在大多數(shù)情況下，其影響不如非生物因素的影響強烈［18］。盡管如此，我們考慮生物多樣性對生態(tài)系統(tǒng)功能的影響時可以發(fā)現(xiàn)，植物群落功能多樣性指數(shù)（FDiv，F(xiàn)Dis，F(xiàn)Eve）對生態(tài)系統(tǒng)功能的貢獻值要高于物種多樣性（Shannon指數(shù)），并且通過逐步回歸分析發(fā)現(xiàn)，植物群落功能多樣性顯著影響著AGB和STP，說明相比植物群落物種變化對生態(tài)系統(tǒng)功能的改變，植物群落功能多樣性即植物性狀的改變對生態(tài)系統(tǒng)功能的影響更大。這也表明在黃土丘陵區(qū)，生態(tài)系統(tǒng)功能的影響因素除了環(huán)境因子之外，較大程度取決于所研究植物群落的功能多樣性，相較于物種層面的改變，物種本身功能性狀對于環(huán)境的適應(yīng)策略也是影響生態(tài)系統(tǒng)功能的原因。

4結(jié)論

本研究綜合分析了黃土丘陵區(qū)植物群落物種多樣性、功能多樣性和生態(tài)系統(tǒng)功能隨環(huán)境條件變化的響應(yīng)，通過冗余分析、全因子回歸的方法，綜合分析了環(huán)境因子、生物多樣性與生態(tài)系統(tǒng)功能的關(guān)系。結(jié)果表明，植物群落物種多樣性、功能多樣性、生態(tài)系統(tǒng)功能在不同環(huán)境條件下存在顯著差異，年均降雨量Pa主要影響著本研究中的生態(tài)系統(tǒng)功能且呈顯著正相關(guān)關(guān)系，另外功能多樣性與AGB和STP的存在顯著相關(guān)性，相比物種多樣性對生態(tài)系統(tǒng)功能的貢獻值更大。

參考文獻

［1］FU B J，WANG S，LIU Y，et al. Hydrogeomorphic Ecosystem Responses to Natural and Anthropogenic Changes in the Loess Plateau of China［J］. Annual Review of Earth and Planetary Sciences，2017，45（1）：223-243

［2］CHEN Y P，WANG K B，LIN Y S，et al. Balancing green and grain trade［J］. Nature Geoscience，2015，8（10）：739-741

［3］FENG X M，F(xiàn)U B J，PIAO S L，et al. Revegetation in China's Loess Plateau is approaching sustainable water resource limits［J］. Nature Climate Change，2016，6（11）：1019

［4］DIAZ S，CABIDO M. Vive la difference：plant functional diversity matters to ecosystem processes［J］. Trends in Ecology amp; Evolution，2001，16（11）：646-655

［5］張超，劉國彬，薛萐，等. 黃土丘陵區(qū)不同植被類型根際微生物群落功能多樣性研究［J］. 草地學(xué)報，2015，23（4）：710-717

［6］韓慶功，彭守璋. 黃土高原潛在自然植被空間格局及其生境適宜性［J］. 水土保持學(xué)報，2021，35（5）：188-193

［7］馮劍豐，李宇，朱琳. 生態(tài)系統(tǒng)功能與生態(tài)系統(tǒng)服務(wù)的概念辨析［J］. 生態(tài)環(huán)境學(xué)報，2009，18（4）：1599-1603

［8］TILMAN D，KNOPS J，WEDIN D，et al. The influence of functional diversity and composition on ecosystem processes［J］. Science，1997，277（5330）：1300-1302

［9］HILLEBRAND H，MATTHIESSEN B. Biodiversity in a complex world：consolidation and progress in functional biodiversity research［J］. Ecology Letters，2009，12（12）：1405-1419

［10］ULRICH B，HELMUT H. Biodiversity and ecosystem functioning in dynamic landscapes［J］. Philosophical Transactions of the Royal Society B-Biological Sciences，2016，371（1694）：267

［11］SCHWARTZ M W，BRIGHAM C A，HOEKSEMA J D，et al. Linking biodiversity to ecosystem function：implications for conservation ecology［J］. Oecologia，2000，122（3）：297-305

［12］BARNES A D，WEIGELT P，JOCHUM M，et al. Species richness and biomass explain spatial turnover in ecosystem functioning across tropical and temperate ecosystems［J］. Philosophical Transactions of the Royal Society B-Biological Sciences，2016，371（1694）：279

［13］MINDEN V，SCHERBER C，PIQUERAS M A C，et al. Consistent drivers of plant biodiversity across managed ecosystems［J］. Philosophical Transactions of the Royal Society B-Biological Sciences，2016，371（1694）：284

［14］FLORES-MORENO H，REICH P B，LIND E M，et al. Climate modifies response of non-native and native species richness to nutrient enrichment［J］. Philosophical Transactions of the Royal Society B-Biological Sciences，2016，371（1694）：273

［15］張建貴，王理德，姚拓，等. 祁連山高寒草地不同退化程度植物群落結(jié)構(gòu)與物種多樣性研究［J］. 草業(yè)學(xué)報，2019，28（5）：15-25

［16］CADOTTE M W，CAVENDER-BARES J，TILMAN D，et al. Using Phylogenetic，F(xiàn)unctional and Trait Diversity to Understand Patterns of Plant Community Productivity［J］. Plos One，2009，4（5）：5695

［17］徐煒，馬志遠，井新，等. 生物多樣性與生態(tài)系統(tǒng)多功能性：進展與展望［J］. 生物多樣性，2016，24（1）：55-71

［18］VAN DER PLAS F. Biodiversity and ecosystem functioning in naturally assembled communities［J］. Biological Reviews，2019，94（4）：1220-1245

［19］CARDINALE B J，DUFFY J E，GONZALEZ A，et al. Biodiversity loss and its impact on humanity ［J］. Nature，2012，486（7401）：59-67

［20］HOOPER D U，ADAIR E C，CARDINALE B J，et al. A global synthesis reveals biodiversity loss as a major driver of ecosystem change［J］. Nature，2012，486（7401）：105-108

［21］HISANO M，SEARLE E B，CHEN H Y H. Biodiversity as a solution to mitigate climate change impacts on the functioning of forest ecosystems［J］. Biological Reviews，2018，93（1）：439-456

［22］PIRES A P F，SRIVASTAVA D S，F(xiàn)ARJALLA V F. Is Biodiversity Able to Buffer Ecosystems from Climate Change？ What We Know and What We Don’t［J］. Bioscience，2018，68（4）：273-280

［23］BLOWES S A，SUPP S R，ANTAO L H，et al. The geography of biodiversity change in marine and terrestrial assemblages［J］. Science，2019，366（6463）：339-345

［24］SPAAK J W，BAERT J M，BAIRD D J，et al. Shifts of community composition and population density substantially affect ecosystem function despite invariant richness［J］. Ecology Letters，2018，21（3）：1315-1324

［25］GILING D P，BEAUMELLE L，PHILLIPS H R P，et al. A niche for ecosystem multifunctionality in global change research［J］. Global Change Biology，2019，25（3）：763-774

［26］黃小波，李帥鋒，蘇建榮，等. 云南松天然次生林物種豐富度與生態(tài)系統(tǒng)多功能性的關(guān)系［J］. 生物多樣性，2017，25（11）：1182-1191

［27］熊定鵬，趙廣帥，武建雙，等. 羌塘高寒草地物種多樣性與生態(tài)系統(tǒng)多功能關(guān)系格局［J］. 生態(tài)學(xué)報，2016，36（11）：3362-3371

［28］董六文，任正煒，張蕊，等. 功能多樣性比物種多樣性更好解釋氮添加對高寒草地生物量的影響［J］. 植物生態(tài)學(xué)報，2022，46（8）：871-881

［29］孟婷婷，倪健，王國宏. 植物功能性狀與環(huán)境和生態(tài)系統(tǒng)功能［J］. 植物生態(tài)學(xué)報，2007（1）：150-165

［30］劉旻霞，張國娟，南笑寧，等. 甘南高寒草甸坡向梯度對植物群落功能多樣性的影響［J］. 西北植物學(xué)報，2020，40（8）：1414-1423

［31］HUANG X B，SU J R，LI S F，et al. Functional diversity drives ecosystem multifunctionality in a Pinus yunnanensis natural secondary forest［J］. Scientific Reports，2019，9（1）：6979

［32］高海東，龐國偉，李占斌，等. 黃土高原植被恢復(fù)潛力研究［J］. 地理學(xué)報，2017，72（5）：863-874

［33］巫翠華. 刺槐人工林群落林下植物多樣性研究［D］.楊凌：西北農(nóng)林科技大學(xué)，2021：18-19

［34］龔時慧，溫仲明，施宇. 延河流域植物群落功能性狀對環(huán)境梯度的響應(yīng)［J］. 生態(tài)學(xué)報，2011，31（20）：6088-6097

［35］王楠. 寧南山區(qū)檸條人工林生物多樣性和生態(tài)系統(tǒng)多功能性的關(guān)系研究［D］.楊凌：西北農(nóng)林科技大學(xué)，2021：23-26

［36］赫曉慧，溫仲明. 小流域地形因子影響下的土壤水分空間變異性研究［J］. 水土保持研究，2008（2）：80-83

［37］VALENCIA E，GROSS N，QUERO J L，et al. Cascading effects from plants to soil microorganisms explain how plant species richness and simulated climate change affect soil multifunctionality［J］. Global Change Biology，2018，24（12）：5642-5654

［38］劉旻霞，張國娟，李亮，等. 甘南高寒草甸海拔梯度上功能多樣性與生態(tài)系統(tǒng)多功能的關(guān)系［J］. 應(yīng)用生態(tài)學(xué)報，2022，33（5）：1291-1299

［39］羅永開，方精云，胡會峰.山西蘆芽山14種常見灌木生物量模型及生物量分配［J］. 植物生態(tài)學(xué)報，2017，41（1）：115-125

［40］魏勝利. 小興安嶺主要闊葉樹種樹干生物量模型研究［J］. 防護林科技，2018（5）：24-26

［41］張欽弟，衛(wèi)偉，陳利頂，等. 黃土高原草地土壤水分和物種多樣性沿降水梯度的分布格局［J］. 自然資源學(xué)報，2018，33（8）：1351-1362

［42］南笑寧. 甘南高寒草甸不同坡向植物群落功能多樣性研究［D］.蘭州：西北師范大學(xué)，2020：26-27

［43］朱云云，王孝安，王賢，等. 坡向因子對黃土高原草地群落功能多樣性的影響［J］. 生態(tài)學(xué)報，2016，36（21）：6823-6833

［44］VILLéGER S，MASON N W H，MOUILLOT D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology［J］. Ecology，2008，89（8）：2290-2301

［45］TILMAN D. Distinguishing between the effects of species diversity and species composition［J］. Oikos，1997，80（1）：185

［46］張中華，馬麗，周秉榮，等. 高寒草甸優(yōu)勢種功能多樣性對增溫和模擬放牧的響應(yīng)［J］. 草地學(xué)報，2021，29（S1）：225-232

［47］RATCLIFFE S，WIRTH C，JUCKER T，et al. Biodiversity and ecosystem functioning relations in European forests depend on environmental context［J］. Ecology Letters，2017，20（11）：1414-1426

［48］WU X，WANG X P，TANG Z Y，et al. The relationship between species richness and biomass changes from boreal to subtropical forests in China［J］. Ecography，2015，38（6）：602-613

［49］袁小強，王敏，張鮮花，等. 烏魯木齊周邊典型荒漠草地群落物種多樣性及生物量關(guān)系研究［J］.草地學(xué)報，2022，30（5）：1253-1262

［50］周宸宇，楊曉淵，邵新慶，等. 不同退化程度高寒草甸植物物種多樣性與生態(tài)系統(tǒng)多功能性關(guān)系［J］.草地學(xué)報，2022，30（12）：3410-3422

［51］王曉芬，馬源，張格非，等. 高寒草甸退化階段植物群落多樣性與系統(tǒng)多功能性的聯(lián)系［J］.草地學(xué)報，2021，29（5）：1053-1060

（責(zé)任編輯 劉婷婷）