祁連山高寒草甸退化對植被群落的影響

摘要：為探討高寒草甸退化對植被群落的影響，本研究以高寒草甸逐漸退化并演替為高寒草原的過程為研究背景，分析了植被群落的生境、生物量、多樣性和群落結構差異，并探索了多樣性和群落差異的驅動因子。結果表明：從高寒草甸退化為高寒草原，降水顯著減少（ ，氣溫顯著升高（ ，土壤 值顯著升高（ ？ lt; 0 . 0 5 ） ，土壤養分流失；生物量顯著下降 ;植被物種豐富度顯著下降 ，群落差異顯著 。同時，年平均降水對植被多樣性有顯著促進作用 ，而放牧對其有顯著負面影響 。此外，生境解釋了植被群落變化的 64 % ，多元回歸和層次分割法表明，放牧、海拔、降水及土壤有機質含量是植被群落變化的主要調節因子 。綜上，調控放牧壓力、加強水資源保護以及改善土壤狀況，是防止草地退化和保持生物多樣性的重要策略。

中圖分類號：S812 文獻標識碼：A 文章編號：1007-0435（2025）04-1249-09

Abstract：To investigate the impact of alpine meadow degradation on vegetation community，in this study we analyzed habitat，biomass，diversity，and community structure differences within the context of alpine meadow succession to alpine steppe.Additionally，we explored the driving factors of diversity and community differences.The results showed that as alpine meadows were degraded into alpine steppes，precipitation significantly was decreased （ ），temperature significantly was increased （ ，soil value significantly raised ），and soil nutrients were depleted. Biomass significantly was declined （ ，species richness of vegetation significantly were decreased ，and community differences became pronounced （ . Furthermore，multi-year average precipitation significantly promoted vegetation diversity （ P lt; O. O1），while grazing had a significant negative impact （ .Habitat accounted for 64 % of the changes in vegetation communities.Multiple regresson and hierarchical partitioning indicated that grazing，altitude，pre cipitation，and soil organic mater were the main factors regulating changes of vegetation communities （ P lt; 0.01）.In summary，adjusting grazing management，protecting water resources，and improving soil conditions are crucial measures to maintain biodiversity in grassland ecosystems.

Key words： Alpine meadow；Degradation；Multiple regression； Hierarchical partitioning

青藏高原作為全球重要的生物多樣性保護地區，被譽為\"世界第三極”1-2]。其中，高寒草甸是該地區物種多樣性的關鍵生態系統[3]。然而，高寒草甸正面臨顯著的穩定性挑戰，且抗干擾能力不足[4]。由于過度放牧和人類活動的干擾，高寒草甸普遍出現退化現象，目前退化草甸面積已占青藏高原自然草地總面積的三分之一[5]，導致草地的生產力下降和物種多樣性喪失，對社會經濟發展及生態可持續性造成了深遠的影響[6-7]。

大量研究已圍繞高寒草甸的退化進行了深入探討，涵蓋了其空間分布[8]、生產力、物種組成[9]、土壤性質[10]、功能多樣性、構建機制[11]以及修復措施[12等多個方面。然而，由于青藏高原地形的復雜性和空間的高度異質性，加之研究時間與技術的限制，先前的研究主要集中于小范圍的植被調查[13]，并對同類型草地的梯度退化進行了詳細論述。近年來，隨著研究的空間尺度擴大，研究者發現高寒草甸的退化不僅涉及沼澤化，還包括灌叢化的演替過程[14-15]。在全球氣候變化的背景下，青藏高原的高寒草甸正逐步向高寒草原轉變[16-18]，草甸退化的后期階段，植被的建群種開始由草甸和草原物種組成的混合群落向草原物種主導的群落演替[19]。最近，借助遙感技術和云計算技術的應用，研究揭示了過去40年間青藏高原地區高寒草甸與高寒草原之間的持續演替現象[20]，這些發現為高寒草甸退化研究提供了堅實的理論基礎。

鑒于此，明確高寒草甸到高寒草原退化演替下，植被群落的生境、生產力、物種多樣性和群落變化及其驅動因素對高寒草地生態系統的物種多樣性保護具有重要意義。本研究著重探討高寒草甸到草原退化演替下植物生境、生物量、物種多樣性以及群落變化，旨在明確從草甸到草原的演替中生境、生產力、物種多樣性和群落演變規律，并揭示驅動物種多樣性和群落變化的關鍵因子。

1 材料與方法

1.1研究區概況、植被調查和土樣采集

祁連山是青藏高原重要的組成部分，其草地生態系統占據了總面積的 57 % ，其中 7 2 . 4 % 的草地正在經歷退化[21-23]。本研究于2021年7月下旬至

8月初（生物量達到峰值的時期）在祁連山區域進行了采樣，涵蓋了 的樣帶 ， ，海拔 m）。研究主要關注高寒草甸和高寒草原兩種生態類型，共設置了14個采樣點（每種類型各7個），采樣區域間最小距離為 （圖1）。

在每個采樣區域內，隨機設立了一個 的標準地塊。在每個地塊中，進一步隨機設置15個樣方 。在每個樣方內，記錄了每種植物物種的存在情況（0或1）、豐度（樣方內植物的總數）及蓋度（采用針刺法測定）。對于現場未能鑒定的植物，采集標本并進行拍照，隨后依據中國植物圖像庫（https：//ppbc.iplant.cn/）和《中國植物志》（https：//www.iplant. 進行詳細鑒定。此外，收集了每個樣方內所有植物的地上生物量，并將其在 的條件下干燥 ，以測得干燥后的地上生物量重量。最后用直徑為 的土鉆采集 王層的土壤樣品，每個樣方采集5個樣本并混合成一個樣本，用于土壤理化性質測定。

1.2空間、氣候、海拔和放牧指標

為了充分表征采樣站點的地理結構并考慮空間自相關性，本研究將經緯度信息轉化為地理特征向量[24]。考慮到地理特征向量眾多，選擇了一個代表較大尺度地理特征的空間向量（MEM3； P lt; 0.01）進行分析。空間向量（即Moran'seigenvectormaps，MEM）的計算是通過R軟件包adspatial中的dbmem函數完成[25]。

1960一2020年間，祁連山國家公園的平均氣溫和降水量均呈現顯著的增加趨勢（ P lt; 0.05）[26]，因此可以將年平均溫度和降水量作為該區域氣候特征的代表。基于相關研究[27]，使用全球定位系統（GPS）記錄每個采樣區域的地理坐標（經緯度）和海拔。隨后，利用時間跨度為2000一2021年、空間分辨率為 的氣候柵格數據（該分辨率小于采樣區域間最小距離 ，滿足精度要求）通過ArcGIS10.8軟件提取了采樣區域的年平均氣溫和年平均降水量，以此代表研究區域的主要氣候特征。氣候柵格數據來源于國家科技基礎條件平臺—國家地球系統科學數據共享服務平臺-黃土高原科學數據中心（網址：http：//loess. geodata.cn）。

羊密度的空間分布可反映區域的放牧壓力[27]。本研究根據采樣地的經緯度，并結合實地觀察及柵格數據（包括冷或暖季牧場的信息），確定采樣地為冷季或暖季牧場。使用ArcGIS10.8軟件，基于2020年高分辨率牲畜柵格數據28（該數據得到第二次青藏高原科學考察與研究計劃的支持，空間分辨率為 ，小于采樣區域間最小距離 ，滿足精度要求），根據采樣區域經緯度提取該區域的羊、牛密度信息。通過計算天然草地的合理載畜量，將不同牲畜轉換為標準羊單位（1只牛 = 5 只綿羊單位），并用單位面積的羊單位表示放牧壓力。相關數據可通過Zenodo鏈接訪問：https：//doi.org/10.5281/zenodo.7692064。

1.3 土壤理化測定

土壤全氮（Totalnitrogen，TN）采用凱氏定氮法測定，有機質（Organicmatter，OM）通過重鉻酸鉀氧化-硫酸亞鐵滴定法進行測定，土壤pH值采用pH計（pHS-25）測定[29]。測定方法與姚世庭等[30]一致。

1.4 統計分析

首先，計算每個物種的重要值，使用以下公式[31]：

式中， 為物種 i 在群落中的重要值， 為物種 i 在群落中的相對高度； 為物種 i 在群落中的相對蓋度； 為物種 i 在群落中的相對生物量[29]。

物種多樣性采用Patrickindex（S），Shannonindex指數（H），Pielou指數（E）來表示，計算公式如下[3]：

其中 H 為多樣性指數， S 為豐富度指數， 為樣地中存在的物種數量， 為樣地中物種 i 的重要值，N 為樣地中所有物種的重要值之和， 為物種 i 的相對重要值。

環境參數、生產力和物種多樣性指標間的差異，采用非參數的Wilcoxonrank-sumtest方法進行檢驗。為探究植被群落變化，基于物種重要值和Bray-Curtis距離，進行主坐標分析（principalcoor-dinatesanalysis，PCoA）和非度量多維標度（NMDS），并使用置換多元方差分析（Permutationalmultivariateanalysisofvariance， PER-MANOVA）對群落差異性進行了檢驗，均在vegan包中進行[32]。為確定影響退化過程中物種多樣性的關鍵因子，以Patrickindex和Shannonindex為響應變量，以空間、氣候、放牧、土壤因子為預測變量，先進行mantel檢驗確定與物種多樣性相關的因子。為進一步揭示各因子對物種多樣性的具體影響，首先為了排除因子的共線性，采用方差膨脹因子（Varianceinflationfactor，VIF）進行計算，排除值大于5的因子，然后使用 包，通過全子集回歸，并判斷其AIC值來選擇最優模型[33]，預測因子對總方差的相對效應通過glmm.hp包對 進行層次分割計算獲得[34]，表達為調整后的 的百分比。為了確定影響退化過程中植被群落變化的因子，以植被群落為響應變量，以空間、氣候、放牧、土壤因子為預測變量，進行冗余分析（Redundancyanalysis，RDA）。

為了進一步確定驅動因子對植被群落的具體影響，基于Bray-curtis相異度，進行了多元回歸（Multipleregressionmatrices，MRM），隨后對MRM的 通過層次分割計算[34]，明確各因子的具體影響。隨后使用ggplot2包進行繪圖。所有統計分析和繪圖均在R4.2.3（RDevelopmentCoreTeam，2023）中完成。

2 結果與分析

2.1 生境差異

從高寒草甸（Alpinemeadow，AM）到高寒草原（Alpinesteppe，AS），植被生境指標顯示出顯著差異。具體來說，AM的年平均降水量顯著高于AS P lt; 0.05），而多年平均溫度和土壤pH值則顯著低于AS中 ，總氮和有機質含量呈現下降趨勢（表1）。結果表明，從高寒草甸到高寒草原的退化演替過程中，氣候趨向于暖干化，并伴隨著土壤養分的流失。

2.2生物量和植被物種多樣性

對高寒草甸到高寒草原的地上生物量和植被物種多樣性進行分析。結果顯示，從AM到AS，地上生物量和物種豐富度顯著降低（圖2a和2b； P lt; 0.05），均勻度顯著降低（圖 ，香農指數極顯著降低（圖2d； 。這些結果表明，從高寒草甸到高寒草原的退化演替，植被物種多樣性減少，生產力下降。

2.3植被群落差異

通過主坐標分析（PCoA）和非度量多維尺度（NMDS）分析方法，基于Bray-Curtis距離對高寒草甸與高寒草原植被群落的差異進行了分析，同時，利用基于Bray-Curtis距離的置換多元方差分析法（PERMANOVA）檢驗了群落間的顯著性差異。結果顯示，盡管AM和AS的植被群落部分重疊，但兩者間的群落差異顯著（圖 ，表明退化演替顯著影響了植被群落的結構。此外，NMDS分析同樣證實了AM與AS植被群落間存在顯著差異（圖3b； 。這些結果表明，高寒草甸的退化顯著促進了植被群落的分化，并增加了群落間的差異性。

2.4植被物種多樣性的驅動因子

對植被 α 多樣性與氣候、放牧和土壤等因子的Manteltest分析結果表明，植被 α 多樣性與空間因子（MEM3）、放牧和年平均降水顯著相關（圖 4 ； P lt; 0.05），說明植被物種多樣性受到地理空間結構、放牧和降水影響較大。

在高寒草甸向高寒草原的退化演替過程中，為了明確哪些因素驅動植被 α 多樣性的變化以及這些因素具體的影響程度，本研究采用全子集回歸方法，并通過判斷AICc值來選擇最優模型。通過對模型的調整 進行層次分割，我們詳細解析了空間、氣候、放牧和土壤因素的具體影響。結果顯示，生境因素對植被 α 多樣性變化的解釋率達到 59 . 7 % （圖5； ，其中氣候因子解釋了 3 5 . 1 5 % ，地理因子解釋了 1 8 . 0 1 % ，放牧因子解釋了 34 . 0 9 % ，土壤因子解釋了 12 . 7 5 % 。在這些因子中，年平均降水對 α 多樣性有顯著的正向影響 ，而放牧則對 α 多樣性產生了顯著的負向影響 。

2.5植被群落的驅動因子

在高寒草甸向高寒草原的退化演替過程中，為了明確驅動植被群落的變化的因子及具體影響程度，本研究，先采用冗余（RDA）分析明確了生境因子與群落的關系（圖6a），然后采用多元回歸方法，通過對模型的調整 進行層次分割，我們詳細解析了空間、氣候、放牧和土壤因素的具體影響，同時也明確了關鍵因子的具體貢獻（圖6b）。結果顯示，生境因素與植被群落顯著相關（圖6a； P lt; 0.001），進一步多元回歸表明，生境解釋了植被群落變化的 64 % （圖 ，其中氣候因子解釋了 2 6 . 8 6 % ，地理因子解釋了2 5 . 2 5 % ，放牧因子解釋了 2 2 % ，土壤因子解釋了2 5 . 8 5 % 。其中，放牧顯著解釋 $2 2 ＼% （ P { lt; } 0 . 0 1 ）$ ，海

Fig.2Biomass and vegetation species diversity 注：（a）生物量；（b）豐富度；（c）均勻度；（d）香農指數；AM為高寒草 甸；AS為高寒草原； . ？ ？ 0 . 0 0 1 lt; P lt; 0 . 0 1 . 0.001 Note：（a）Biomass；（b）Patrick index；（c）Piclou index；（d）Shannon index；AM represents alpine meadow；AS represents alpine steppe； . ？ ？ 0 . 0 0 1 lt; P lt; 0 . 0 1 ： （2

拔顯著解釋 1 9 . 7 5 % 的變異（ ，年平均降水顯著解釋 1 7 . 1 6 % 的變異（ ，土壤有機質顯著解釋 1 6 . 5 9 % 的變異（ ，年平均溫度解釋了 9 . 7 0 % 的變異 ，土壤 值解釋了9 . 2 8 % 的變異 ，空間因子解釋了 5 . 5 0 % 的變異（ 。

注：（a）基于主坐標分析；（b）基于非度量多維尺度分析；AM表示高寒草甸；AS表示高寒草原； Note：（a）Basedonpricipalcoordinatesanalysis（PCoA）；（b）Basedonnon-metricmultidimensionalsaling（NMS）；AMrepresentsalpinemeadow；AS represents alpine steppe；

3討論

在青藏高原，高寒草甸退化的后期階段，植被的建群種開始由草甸物種和草原物種組成的混合群落，逐漸演替為以草原物種為主導的群落[19]。本研究中，從高寒草甸（AM）退化到高寒草原（AS），氣候暖干化，土壤養分流失。這說明從AM到AS的退化演替，其環境向著一個更加干燥、溫度升高、土壤養分貧瘠的趨勢發展。這與高寒區植物群落傾向于退化為更干燥的生態系統，即從高山沼澤草甸到草甸、草原，沙漠的觀點一致[35-36]。降水變化是草地生產力的一個關鍵決定因素，并直接加劇草地退化[37-38]。本研究發現，從AM到AS的退化演替中，植被生產力顯著降低 ，這說明，降雨量的減少、氣溫升高、土壤養分流失、環境壓力的增加直接影響了植物的生長，限制了生態系統的生物量積累和能量流動，從而降低了整個生態系統的生產力。

草地退化是一個復雜的生態過程，其顯著特征之一是物種多樣性的喪失，這不僅減少了生態系統內的物種數量，還破壞了生物群落的結構和功能[39]。當草地退化時，原本豐富多樣的植物群落通常被少數耐旱和耐貧瘠的物種所替代，削弱了其對環境變化的抵抗力和適應能力，導致生態系統的復原力和穩定性下降[40]。因此，退化草地中的物種多樣性一直是研究的熱點領域。本研究針對從高寒草甸（AM）到高寒草原（AS）的退化演替過程，發現Patrick指數、Shannon指數和Pielou指數均顯著下降，表明植被物種多樣性和均勻度明顯降低。這一結果與本研究中AM到AS環境中降雨量的減少、溫度的升高、土壤養分的流失以及環境壓力的增加，直接影響群落生產力降低的觀察結果相對應。這與Bai等在內蒙古對不同退化程度的溫帶草原植被進行的研究相符，他們發現植被物種多樣性與退化梯度（退化加劇）呈負相關[40]。同樣，Wang等[41]也觀察到從AM到AS的植被物種多樣性和均勻度呈遞減趨勢。此外，本研究發現，在AM到AS的退化演替過程中，物種多樣性的變化受到生境因子（氣候、地理、放牧、土壤）的顯著影響 ： 。其中，年平均降水對 α 多樣性具有顯著的正向作用 ，而放牧則對 α 多樣性產生了顯著的負向影響 。這一發現與全球多地區的研究結果一致，即水分供給增加通常會增加植物可用資源，從而促進更多種類的植物生存并增加生物多樣性[42]。與此相對應，放牧通過植被的移除、土壤壓實和可能引入的入侵種等機制，直接破壞了生態系統的自然狀態，減少了物種的多樣性[43]。本研究結果強調了在管理退化草地時，應該關注調整放牧強度和改善水資源管理，以保護和恢復草地植被物種多樣性。

植被群落是物種多樣性的綜合表現，通過物種之間相互作用和對環境因子的響應，共同構建了生態系統的多樣性框架[44]。本研究中，主坐標分析（PcoA）和非度量多維尺度排序分析（NMDS）結果表明，AM群落聚集緊密，而AS群落則相對分散，且AM與AS群落間差異極顯著（ 。這說明AM與AS群落間差異極顯著，但AM群落在物種組成上展現出較高的相似度，揭示了其群落內部具有較強的一致性，這種一致性可能是相似的環境條件和生態位，使得適應該環境的物種得以繁盛[45]。相反，AS群落一致性較差，這可能反映了更為復雜和多變的環境壓力。正如本研究中發現AS環境中降水量的減少和溫度升高可能直接影響植物可利用水分，加之土壤養分的流失，進一步影響了植被群落。在較大的空間尺度上，植物群落的結構主要受氣候和母質主導[46]。此外，本研究發現，在AM到AS的退化演替過程中，植被群落的變化受到生境因子的影響 。其中放牧、海拔和多年平均降水為主要的影響因子。放牧活動顯著解釋了 2 2 % 的變化 ，放牧不僅通過直接的植被消耗改變群落組成，還可能通過土壤壓實和植被覆蓋的減少間接影響生態系統的其他組成部分。海拔顯著解釋了 1 9 . 7 5 % 的群落變化 ？ P lt; 0.01），說明海拔梯度對植物生長條件如溫度和濕度的基本限制，從而直接影響物種的生存和分布。多年平均降水則顯著解釋了 1 7 . 1 6 % 的變化 （ P lt; 0.01），降水量直接影響植物水分的可用性，進而影響植物的生長和物種間的競爭關系。土壤有機質顯著解釋了 1 6 . 5 9 % 的變化（ ，土壤有機質可以提高土壤肥力，支持更多物種的生長。這些結果強調了在草地管理和保護策略制定時，必須考慮這些關鍵的生境因子，特別是在調整放牧管理、保護水資源和改善土壤條件方面，以維持或恢復草地生態系統的健康。

4結論

本研究專注于祁連山高寒草甸向高寒草原退化的演替過程，探討了植被的物種多樣性和群落差異及驅動因子。結果表明，在退化過程中，氣候趨向暖干化，土壤養分流失，生產力和物種多樣性下降，群落結構分化加劇；降水促進了植被多樣性，而放牧則產生了負面影響；放牧、海拔、降水和土壤有機質是影響植被群落的主要因子。本研究指出，適當調控放牧強度、有效保護水資源及提升土壤質量是維護草地生態系統生物多樣性并防止草地退化的核心措施。

參考文獻

[1]MYERSN，MITTERMEIERRA，MITTERMEIERCG，etal.Biodiversity hotspots for conservation priorities[J].Nature，2000，403（6772）：853-858

[2]CHEN H，ZHU QA，PENG C H，et al. The impacts of cli-mate change and human activities on biogeochemical cycles ontheQinghai-TibetanPlateau[J].Global ChangeBiology，2013，19（10）：2940-2955

[3]董世魁，湯琳，張相鋒，等．高寒草地植物物種多樣性與功能多樣性的關系[J].生態學報，2017，37（5）：1472-1483

[4]謝飆，熊成品，劉尋續．三江源地區水土流失成因、特點及防治對策[J].中國水土保持，2007（11）：19-21

[5]董世魁，鮑永清，徐永春．青藏高寒草原高大山脈的綠帶［J].森林與人類，2020（S1）：70-77

[6]尚占環，董全民，施建軍，等．青藏高原“黑土灘\"退化草地及其生態恢復近10年研究進展——兼論三江源生態恢復問題[J]．草地學報，2018，26（1）：1-21

[7]張騫，馬麗，張中華，等．青藏高寒區退化草地生態恢復：退化現狀、恢復措施、效應與展望[J].生態學報，2019，39（20）：7441-7451

[8]李積蘭，王苑，李希來，等．黃河源區高寒退化草地空間分布特征[J].西北農業學報，2024，33（1）：108-120

[9]趙帥，楊文權，藺寶珺，等．祁連山國家公園不同退化高寒草甸植物與土壤特性研究[J].草地學報，2023，31（5）：1530-1538

[10］馬源，李林芝，張德罡，等．高寒草甸根際土壤化學計量特征對草地退化的響應[J].應用生態學報，2019，30（9）：3039-3048

[11］王丹，王孝安，郭華，等．環境和擴散對草地群落構建的影響[J].生態學報，2013，33（14）：4409-4415

[12］陳晨旭，張力，趙巍，等．中國西部退化草地對不同類型恢復措施響應的Meta分析[J].草業科學，2023，40（12）：3048-3061

[13]LIU S Y，ZHANG C，ZHANG Y H，et al.Miniaturized spec-tral imaging for environment surveillance based on UAV plat-form[C]//AOPC 2017：Optical Spectroscopy and Imaging，2017：420-429

[14］候蒙京，高金龍，葛靜，等.青藏高原東部高寒沼澤濕地動態變化及其驅動因素研究[J].草業學報，2020，29（1）：13-27

[15］崔增.群落異向演替下高寒草甸多功能性變化及對短期禁牧的響應[D].西安：中國科學院大學，2023：98-99

[16］王根緒，胡宏昌，王一博，等．青藏高原多年凍土區典型高寒草地生物量對氣候變化的響應[J]．冰川凍土，2007，29（5）：671-679

[17］王謀，李勇，白憲洲，等．全球變暖對青藏高原腹地草地資源的影響[J].自然資源學報，2004，19（3）：331-336

[18]王一博，王根緒，沈永平，等．青藏高原高寒區草地生態環境系統退化研究[J].冰川凍土，2005，27（5）：633-640

[19］王俊峰，王根緒，吳青柏．沼澤與高寒草甸退化對 和 通量的影響[J].中國環境科學，2009，29（5）：474-480

[20] WANG Y F，XUE K，HU R H，et al. Vegetation structuralshift tells environmental changes on the Tibetan Plateau over 40years[J]. Science Bulletin，2023，68（17）：1928-1937

[21］劉宥延．生態優先背景下祁連山草地\"三生功能\"權衡協同及空間格局優化研究[D].蘭州：蘭州大學，2023：1-2

[22］蔣興國，鄭杰，許登奎．祁連山山水林田湖草保護修復調查研究之二——祁連山生態環境與可持續發展存在的問題[J].邊疆經濟與文化，2018（3）：31-34

[23]晉王強，郝春旭，妙旭華，等．甘肅祁連山生態文明示范區建設路徑研究[J].環境保護，2019，47（14）：33-36

[24] BORCARD D，LEGENDRE P. All-scale spatial analysis ofecological data by means of principal coordinates of neighbourmatrices[J].Ecological Modelling，2002，153（1/2）：51-68

[25]DRAY S， BAUMAN D，BLANCHET G，et al. adespatial：Multivariate multiscale spatial analysis.R package version 0.3-7[J].Ecological Monographs，2019，2（3）

[26］彭擎．氣候變化背景下祁連山國家公園生態環境模擬及預測[D].南京：南京信息工程大學，2022：61-62

[27] SARUUL K，LI J W，NIU JM，et al. Typical steppe ecosys-tems maintain high stability by decreasing the connectionsamong recovery，resistance，and variability under high grazingpressure[J]. Science of the Total Environment，2O19，659：1146-1157

[28] ZHAN N，LIU W，YE T，et al. High-resolution livestock sea-sonal distribution data on the Qinghai-Tibet Plateau in 2020[J].Scientific Data，2023，10（1）：142

[29]劉玉禎，趙新全，董全民，等．放牧對草地生態系統結構與功能影響的研究進展[J]．草地學報，2023，31（8）：2253-2262

[30］姚世庭，蘆光新，鄧曄，等．模擬增溫對土壤真菌群落組成及多樣性的影響[J].生態環境學報，2021，30（7）：1404-1411

[31] ZHANG R，LIU W，XUH L. Effects of plateau pika（Ochotona curzoniae） on alpine meadow phytocoenosium andanalysis of the strategy it uses to expand habitat[J].Acta Eco-logica Sinica，2018，38（1）：48-52

[32] OKSANEN J， BLANCHET F G，KINDT R，et al. Package‘vegan’[J].Community EcologyPackage，Version，2013，2（9）：1-295

[33]GARCIA-PALACIOSP，GROSSN，GAITANJ，etal. Cli-mate mediates the biodiversity- ecosystem stability relationshipglobally[J]. Proceedings of the National Academy of Sciencesof the United States of America，2018，115（33）：8400-8405

[34]LAI JS，ZOU Y，ZHANG JL，et al. Generalizing hierarchicalandvariation partitioning in multiple regression and canonicalanalysesusing therdacca.hpRpackage[J].Methodsin Ecol-ogyandEvolution，2022，13（4）：782-788

[35]JINXY，JINHJ，IWAHANAG，etal.Impactsofclimate-induced permafrost degradation on vegetation：A review[J].Advances in ClimateChangeResearch，2021，12（1）：29-47

[36］秦大河，丁永建.冰凍圈變化及其影響研究——現狀、趨勢及關鍵問題[J].氣候變化研究進展，2009，5（4）：187-195

[37]ZHANGZ，ZHANGZH，HAUTIERY，etal.Effects of intra-annual precipitation patterns on grassland productivitymoder-atedbythedominantspeciesphenology[J].Frontiersin PlantScience，2023，14：1142786

[38]XUYD，DONG SK，SHENH，etal.Degradation signifi-cantlydecreased the ecosystemmultifunctionalityof threealpinegrasslands：evidences from a large-scale survey on the Qinghai-TibetanPlateau[J].Journal ofMountain Science，2O21，18（2）：357-366

[39]CUIHW，WAGGC，WANGXT，etal.The loss of above-and belowground biodiversityin degraded grasslandsdrives thedecline of ecosystemmultifunctionality[J].Applied Soil Ecol-ogy，2022，172：104370

[40]BAIXH，ZHAO WW，WANGJ，etal.Reducing plant com-munityvariabilityand improvingresilience for sustainableresto-ration of temperate grassland[J].Environmental Research，2022，207：112149

[41]WANG YC，DANGN，FENG K，et al.Grass-microbial inter-domain ecological networks associated with alpine grasslandproductivity[J].Frontiers in Microbiology，2023，14：1109128

[42]PEREIRAHM，LEADLEYPW，PROENCA，etal.Sce-nariosforglobal biodiversityin the 2lst century[J].Science，2010，330（6010）：1496-1501

[43]CARBAJAL-MORON N A，MANZANO M G，MATA-GONZALEZ R. Soil hydrology and vegetation as impacted bygoatgrazing in Vertisolsand Regosols in semi-arid shrublandsofnorthern Mexico[J].The Rangeland Journal，2O17，39（4）：363

[44]COVERDALETC，DAVIESAB. Unravellingtherelation-ship between plant diversity and vegetation structural complex-ity：A review and theoretical framework[J].Journal of Ecol-ogy，2023，111（7）：1378-1395

[45]TRISOS C H，MEROW C，PIGOT A L.The projectedtiming of abrupt ecological disruption from climate change[J].Nature，2020，580（7804）：496-501

[46]HAAFD，SIXJ，DOETTERL S. Global patternsof geo-eco-logical controls on the response of soil respiration to warming[J].Nature Climate Change，2021，11：623-627

（責任編輯閔芝智）