摘 要：為了探究祁連山青海云杉森林土壤質(zhì)量情況，運(yùn)用經(jīng)典統(tǒng)計(jì)方法，分析了祁連山青海云杉森林土壤質(zhì)量各項(xiàng)指標(biāo)的均值與變化范圍。通過主成分分析法，最小數(shù)據(jù)集篩選出5個(gè)指標(biāo)，計(jì)算得到基于最小數(shù)據(jù)集的土壤質(zhì)量指數(shù)平均值為0.583，評(píng)價(jià)結(jié)果表明祁連山青海云杉森林土壤質(zhì)量以低水平和較低水平為主。

關(guān)鍵詞：祁連山青海云杉森林；土壤質(zhì)量評(píng)價(jià)；主成分分析；最小數(shù)據(jù)集

中圖分類號(hào)：S158 文獻(xiàn)標(biāo)志碼：B 文章編號(hào)：2095–3305（2024）08–0-03

作為森林生態(tài)系統(tǒng)的主要組成成分，森林土壤面積在全球范圍內(nèi)占有較大比例，為森林生態(tài)系統(tǒng)提供了生存所需的水分、養(yǎng)分、空氣、熱量等物質(zhì)，對(duì)森林的健康發(fā)展具有重要作用[1]。近年來，隨著森林生態(tài)系統(tǒng)在森林資源經(jīng)營(yíng)管理中的地位日益突出，森林土壤在森林生態(tài)系統(tǒng)中的作用也越來越受到關(guān)注[2]。因此，在今后的研究中，應(yīng)深入探討不同林型的土壤質(zhì)量，為林業(yè)的健康管理和可持續(xù)發(fā)展提供理論基礎(chǔ)。

當(dāng)前，國(guó)內(nèi)外已經(jīng)研究出多種評(píng)估森林土壤質(zhì)量的方法，主要有層次分析法、最小數(shù)據(jù)集法（MDS）、GIS結(jié)合數(shù)學(xué)模型技術(shù)和灰色關(guān)聯(lián)法。劉湘君等[3]研究了基于土壤理化性質(zhì)和剖面特征的黃淮海干旱農(nóng)業(yè)區(qū)耕地土壤質(zhì)量特征，利用主成分分析法，對(duì)黃淮海干旱農(nóng)業(yè)區(qū)耕地土壤質(zhì)量進(jìn)行了評(píng)價(jià)。當(dāng)前應(yīng)用MDS方法所做的品質(zhì)評(píng)估多集中在南方林地，而在我國(guó)北方鮮有相關(guān)研究。

1 材料和方法

1.1 試驗(yàn)區(qū)概況

位于甘肅省河西走廊中段的馬蹄林區(qū)是研究區(qū)馬蹄自然保護(hù)站的一部分，坐落于祁連山國(guó)家級(jí)自然保護(hù)區(qū)的試驗(yàn)區(qū)，年均氣溫1～3 ℃，日照時(shí)間約120 d，年平均降水量在360～490 mm之間[4]。從土壤類型來看，林區(qū)主要為山地灰褐土，優(yōu)勢(shì)樹種為青海云杉和祁連圓柏，植被類型和土壤類型具有明顯的垂直分布特征[5]。

1.2 研究方法

1.2.1 主成分分析法

主成分分析方法（PCA）是確定主要數(shù)據(jù)集的一種常用方法，也是評(píng)估土壤質(zhì)量的重要手段。全面考慮土壤質(zhì)量參數(shù)對(duì)土壤質(zhì)量的影響[6]，利用PCA方法對(duì)最小數(shù)據(jù)集進(jìn)行分析，以篩選出對(duì)土壤特性變化敏感的參數(shù)。針對(duì)同一主成分的選定測(cè)量指標(biāo)進(jìn)行相關(guān)性分析，如果指標(biāo)之間存在顯著的正相關(guān)性，即在保留主要特征的同時(shí)，可以消除冗余信息，并且相關(guān)系數(shù)＞0.5，則將正常值作為參考標(biāo)準(zhǔn)，將其納入最小數(shù)據(jù)集[7]。

Nik=（1）

1.2.2 土壤質(zhì)量指數(shù)計(jì)算

土壤品質(zhì)綜合評(píng)價(jià)的計(jì)算包括3個(gè)步驟：首先，對(duì)指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理；其次，利用主成分分析（PCA）確定指標(biāo)的權(quán)重；最后，通過加權(quán)求和得出土壤品質(zhì)指數(shù)（（SQI））[8-9]。

升型函數(shù)（2）和降型函數(shù)（3）公式為：

（2）

（3）

評(píng)價(jià)指標(biāo)的得分f（x）與各指標(biāo)的測(cè)量值x相關(guān)，其中，最大值為U，最小值為L(zhǎng)。

Wi=Ci/Ci（4）

式（4）中，Wi是各指標(biāo)的權(quán)重；Ci是各評(píng)價(jià)指標(biāo)的公因子方差；n是MDS中包含的指標(biāo)數(shù)。

SQI=Wi×fi（5）

式（5）中，SQI為土壤質(zhì)量指數(shù)，fi為各評(píng)價(jià)指標(biāo)的得分值，n為各個(gè)數(shù)據(jù)集中的土壤質(zhì)量指標(biāo)數(shù)量。

1.2.3 數(shù)據(jù)處理

利用Excel軟件進(jìn)行數(shù)據(jù)和統(tǒng)計(jì)分析，借助統(tǒng)計(jì)軟件SPSS 26軟件進(jìn)行主成分分析、方差分析和LSD數(shù)據(jù)的多重比較，使用Origin 2022軟件對(duì)處理后的數(shù)據(jù)繪制圖表。

2 土壤質(zhì)量評(píng)價(jià)

2.1 土壤質(zhì)量指標(biāo)統(tǒng)計(jì)分析

在參考袁程昱[10]等對(duì)土壤質(zhì)量指標(biāo)選擇的基礎(chǔ)上，選取土壤有機(jī)質(zhì)（SOM）、孔隙度（POR）、土壤容重（BD）、全磷（TP）、全鉀（TK）、全氮（TN）、水解氮（AN）、速效鉀（AK）、速效磷（AP）和pH值10個(gè)因子作為土壤質(zhì)量評(píng)價(jià)指標(biāo)（表1）。

由表1結(jié)果可知，祁連山青海云杉森林土壤質(zhì)量指標(biāo)空間變異性具有明顯差距，其變化范圍在1.18%～65.31%之間，變異系數(shù)從大到小依次為全鉀、水解氮、全磷、土壤有機(jī)質(zhì)、速效鉀、全氮、速效磷、土壤容重、孔隙度、土壤pH值。根據(jù)變異系數(shù)分級(jí)[11]，土壤容重、土壤pH值和孔隙度是不敏感指標(biāo)；土壤有機(jī)質(zhì)、全氮、全磷、水解氮、速效磷、速效鉀為低敏感指標(biāo)；全鉀是中敏感指標(biāo)。根據(jù)第二次全國(guó)土壤普查的養(yǎng)分評(píng)價(jià)標(biāo)準(zhǔn)[12]，祁連山青海云杉森林土壤氮、磷、鉀含量很高，平均全氮含量為5.92 g/kg，這說明土壤基礎(chǔ)肥力較高。

2.2 最小數(shù)據(jù)集指標(biāo)的篩選

經(jīng)過對(duì)10個(gè)土壤指標(biāo)的主成分分析，發(fā)現(xiàn)各土壤質(zhì)量指標(biāo)的因子載荷表現(xiàn)出4個(gè)主成分，其特征值分別為2.825、1.969、1.314、1.124，均大于1。這些主成分的累計(jì)貢獻(xiàn)率達(dá)到72.329%，表明這些指標(biāo)能夠反映土壤質(zhì)量70%左右的信息。

根據(jù)表2，從主成分中挑選出旋轉(zhuǎn)因子載荷絕對(duì)值≥0.5的土壤質(zhì)量評(píng)價(jià)指標(biāo)。在主成分1中，發(fā)現(xiàn)土壤有機(jī)質(zhì)、全氮和速效鉀這3個(gè)指標(biāo)的載荷值達(dá)到了0.5及以上，歸為第一組指標(biāo)；在主成分2中，土壤容重和孔隙度的因子載荷值比主成分1中的指標(biāo)更高，劃分為第二組，主成分2中的指標(biāo)包括土壤容重、全鉀和孔隙度；主成分3中的指標(biāo)則包括全磷、水解氮和速效磷，為第三組指標(biāo)；主成分4中僅有土壤pH值這一個(gè)指標(biāo)，歸為第四組。

為進(jìn)一步減少信息冗余，通過計(jì)算得出各土壤指標(biāo)的正常值。結(jié)果表明：在第一組中，土壤有機(jī)質(zhì)Norm值最大，為2.116，其他指標(biāo)不在土壤有機(jī)質(zhì)Norm值的10%范圍內(nèi)因?yàn)楸惶蕹?；第二組中Norm值最大的是土壤容重，為1.304，孔隙度在其值的10%范圍內(nèi)，且與其顯著相關(guān)（P＞0.01），第二組土壤容重與孔隙度被選入最小數(shù)據(jù)集；第三組全磷Norm值最大，為0.826，水解氮與速效磷不在其Norm值10%的范圍內(nèi)，被剔除；第四組僅包含土壤pH值，入選最小數(shù)據(jù)集。綜上，祁連山青海云杉森林土壤質(zhì)量指標(biāo)包括土壤有機(jī)質(zhì)、土壤容重、孔隙度、全磷和土壤pH值。

2.3 土壤質(zhì)量評(píng)價(jià)

采用主成分分析方法得到的共同因素方差之和，求出各個(gè)數(shù)據(jù)集合中各指標(biāo)的權(quán)重。利用公式（2）和（3）將各個(gè)指標(biāo)在主分量上的因素負(fù)荷的正與負(fù)來計(jì)算隸屬度，并將所得的隸屬度與加權(quán)值乘以，再由公式（4）算出以整個(gè)土壤質(zhì)量指數(shù)為基礎(chǔ)的土壤質(zhì)量指數(shù)SQI-TDS。

祁連山青海云杉森林基于最小數(shù)據(jù)集的土壤質(zhì)量指數(shù)平均值為0.583，變化范圍在0.136～0.758之間，變異系數(shù)24.361%，變異系數(shù)較大，屬于中等變異（表4）。

根據(jù)土壤質(zhì)量指數(shù)分級(jí)表和上述2個(gè)土壤質(zhì)量指數(shù)的平均值和范圍，可以看出祁連山青海云杉森林土壤質(zhì)量主要處于低水平和較低水平。

3 結(jié)論

本研究采用傳統(tǒng)統(tǒng)計(jì)方法，研究了祁連山馬蹄自然保護(hù)站范圍內(nèi)青海云杉森林土壤質(zhì)量的各項(xiàng)參數(shù)變化。通過應(yīng)用最小數(shù)據(jù)集（MDS）評(píng)估土壤質(zhì)量，得出以下結(jié)論：

（1）祁連山青海云杉森林土壤質(zhì)量基于最小數(shù)據(jù)集的土壤質(zhì)量指數(shù)平均值為0.583，變化范圍在0.136～0.758之間。評(píng)價(jià)結(jié)果表明，祁連山青海云杉森林土壤質(zhì)量以低水平和較低水平為主。

（2）利用主成分析法篩選出的最小數(shù)據(jù)集評(píng)估土壤品質(zhì)是一個(gè)科學(xué)、實(shí)用的方法。研究成果為該地區(qū)土地質(zhì)量評(píng)估提供了有價(jià)值的信息，在此基礎(chǔ)上優(yōu)化和篩選出精確的土壤指標(biāo)，從而構(gòu)建全面的土壤質(zhì)量量化評(píng)估指標(biāo)。此項(xiàng)研究將為馬蹄自然保護(hù)站下屬地區(qū)的青海云杉森林在土壤和土地管理及規(guī)劃方面提供有力的數(shù)據(jù)基礎(chǔ)。

參考文獻(xiàn)

[1] 張繼平，張林波，王風(fēng)玉，等.井岡山國(guó)家級(jí)自然保護(hù)區(qū)森林土壤養(yǎng)分含量的空間變化[J].土壤，2014，46（2）：262-268.

[2] 方偉東.長(zhǎng)白山地區(qū)四種森林類型土壤理化性質(zhì)及水源涵養(yǎng)功能[D].北京：北京林業(yè)大學(xué)，2011.

[3] 劉湘君，喬冠宇，郭豐浩，等.基于最小數(shù)據(jù)集的黃淮海旱作區(qū)耕層土壤質(zhì)量評(píng)價(jià)及障礙分析[J].農(nóng)業(yè)工程學(xué)報(bào)， 2023，39（12）：104-113.

[4] 趙維俊，楊永華，白登忠.祁連山排露溝流域不同海拔梯度青海云杉林生長(zhǎng)特征[J].林業(yè)科技通訊，2017（12）：15-18.

[5] 劉濟(jì)萌，劉賢德，馬瑞，等.祁連山排露溝流域青海云杉天然更新空間關(guān)聯(lián)性研究[J].干旱區(qū)資源與環(huán)境，2024，38（1）： 161-170.

[6] 解雪峰，濮勵(lì)杰，朱明，等.基于MDS與TOPSIS模型的濱海灘涂圍墾區(qū)土壤質(zhì)量評(píng)價(jià)[J].環(huán)境科學(xué)，2019，40（12）： 5484-5492.

[7] 趙賀，王緒奎，劉紹貴，等.基于水稻產(chǎn)量的江蘇省稻麥輪作區(qū)土壤質(zhì)量評(píng)價(jià)[J].土壤，2020，52（6）：1230-1238.

[8] 王安寧，劉歌暢，徐學(xué)華，等.鐵尾礦廢棄地不同復(fù)墾年限土壤質(zhì)量評(píng)價(jià)[J].北京林業(yè)大學(xué)學(xué)報(bào)，2020，42（1）：104-113.

[9] 李鵬飛，張興昌，郝明德，等.基于最小數(shù)據(jù)集的黃土高原礦區(qū)復(fù)墾土壤質(zhì)量評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào)，2019，35（16）：265 -273.

[10] 袁程昱，李發(fā)永，胡雪菲，等.基于最小數(shù)據(jù)集的高寒高海拔礦區(qū)土壤質(zhì)量評(píng)價(jià)[J].湖北理工學(xué)院學(xué)報(bào)，2021，37 （5）：17-23.

[11] 史東梅，江娜，蔣光毅，等.紫色土坡耕地耕層質(zhì)量影響因素及其敏感性分析[J].農(nóng)業(yè)工程學(xué)報(bào)，2020，36（3）：135-143.

[12] 宋洪炳，郭偉玲.基于最小數(shù)據(jù)集的江淮丘陵區(qū)土壤質(zhì)量綜合評(píng)價(jià)[J].農(nóng)業(yè)與技術(shù)，2023，43（12）：93-98.