九峰山不同林分類型生態恢復植被-土壤系統耦合關系評價

王皓月，郭月峰*，徐雅潔，祁偉, ，卜繁靖，祁慧娟

1. 內蒙古農業大學沙漠治理學院，內蒙古 呼和浩特 010010；2. 內蒙古水利水電勘測設計院，內蒙古 呼和浩特 010020

九峰山是陰山山脈的主體—大青山的主峰，最高峰海拔為2338 m（李利平，1993）。該地區植被類型豐富，物種多樣，植被保存比較完整。近年來，隨著旅游業和采礦業的興盛，對當地生態環境造成巨大破壞，如何快速有效地進行植被恢復是促進當地生態環境建設的首要問題。

林下植被是構成人工林生態系統的重要組成部分，也是人工林生態系統中物種多樣性的主要組成部分。越來越多的研究表明，林下植被對維持人工林生態系統結構和功能發揮著至關重要的作用，尤其在維持人工林生態系統穩定性、調控地上地下養分循環及能量轉化、促進生物多樣性發展等方面發揮著不可替代的作用（費玲等，2016；溫晶等，2019）。土壤理化性質是影響土壤肥力的重要因素，也是土壤質量組成的綜合反映（馬國飛等，2017），決定著林地生產力的大小。植被類型差異是影響理化性質最直接的因素（劉俊廷，2020）。一方面植被重建可以改良土壤結構，提高土壤肥力；另一方面，土壤理化性質的好壞又可以直接反作用于植被的生長發育。土壤-植被兩系統之間存在著天然的耦合關系，耦合度能有效地反映生態恢復過程中植被、土壤系統間的互動關系，評判兩者間的耦合狀態，耦合度模型分析有助于我們更深刻地認知植物-土壤反饋在全球環境變化下對生態系統影響的貢獻能力（李霖，2019）。近年來學者關于九峰山生態環境的研究主要集中于景觀資源評價、植被區系組成分析等（湯佳，2016），而對于土壤-植被系統耦合關系的研究較為缺乏。因此本研究利用分析九峰山不同林分類型土壤-植被耦合關系，以期為揭示九峰山生態恢復過程中植被和土壤之間存在的相互作用機制提供科學依據。

1 研究方法

1.1 研究區概況

九峰山地處陰山山脈中段，位于包頭市土默特右旗境內，110°25′—110°48′N，40°34′—40°51′N 之間。研究區大部分是山地，且多為25—45°的陡坡，九峰山山系為東西走向，整個山地由中山與低山丘陵組成，南低北高，相對高度為800—1000 m，為明顯的構造斷塊地形。據內蒙古氣候帶和氣候區的劃分，九峰山處于冬季溫帶和溫帶氣候帶，東勝氣候區。該區屬大陸性季風氣候即冬季受西伯利亞寒流侵入夏季受東南季風的影響，年均氣溫5.2 ℃，年降水量380—435 mm，無霜期117 d。九峰山的土壤在內蒙古自治區土壤分區中，屬于山地栗鈣土帶、山地淡栗鈣土亞帶。九峰山作為半干旱地區內，中山山體的一部分，是暖溫型和中溫型草原的天然屏障。九峰山植被類型豐富，其中森林植被包括針葉林和闊葉林，由油松（Pinus tabuliformis Carriere）、白扦（Picea meyeri Rehder & E. H. Wilson）、遼東櫟（Quercus wutaishanica Blume）、白樺（Betula platyphylla Suk.）構成。

1.2 試驗設計和測定方法

本研究于2020年8月7—18日，在內蒙古包頭市九峰山選取6種不同類型人工林地，分別為：側柏（Platycladus orientalis (L.) Franco）、檸條（Caragana korshinskii Kom.）、黃刺玫（Rosa xanthina Lindl.）、虎榛子（Ostryopsis davidiana Decaisne）、長梗扁桃（Amygdalus pedunculata Pall）、山楊（Populus davidiana Dode），在每種林分中設置1塊調查樣地，面積為1500 m2（表1）。為減少地形、坡度、坡向以及生態演替等因素的影響，選取樣地均為相同或相似的立地條件，并在樣地附近選取撂荒地作為對照。本次研究共分為兩個子目標層，分別為土壤綜合子系統與植物綜合子系統。其中土壤綜合子系統包含的下級指標為：土壤含水量、土壤容重、土壤孔隙度、土壤全氮、土壤全磷、土壤全鉀以及土壤有機碳。在采樣過程中考慮到結果的準確性，每種林分模式選取3塊樣地，并分別沿樣地對角線布設3個典型樣方（喬木樣方20 m×20 m、灌木樣方5 m×5 m、草本樣方1 m×1 m）進行植被調查，并利用剖面法和環刀法采集原狀土壤樣品，考慮到九峰山土層較薄，且人工林栽植年限較短，土壤水分變化深度較淺，故采樣深度為60 cm，分3層，每20 cm為一層，每層取樣兩次。土壤樣品經剔除石塊和動植物殘體等雜質，風干，磨碎，過篩后用于測定土壤的理化性質。土壤水分采用烘干法測定、土壤容重采用環刀法測定、全氮采用半微量凱式法、全磷采用NaOH熔融-鉬銻抗比色法、全鉀采用NaOH熔融，火焰光度計法、有機碳采用重鉻酸鉀外加熱法測定。

表1 樣地信息Table 1 Basic information of sampling plots

1.3 指標體系構建及權重計算

為了準確評價自然恢復過程中不同林分土壤和植被的關系，在設置評價指標體系時，借鑒相關研究成果（彭晚霞等，2011；李豪等，2019），按照目的科學性和系統完整性原則，兼顧指標體系的可操作性、易比較性、代表性和協調性，構建能夠反映不同人工林分類型生態恢復植被土壤耦合協調評判的2級次指標。本研究共選取植被指標5個，土壤指標7個，共計12個指標。

為了消除指標單位不同所帶來的誤差，使其數據具有可比性，采用極差法對研究數據進行標準化處理，公式為（朱乾隆等，2019）：

式中：

Yij——各指標的標準化值；

xij——各指標的原始數據值；

xmax和 xmin——第 i個指標在研究時段內的最大值和最小值。公式（1）和公式（2）分別表示正向指標與負向指標。

各個評價指標權重的確定在耦合模型評判中占有非常重要的位置，權重的大小對評估結果十分重要，它反映了各指標的相對重要性?？紤]到客觀賦權法過于強調各項評價指標數據之間的內部變化，缺乏對實際情況的針對性分析（唐李斌等，2020），所以在各指標進行標準化處理后，采用可將多維的非線性問題簡化成一維線性疊加的層次分析法（焦菊英等，2005）。最底層（要素層）是隸屬于上層的各個評價指標，每一層里2個元素進行相對于上層某元素的重要性相比，建立起每層的對比矩陣，并且通過一致性檢驗后（陳萍等，2011），計算得出各個因子權重值（表2）。

表2 土壤-植被系統耦合協調度指標及權重Table 2 The index and weight of coupling coordination degree of soil-vegetation system

1.4 模型建立

1.4.1 耦合度模型

耦合度是度量系統或要素之間相互作用和影響程度的定量指標，它刻畫了某一時點區域系統之間、系統內各要素之間交互脅迫、交互依存關系的演進態勢或趨向（王明全等，2009）。假設正數x1,x2, x3……xp為描述植被子系統特征的p個指標，設正數y1, y2, y3……yq為描述土壤子系統特征的q個指標，由此構建植被與土壤的耦合度模型為：

式中：

C——不同人工林分類型生態恢復植被土壤耦合度；

k——調節系數，一般為2≤k≤5，本文k=2。耦合度C是反映P(x)與S(y)之和在一定條件下，植被-土壤協調數量的程度，其值在0—1之間，當C趨近于0時，表明系統要素之間的耦合度極小，近乎處于無關狀態，當C趨近于1時，表明系統要素之間呈良性耦合狀態。P(x)為植物綜合評價函數，S(x)為土壤綜合評價指數。

式中：

i——植被特征指標；

j——土壤特征指標的個數；

ai——第i個植物指標；

bj——第j個土壤指標的權重值；

xi——第i個植物指標；

yj——第j個土壤指標的標準化值。由綜合函數計算得出的指數越高，意味著植被恢復或土壤環境越好，反之則越差。

1.4.2 耦合度協調模型

由于僅依靠耦合度很難全面反映出系統間或系統內各要素之間對植被-土壤系統的協調程度，甚至會對評價結果的準確性產生誤導。為了進一步反映不同人工林分類型生態恢復植物-土壤交互耦合的協調程度，在耦合度C的基礎上，構建了耦合協調度模型（徐明等，2016）：

式中：

D——不同人工林分類型生態恢復土壤-植被系統耦合協調度，D值在0—1之間，D值越趨近于1，則植被恢復與土壤環境的整體水平越高，反之則越低；

T——植被-土壤系統綜合協調指數；

α、β——待定系數，考慮到系統中植被恢復與土壤環境的改善同樣重要，所以α與β均為1/2，故：

1.5 植被-土壤系統協調發展類型的判斷標準

由于目前耦合協調度等級尚無統一的劃分標準，本文所采用的植被土壤耦合類型分類體系及評判標準是在參考前人研究結果（南國衛等，2021）的基礎上，結合本研究的相關指標計算結果得出的。根據D值的大小，同時結合植物綜合評價函數P(y)與土壤綜合評級函數 S(x)，得出耦合協調類型及評判標準（表3）。

表3 土壤-植被系統耦合類型劃分Table 3 Classification of coupling types of soil-vegetation systems

2 結果與分析

2.1 不同人工林分類型指標層的差異性

土壤指標層中，檸條與撂荒地的土壤全磷、土壤容重存在顯著性差異（P＜0.05）；側柏與撂荒地的土壤有機碳之間存在顯著性差異（P＜0.05）；山楊與撂荒地的土壤全氮存在顯著性差異（P＜0.05）；側柏、檸條、長梗扁桃、虎榛子與撂荒地的土壤孔隙度之間顯著性差異（P＜0.05）；土壤含水率及土壤全鉀在不同林分間無顯著差異（表4）。植物指標層中，虎榛子與撂荒地的Margalef指數存在顯著性差異（P＜0.05）；撂荒地與其他人工林分的Simpson指數存在顯著性差異（P＜0.05）；側柏、檸條、虎榛子、黃刺玫與撂荒地的Pielou指數存在顯著性差異（P＜0.05）；側柏、山楊與其他人工林分及撂荒地的生物量之間存在顯著性差異（P＜0.05），Shannonwiener指數在不同林分間無顯著差異（表5）。

表4 不同人工林分類型土壤指標層的差異性Table 4 The different of soil index layer in different artifical stand types

表5 不同人工林分類型植被指標層差異性Table 5 The different vegetation index layer among different artificial stand types

2.2 不同人工林分類型生態恢復土壤-植被系統的耦合關系評價

各不同人工林分生態恢復植被綜合指數P(y)、土壤綜合指數S(x)、植被土壤系統耦合度C以及耦合協調度D指標的計算結果如下（表6）?？傮w上，檸條具有較高的植被指數P(y)、耦合度指數C以及耦合協調指數D；虎榛子具有較高的土壤指數S(x)。植被指數具體表現為：檸條＞側柏＞虎榛子＞山楊＞長梗扁桃＞黃刺玫＞撂荒地，說明所有的人工林分的植物水平的生長水平均較高于撂荒地，特別是黃刺玫與側柏；土壤指數具體表現為：虎榛子＞長梗扁桃＞側柏＞山楊＞黃刺玫＞檸條＞撂荒地，說明所有的6種人工林分均能在不同程度上改良土壤質量，特別是黃刺玫與長梗扁桃。由于人工林栽植年限較短，生態恢復效果還未完全體現，以至于6種人工林分的耦合度與耦合協調度均不高，耦合度表現為：檸條＞側柏＞山楊＞虎榛子＞長梗扁桃＞黃刺玫＞撂荒地；耦合協調度表現為：檸條＞側柏＞黃刺玫＞山楊＞長梗扁桃＞虎榛子＞撂荒地。黃刺玫與撂荒地為嚴重失調衰退類植被損益型；側柏、虎榛子、長梗扁桃以及山楊為瀕臨失調衰退類植被損益型；檸條表現為瀕臨失調衰退類植被土壤共損型。由此來看，在九峰山人工林生態恢復過程中，檸條的耦合協調度更好，除撂荒地與黃刺玫外，人工林分的土壤-植被耦合協調系統發展趨勢基本處于瀕臨失調衰退型，這說明人工林分較撂荒地在生態恢復中起到一定的作用，但整體仍處于較低水平。

表6 不同人工林分生態恢復土壤-植被系統耦合協調關系評價Table 6 Evaluation of coupling coordination relationship between soil and vegetation systems in different artificial stands for ecological restoration

3 討論

本文以植被-土壤系統耦合理論為基礎，通過建立植被與土壤耦合模型，將植被與土壤兩個系統通過各自的調節權數耦合在一起，并進一步闡釋若干子系統之間的相互關系，完善評價過程（張艷等，2013；王淑佳等，2021），綜合客觀的評價九峰山生態恢復的效益。

植被-土壤是一個相互作用，相互影響的有機整體，植物的生長需要從土壤中吸取養分和水分，并需要土壤作為基礎支撐；地上部分植物生長及其覆蓋也在不斷改變著土壤理化性狀及土壤微生態環境（余軒等，2021）。通過耦合協調關系評價發現，不同人工林分的耦合度和耦合協調度并不一一對應，植被與土壤的綜合指數在不同人工林分類型中的位置也并不完全一致。如植被綜合指數表現為：檸條＞側柏＞虎榛子＞山楊＞長梗扁桃＞黃刺玫＞撂荒地，土壤綜合指數表現為：虎榛子＞長梗扁桃＞側柏＞山楊＞黃刺玫＞檸條＞撂荒地，這可能與植被和土壤的互作機制有關，植被生長消耗土壤養分和水分的速率與植被改良土壤理化性質的速率不一致所導致的。雖然通過長時間的人工造林，對生態恢復起到了一定作用，但不同林分的植被土壤耦合協調狀況并未達到理想狀態。在九峰山的長期造林過程中，不合理的栽植密度以及同齡純林育林模式，忽視了人工林的演替規律，在林地選擇時僅用水分、植被類型等因子作為依據，使得部分植被出現不同程度的衰退，進而導致九峰山不同人工林分類型生態恢復植被—土壤系統基本處于瀕臨失調衰退型發展水平，整體水平較低。在上述人工林分類型中，植被—土壤耦合協調度較高的是檸條，除檸條外，其余林分均表現為不同程度的失調衰退類土壤損益型。這可能與當地降雨無法滿足植物的正常生長有關，隨著植物根系不斷地生長，植物需要將根系伸進深層土壤來獲得水資源（王云強等，2012），深層土壤中水分的損失，進一步加劇土壤干旱，最終形成土壤干層（Yang et al.，2011；Mazzacavallo et al.，2017）。在土壤干層形成過程中，土壤理化性質也隨之下降，最終反作用于植物。但檸條具有耐旱和固氮等優良特性，并且灌叢根系深廣，可以吸收和消耗周圍較大面積和深層土壤中的水分（張建華等，2011），加上微環境的影響，使土壤蒸發減少，檸條灌叢土壤水分含量隨著冠幅增大而增大（劉佳楠等，2019），間接增加土壤養分和減少水分的虧缺，從而使檸條較其他林分存在一定優勢。在水分虧缺時，抗旱植被能以高效的水分利用方式來抵御干旱的脅迫（閆海龍等，2010），而旱生植物具有發達機械組織以及表皮有多層細胞等特征來減少水分散失，但黃刺玫相較于其他林分，它的海綿組織厚度和葉片組織結構疏松度最低（王樹森等，2020），因此在養分輸送功能和保水、貯水作用方面表現較差，導致較差的抗旱能力；其次，黃刺玫的根系基本上沒有垂直伸展的主根，但具有較發達的水平根和大量的斜生根，集中分布在0—40 cm土層中（王芳等，2006），這就導致土壤淺層水分被消耗到一定程度時，黃刺玫無法吸收土壤深層水分，并且其保水儲水能力較差，這也可能是導致黃刺玫評價為輕度失調衰退類土壤損益型的原因之一。

4 結論

（1）通過顯著性差異分析，植被指標層中除Shannon-wiener指數外，其余指標在不同林分間均存在顯著差異性；土壤指標層中除土壤含水率與土壤全鉀外，其余指標在不同林分間均存在顯著性差異。

（2）經過對不同人工林分類型的體系構建與評價，不同人工林分類型生態恢復植被-土壤耦合協調關系的差異較為小，整體上呈現不同程度的失調衰退類型。

（3）不同人工林分類型的植被恢復較差，土壤發展狀況優于植被發展狀況。

（4）檸條能夠改善植被-土壤系統中植被恢復與土壤環境的耦合協調性，提高生態恢復效益。