朱家晸，秦富倉*，李 龍，楊振奇，芳 菲，趙 琦

（1．內(nèi)蒙古農(nóng)業(yè)大學(xué)沙漠治理學(xué)院，內(nèi)蒙古 呼和浩特 010018；2．荒漠生態(tài)系統(tǒng)保護與修復(fù)國家林業(yè)和草原局重點實驗室，內(nèi)蒙古 呼和浩特 010000；3．阿拉善左旗科學(xué)技術(shù)和林業(yè)草原局，內(nèi)蒙古 阿拉善盟 750300）

生態(tài)人工林主要為滿足保護和改善人類生存環(huán)境、維護生態(tài)平衡、保存物種資源、科學(xué)試驗等需要，其主要作用是增強森林涵養(yǎng)水源能力、調(diào)節(jié)當(dāng)?shù)匦夂蛞约氨３稚鷳B(tài)平衡等［1］。根據(jù)第九次全國森林資源清查內(nèi)蒙古自治區(qū)清查結(jié)果，內(nèi)蒙古現(xiàn)有森林面積2614.85 萬hm2，其中人工林面積600.01萬hm2，占22.94%，面積巨大。這其中，土壤作為森林生態(tài)系統(tǒng)的重要組成部分，對森林生態(tài)系統(tǒng)的養(yǎng)分循環(huán)及健康發(fā)育均有著重要影響［2］。

林下土壤作為供應(yīng)樹木成長所需營養(yǎng)的重要基質(zhì)，對林木的生長和分布狀況有著重要作用［3］。林下土壤肥力是森林土壤性能的重要指標。國內(nèi)外研究表明，土壤肥力主要受土壤的物質(zhì)組成、結(jié)構(gòu)、環(huán)境等因素的影響，是土壤理化性質(zhì)的綜合反映［4］。土壤肥力狀況能夠反映植物從土壤中吸取養(yǎng)分的潛力，主要包括土壤的機械組成、有機質(zhì)含量、氮、磷和鉀的全量和有效量，這些指標對樹木的生長有著重要影響。不同林分由于樹種、造林模式、郁閉度等組成條件的不同，會導(dǎo)致土壤理化性質(zhì)發(fā)生變化，進而影響森林生態(tài)系統(tǒng)結(jié)構(gòu)變化，土壤抵抗鹽堿化的能力也隨之改變，并可能最終改變森林演替軌跡［5-6］。因此，深入研究不同林分類型的土壤肥力狀況，對評價土壤質(zhì)量狀況、科學(xué)利用土地資源、合理進行森林經(jīng)營管理以及維護森林生態(tài)系統(tǒng)平衡均具有重要意義。

目前常用的傳統(tǒng)土壤肥力評價方法主要分為經(jīng)驗權(quán)重法和統(tǒng)計權(quán)重法［7-10］。不同評價方法采用的體系和側(cè)重點不同［11-12］，其最終得出的結(jié)果也就有所不同。本文主要采用更側(cè)重于數(shù)據(jù)本身體現(xiàn)的統(tǒng)計學(xué)意義和數(shù)理統(tǒng)計方法的統(tǒng)計權(quán)重中的主成分分析法和相關(guān)性分析法。

為研究不同林分類型的土壤肥力狀況，本研究以內(nèi)蒙古自治區(qū)清水河縣公益林區(qū)內(nèi)的油松林、落葉松林、山杏林、喬灌混交林（油松+檸條）、針闊混交林（銀白楊+落葉松）5 種主要林分人工林林下土壤為研究對象，采用相關(guān)性分析法［10］以及主成分分析法［13］對土壤肥力進行綜合評價。旨在探討對比不同林分類型下人工林土壤肥力的狀況。為黃土丘陵區(qū)生態(tài)公益林管理和對森林土壤資源的科學(xué)管理提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

清水河縣位于內(nèi)蒙古自治區(qū)呼和浩特市最南端（111°21′～112°07′E，39°35′～40°11′N），海拔960 ～1837 m?？h境內(nèi)平川95 km2，山地974 km2，丘陵1499 km2，構(gòu)成了溝壑和平原并存的地貌類型。林木保存面積8.8 萬hm2，森林覆蓋率30.8%，主要樹種有油松、落葉松、檸條、沙棘、杏等。清水河縣地處中溫帶，屬典型的溫帶大陸性季風(fēng)氣候；雨熱同季，雨量集中；主要集中在6 ～9 月，占年降水量的80%；歷年平均降水量413.8 mm，平均降水日數(shù)為75 d，全縣年蒸發(fā)量達2577.2 mm；年平均氣溫7.1℃。主要有栗鈣土、栗褐土、灰褐土、潮土、風(fēng)沙土、沼澤土、鹽土、石質(zhì)土8 個土類；地帶性土壤為栗鈣土。

1.2 土壤樣品采集

在查閱文獻和實地調(diào)查的基礎(chǔ)上，于2020年7 ～9 月在清水河縣公益林區(qū)內(nèi)，根據(jù)樹種分布現(xiàn)狀和數(shù)量，篩選出立地條件相同、具有代表性的5 種典型林分類型，包括油松林、落葉松林、山杏林、喬灌混交林（油松+檸條）、針闊混交林（銀白楊+落葉松），林齡均在15 年以上，在其中選取具有代表性的區(qū)域，并在每種林分代表區(qū)域中隨機設(shè)置5 個樣地斑塊，共計25塊樣地，每塊樣地面積均為20 m×20 m，在樣地對角線上設(shè)置3 個采樣點，每個采樣點面積為1 m×1 m；收集0 ～20 cm 土層的土壤，去除表層腐殖質(zhì)、殘根等，平均分成4 份共收集1 kg 的樣品，編號后帶回實驗室，風(fēng)干后過篩備用。樣地基本情況見表1。

1.3 土壤化學(xué)性質(zhì)測定

土壤pH 值采用混合指示劑比色法測定；有機質(zhì)含量采用重鉻酸鉀-硫酸溶液氧化法測定；有效磷含量用NaHCO3浸提-鉬銻抗比色法測定；全磷、全氮含量采用HClO4-H2SO4消化法測定；速效鉀含量采用中性乙酸銨浸提法測定；堿解氮含量用堿解擴散法測定［14］。

表1 樣地基本情況

1.4 評價方法

1.4.1 指標標準化處理

在進行評價之前，因為各項待評價指標的量綱數(shù)據(jù)具有差異性，因而需要解決其統(tǒng)一標準化的問題，進而完成各項評價指標的數(shù)據(jù)統(tǒng)一化，使各評價指標間具備可比性［4］。因此，需要創(chuàng)建出各項評價指標的隸屬度函數(shù)，并計算各指標的隸屬度函數(shù)值。按照相關(guān)文獻記載以及不同樹種的實際生長情況，可將隸屬度函數(shù)分為拋物線型和S 型。

拋物線型隸屬度函數(shù)：評價指標體系中具備一定范圍的樹木最佳生長區(qū)間。假設(shè)其指標數(shù)據(jù)離開此區(qū)間，則認定不適于樹木的生長。指標數(shù)據(jù)的偏離程度可以反映其對樹木生長的影響，pH 指標離散程度較低，適用于此函數(shù)。隸屬度函數(shù)值F（x）的計算公式為：式中，x 為測定評價指標的實際值，xa、xb、xc、xd為函數(shù)曲線轉(zhuǎn)折點的取值。

S 型隸屬度函數(shù)：S 型隸屬度函數(shù)的曲線特征是，指標假定在一固定數(shù)值區(qū)間內(nèi)增長時，可以反映出評價指標對樹種的生長具有較大影響，一旦低于或高于此區(qū)間，則評價指標影響較小，適用于數(shù)據(jù)離散程度較大的土壤有機質(zhì)、全氮、全磷、有效磷、速效鉀、堿解氮含量。式中，x 為測定的評價指標的實際值，xa、xb為函數(shù)曲線轉(zhuǎn)折點取值。

綜上，根據(jù)公式（1）和（2）可將指標數(shù)據(jù)轉(zhuǎn)化為0.1 ～1.0 的無綱量化數(shù)值，數(shù)值越接近1.0 表明指標對樹木成長的影響越大。依據(jù)全國第二次土壤普查數(shù)據(jù)［15］和查閱相關(guān)的文獻［16-17］。函數(shù)曲線轉(zhuǎn)折點的取值如表2 所示。

表2 隸屬度函數(shù)的轉(zhuǎn)折點取值

1.4.2 綜合評價方法

采用下列2 種評價方法對5 種林分類型進行土壤肥力綜合評價。

（1）主成分分析法：運用SPSS 22.0 對評價指標實際值進行降維處理，重新組合成一組獨立的綜合指標，采用加權(quán)求和法計算公因子方差，進一步計算各個公因子方差占公因子方差總和的比例［18］，將其作為指標的權(quán)重αi。土壤肥力綜合得分C 計算公式為：

式中，αi為第i 個土壤肥力指標占總體的權(quán)重，n為指標數(shù)量，F(xiàn)（i）為某一林分類型各項指標的隸屬度值。

（2）相關(guān)性分析法：主要研究兩個或兩個以上處于同等地位的隨機變量間的相關(guān)關(guān)系的統(tǒng)計分析方法。在存在多變量的情況下，由于變量之間存在著直接或間接的影響，所以需要消除其他變量的影響。在土壤肥力評價體系中，在控制其他指標的影響后，可根據(jù)各指標間的相關(guān)程度，建立偏相關(guān)系數(shù)矩陣，計算其平均值，求出該平均值占所有指標平均值之和的百分比。最終求得單項肥力指標在表征土壤肥力狀況中的貢獻率，即權(quán)重系數(shù)Wi［19］，最終基于Wi和各指標的隸屬度函數(shù)值可求得綜合指數(shù)IFI［20］：

式中，Wi為第i 個土壤肥力指標的權(quán)重系數(shù)，n 為指標數(shù)量，F(xiàn)（i）為某一林分類型各項指標的隸屬度值。

1.5 數(shù)據(jù)分析

數(shù)據(jù)處理采用Excel 2010，使用拋物線型和S型曲線，對實際觀測數(shù)據(jù)進行統(tǒng)一化處理，使用SPSS 22.0 進行數(shù)據(jù)分析，對評價指標進行描述性統(tǒng)計及相關(guān)性分析，利用相關(guān)性分析法以及主成分分析法分析計算得出不同林分類型林下土壤肥力分值。

2 結(jié)果與分析

2.1 不同林分類型土壤肥力指標分析

根據(jù)表3 顯示的各項土壤肥力指標的描述性統(tǒng)計結(jié)果顯示，數(shù)據(jù)均呈正態(tài)分布，變異系數(shù)由大到小為速效鉀＞全氮＞堿解氮＞有效磷＞全磷＞有機質(zhì)＞pH。其中pH、有機質(zhì)為低于10%的弱變異，數(shù)據(jù)離散程度較小，說明5 種林分土壤pH、有機質(zhì)兩項指標數(shù)據(jù)較為接近；其他6 項指標均屬中等變異，速效鉀為最高，說明5 種林分的土壤速效鉀含量差異較大。

表3 土壤指標描述性統(tǒng)計

不同林分類型的土壤肥力指標測定值如表4 所示，土壤pH 平均值為8.28，表明研究區(qū)土壤為堿性土壤。根據(jù)全國第二次土壤普查分類標準［16］（表5），5 種林分類型中，5 種林分的土壤全磷、有效磷含量均為三級，屬于中上水平；堿解氮含量平均為68.19 ～88.66 mg·kg-1，均為四級中下水平，整體含量偏低；落葉松林、山杏林的有機質(zhì)、全氮含量為四級，油松林、喬灌混交林、針闊混交林為三級；針闊混交林的速效鉀含量為二級較高水平，油松林、山杏林、喬灌混交林為三級，落葉松林為四級。

表4 不同林分類型土壤肥力指標測定值

表5 全國第二次土壤普查養(yǎng)分分級標準

由表4 可得，除pH 外，以落葉松+楊樹混交為主要樹種的針闊混交林的各項肥力指標數(shù)值均高于其他4 種林分，土壤肥力狀況較其他林分更好，5 種林分中，山杏林肥力狀況為最差。且數(shù)據(jù)顯示，混交林的土壤肥力要優(yōu)于純林。5 種林分的土壤全磷、有效磷含量整體較豐富，其中針闊混交林與山杏林的有效磷含量差異顯著，差值達到4.69 mg·kg-1。速效鉀含量整體差異較大，平均值為156.16 ～97.65 mg·kg-1，其中油松林、山杏林和喬灌混交林之間無顯著差異，落葉松林與針闊混交林差異顯著。5 種林分有機質(zhì)、全氮和堿解氮含量平均值分 別 為18.06 ～19.84 g·kg-1、0.82 ～1.00 g·kg-1和68.19 ～88.66 mg·kg-1，整體無顯著差異。

2.2 土壤肥力單項指標評價

根據(jù)雷達圖的幾何意義，坐標點離原點越遠表明指標的狀態(tài)越好，各坐標點圍成的面積越大，則表明其評價對象的整體狀態(tài)更佳。因此，雷達圖可更加直觀地反映出單項土壤肥力指標在土壤中的狀態(tài)及土壤肥力的整體情況［3］，5 種林分各項土壤肥力指標隸屬度函數(shù)值如圖1 所示。5 種林分類型各項土壤肥力指標的隸屬度值最小的均為pH，隸屬度值基本在0.5 左右，說明pH 對林下土壤質(zhì)量的作用分值較小。隸屬度值最高為全氮，隸屬度值基本在0.7 ～0.8 之間，作用分值較大。其他指標隸屬度值基本在0.55 ～0.7 之間，處于中間狀態(tài)。其中，最高的pH 隸屬度值出現(xiàn)在落葉松林，最高的有機質(zhì)、全磷、全氮、有效磷、速效鉀和堿解氮隸屬度值均出現(xiàn)在針闊混交林，說明針闊混交林的整體土壤肥力狀況較其他4 種林分更好。

圖1 中各項指標隸屬度值圍合而成的面積能夠直觀反映出各項指標對土壤肥力的作用影響，總體上，面積越大，林下土壤肥力情況越好。但是，以上數(shù)據(jù)僅是依據(jù)土壤肥力指標反映了土壤肥力的狀態(tài)且是建立在土壤肥力相等的前提下，實際上，由于指標間對土壤肥力的影響具有差異性，有必要結(jié)合單一土壤肥力指標的綜合平均貢獻率進一步驗證。

圖1 各項土壤肥力指標隸屬度函數(shù)值雷達圖

對不同林分類型的土壤肥力指標隸屬度值進行相關(guān)分析，并通過相關(guān)系數(shù)逆矩陣求得各指標的權(quán)重系數(shù)，進而可知各項土壤肥力指標對土壤肥力的貢獻率，結(jié)果如圖 2 所示。不同林分類型各項指標綜合平均貢獻率從大到小為：全氮（0.199）＞全磷（0.163）＞有效磷（0.108）＞堿解氮（0.099）＞有機質(zhì)（0.090）＞速效鉀（0.088）＞pH（0.078）。由此可以看出，全氮與全磷對改善土壤肥力有著重要作用。

圖2 各項土壤肥力指標隸屬度函數(shù)值綜合平均貢獻率雷達圖

2.3 土壤肥力綜合評價

2.3.1 主成分分析法

運用SPSS 22.0 對7 項土壤肥力指標進行檢驗，KMO 測 驗 結(jié) 果 為0.611（＞0.6），Bartlett 球 形檢驗的顯著系數(shù)為0.00（＜0.05），說明數(shù)據(jù)適宜進行主成分分析。經(jīng)過主成分分析得到各主成分數(shù)據(jù)值（表6），共提取出3 個主成分，各主成分的特征值分別為3.152、1.360 和1.054，前3 個主成分的累計方差貢獻率分別達39.562%、63.740%和79.510%。第1 主成分包括有機質(zhì)、全磷、有效磷和堿解氮，第2 主成分包括pH 和速效鉀，第3 主成分為全氮。基于公因子方差求出的權(quán)重系數(shù)表現(xiàn)為全氮（0.161）＞有效磷（0.159）＞全磷（0.153）＞速效鉀（0.139）＞堿解氮（0.135）＞有機質(zhì)（0.127）＞pH（0.123），全氮對土壤肥力的貢獻率最大，其次為有效磷和全磷。

按照各林分類型的隸屬度函數(shù)值（圖1）以及求出各指標的權(quán)重（表6），通過對其加權(quán)計算，可得到的土壤肥力綜合評價值如表7 所示，綜合得分為0.586 ～0.737，大小順序為針闊混交林＞喬灌混交林＞油松林＞落葉松林＞山杏林，針闊混交林分值最高，山杏林最小。

2.3.2 相關(guān)性分析法

基于各項指標的隸屬度函數(shù)值（圖1）和各指標對土壤肥力的貢獻率（圖2）進行加權(quán)計算，對不同林分類型各項土壤肥力指標進行綜合評價，結(jié)果如表7 所示。綜合指數(shù)為0.495 ～0.629，與主成分分析法結(jié)果一致。大小順序為針闊混交林＞喬灌混交林＞油松林＞落葉松林＞山杏林，針闊混交林肥力狀況最好，山杏林最差。

綜上，兩種評價方法進行土壤肥力綜合評價的最終結(jié)果一致，均為針闊混交林＞喬灌混交林＞油松林＞落葉松林＞山杏林，說明數(shù)據(jù)本身對兩種方法體現(xiàn)出統(tǒng)一性，更具有說服力，且兩種方法適用于土壤肥力的綜合評價。二者最終得出的評價分值實際較為接近，但指標權(quán)重有一定差別，例如相關(guān)性分析法得出的全磷平均綜合貢獻率在所有指標中排名第2，為0.163，但主成分分析法得出的權(quán)重為0.153，在指標中位列第3，說明不同評價方法偏重的內(nèi)容有所不同。主成分分析法主要是通過降維將原來具有一定相關(guān)性的變量重新組合成一組新的相互無關(guān)的幾個綜合變量，在盡可能多地保留原始變量信息的情況下，從中導(dǎo)出少數(shù)幾個主成分［21］；而相關(guān)性分析法更多的是考察多個變量間的相關(guān)性，通過指標之間的敏感度大小來確定指標權(quán)重［3］。所以，在方法的選擇上，還應(yīng)根據(jù)指標體系的具體情況來點定。

表6 土壤肥力指標的主成分分析

表7 不同林分類型的土壤肥力綜合評價

3 結(jié)論與討論

本研究結(jié)果表明，研究區(qū)內(nèi)蒙古清水河縣生態(tài)公益林區(qū)土壤pH 平均值為8.35，屬于堿性土壤，這與以往關(guān)于黃土丘陵區(qū)土壤pH 值研究結(jié)果一致［22］。其中，混交林的土壤肥力狀況要高于純林。混交林中針闊混交林（喬木混交）土壤肥力優(yōu)于喬灌混交林，純林中油松林＞落葉松林＞山杏林。

混交林土壤肥力表現(xiàn)較純林更好，基本與前人研究結(jié)果一致［23-24］，分析原因是由于人為對林分結(jié)構(gòu)進行改造，改變了樹木的空間分布格局［25］，使其能夠更加合理地吸收土壤養(yǎng)分，使土壤通透性增加，土壤結(jié)構(gòu)發(fā)生改變，更有利于養(yǎng)分的積累［26］；且混交林林內(nèi)枯落物較純林更為豐富，成分復(fù)雜，為林下土壤微生物的活動提供有利條件［27］，使得土壤微生物更為活躍，加速了枯落物的分解和養(yǎng)分釋放的速度，改良了土壤理化性質(zhì)，進而提高了土壤肥力。5 種林分中，山杏林土壤肥力表現(xiàn)最差，原因主要是山杏成長速度比其余樹種快，在生長過程中對營養(yǎng)的需求量較大，且林下植被稀少，最終導(dǎo)致枯落物蓄積量較少，只是表層養(yǎng)分積累，土壤整體肥力狀況不佳。但也有研究表明，山杏林土壤肥力狀況較油松等針葉林更好［28］，推測可能與當(dāng)?shù)負嵊芾泶胧┯嘘P(guān)。綜上所述，在植被恢復(fù)和撫育管理人工林的過程中，在遵從自然規(guī)律以及因地制宜的原則上，可以適當(dāng)增加混交林所占比例，引入適宜的混交樹種，并根據(jù)土壤肥力狀況進行科學(xué)的土壤養(yǎng)分管理。

由主成分分析法計算出的單項指標權(quán)重及雷達圖所呈現(xiàn)的貢獻率均說明全氮在當(dāng)?shù)赝寥婪柿χ胸暙I最大，其含量的大小與土壤養(yǎng)分環(huán)境密切相關(guān)［29］，可在今后作為指示土壤肥力的第一因素。pH 是土壤化學(xué)性質(zhì)的綜合體現(xiàn)［30-31］，對土壤養(yǎng)分的存在形態(tài)有直接影響［32］。本研究表明，pH 是土壤肥力的第一限制性因子，因此應(yīng)注意調(diào)節(jié)土壤pH，可在土壤中適當(dāng)添加以石膏或磷石膏為主的土壤改良劑用以降低土壤堿性［33］。