朱家晸,秦富倉(cāng)*,李 龍,楊振奇,芳 菲,趙 琦
(1.內(nèi)蒙古農(nóng)業(yè)大學(xué)沙漠治理學(xué)院,內(nèi)蒙古 呼和浩特 010018;2.荒漠生態(tài)系統(tǒng)保護(hù)與修復(fù)國(guó)家林業(yè)和草原局重點(diǎn)實(shí)驗(yàn)室,內(nèi)蒙古 呼和浩特 010000;3.阿拉善左旗科學(xué)技術(shù)和林業(yè)草原局,內(nèi)蒙古 阿拉善盟 750300)
生態(tài)人工林主要為滿足保護(hù)和改善人類生存環(huán)境、維護(hù)生態(tài)平衡、保存物種資源、科學(xué)試驗(yàn)等需要,其主要作用是增強(qiáng)森林涵養(yǎng)水源能力、調(diào)節(jié)當(dāng)?shù)匦夂蛞约氨3稚鷳B(tài)平衡等[1]。根據(jù)第九次全國(guó)森林資源清查內(nèi)蒙古自治區(qū)清查結(jié)果,內(nèi)蒙古現(xiàn)有森林面積2614.85 萬hm2,其中人工林面積600.01萬hm2,占22.94%,面積巨大。這其中,土壤作為森林生態(tài)系統(tǒng)的重要組成部分,對(duì)森林生態(tài)系統(tǒng)的養(yǎng)分循環(huán)及健康發(fā)育均有著重要影響[2]。
林下土壤作為供應(yīng)樹木成長(zhǎng)所需營(yíng)養(yǎng)的重要基質(zhì),對(duì)林木的生長(zhǎng)和分布狀況有著重要作用[3]。林下土壤肥力是森林土壤性能的重要指標(biāo)。國(guó)內(nèi)外研究表明,土壤肥力主要受土壤的物質(zhì)組成、結(jié)構(gòu)、環(huán)境等因素的影響,是土壤理化性質(zhì)的綜合反映[4]。土壤肥力狀況能夠反映植物從土壤中吸取養(yǎng)分的潛力,主要包括土壤的機(jī)械組成、有機(jī)質(zhì)含量、氮、磷和鉀的全量和有效量,這些指標(biāo)對(duì)樹木的生長(zhǎng)有著重要影響。不同林分由于樹種、造林模式、郁閉度等組成條件的不同,會(huì)導(dǎo)致土壤理化性質(zhì)發(fā)生變化,進(jìn)而影響森林生態(tài)系統(tǒng)結(jié)構(gòu)變化,土壤抵抗鹽堿化的能力也隨之改變,并可能最終改變森林演替軌跡[5-6]。因此,深入研究不同林分類型的土壤肥力狀況,對(duì)評(píng)價(jià)土壤質(zhì)量狀況、科學(xué)利用土地資源、合理進(jìn)行森林經(jīng)營(yíng)管理以及維護(hù)森林生態(tài)系統(tǒng)平衡均具有重要意義。
目前常用的傳統(tǒng)土壤肥力評(píng)價(jià)方法主要分為經(jīng)驗(yàn)權(quán)重法和統(tǒng)計(jì)權(quán)重法[7-10]。不同評(píng)價(jià)方法采用的體系和側(cè)重點(diǎn)不同[11-12],其最終得出的結(jié)果也就有所不同。本文主要采用更側(cè)重于數(shù)據(jù)本身體現(xiàn)的統(tǒng)計(jì)學(xué)意義和數(shù)理統(tǒng)計(jì)方法的統(tǒng)計(jì)權(quán)重中的主成分分析法和相關(guān)性分析法。
為研究不同林分類型的土壤肥力狀況,本研究以內(nèi)蒙古自治區(qū)清水河縣公益林區(qū)內(nèi)的油松林、落葉松林、山杏林、喬灌混交林(油松+檸條)、針闊混交林(銀白楊+落葉松)5 種主要林分人工林林下土壤為研究對(duì)象,采用相關(guān)性分析法[10]以及主成分分析法[13]對(duì)土壤肥力進(jìn)行綜合評(píng)價(jià)。旨在探討對(duì)比不同林分類型下人工林土壤肥力的狀況。為黃土丘陵區(qū)生態(tài)公益林管理和對(duì)森林土壤資源的科學(xué)管理提供理論依據(jù)。
清水河縣位于內(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ā)量達(dá)2577.2 mm;年平均氣溫7.1℃。主要有栗鈣土、栗褐土、灰褐土、潮土、風(fēng)沙土、沼澤土、鹽土、石質(zhì)土8 個(gè)土類;地帶性土壤為栗鈣土。
在查閱文獻(xiàn)和實(shí)地調(diào)查的基礎(chǔ)上,于2020年7 ~9 月在清水河縣公益林區(qū)內(nèi),根據(jù)樹種分布現(xiàn)狀和數(shù)量,篩選出立地條件相同、具有代表性的5 種典型林分類型,包括油松林、落葉松林、山杏林、喬灌混交林(油松+檸條)、針闊混交林(銀白楊+落葉松),林齡均在15 年以上,在其中選取具有代表性的區(qū)域,并在每種林分代表區(qū)域中隨機(jī)設(shè)置5 個(gè)樣地斑塊,共計(jì)25塊樣地,每塊樣地面積均為20 m×20 m,在樣地對(duì)角線上設(shè)置3 個(gè)采樣點(diǎn),每個(gè)采樣點(diǎn)面積為1 m×1 m;收集0 ~20 cm 土層的土壤,去除表層腐殖質(zhì)、殘根等,平均分成4 份共收集1 kg 的樣品,編號(hào)后帶回實(shí)驗(yàn)室,風(fēng)干后過篩備用。樣地基本情況見表1。
土壤pH 值采用混合指示劑比色法測(cè)定;有機(jī)質(zhì)含量采用重鉻酸鉀-硫酸溶液氧化法測(cè)定;有效磷含量用NaHCO3浸提-鉬銻抗比色法測(cè)定;全磷、全氮含量采用HClO4-H2SO4消化法測(cè)定;速效鉀含量采用中性乙酸銨浸提法測(cè)定;堿解氮含量用堿解擴(kuò)散法測(cè)定[14]。
表1 樣地基本情況
1.4.1 指標(biāo)標(biāo)準(zhǔn)化處理
在進(jìn)行評(píng)價(jià)之前,因?yàn)楦黜?xiàng)待評(píng)價(jià)指標(biāo)的量綱數(shù)據(jù)具有差異性,因而需要解決其統(tǒng)一標(biāo)準(zhǔn)化的問題,進(jìn)而完成各項(xiàng)評(píng)價(jià)指標(biāo)的數(shù)據(jù)統(tǒng)一化,使各評(píng)價(jià)指標(biāo)間具備可比性[4]。因此,需要?jiǎng)?chuàng)建出各項(xiàng)評(píng)價(jià)指標(biāo)的隸屬度函數(shù),并計(jì)算各指標(biāo)的隸屬度函數(shù)值。按照相關(guān)文獻(xiàn)記載以及不同樹種的實(shí)際生長(zhǎng)情況,可將隸屬度函數(shù)分為拋物線型和S 型。
拋物線型隸屬度函數(shù):評(píng)價(jià)指標(biāo)體系中具備一定范圍的樹木最佳生長(zhǎng)區(qū)間。假設(shè)其指標(biāo)數(shù)據(jù)離開此區(qū)間,則認(rèn)定不適于樹木的生長(zhǎng)。指標(biāo)數(shù)據(jù)的偏離程度可以反映其對(duì)樹木生長(zhǎng)的影響,pH 指標(biāo)離散程度較低,適用于此函數(shù)。隸屬度函數(shù)值F(x)的計(jì)算公式為:式中,x 為測(cè)定評(píng)價(jià)指標(biāo)的實(shí)際值,xa、xb、xc、xd為函數(shù)曲線轉(zhuǎn)折點(diǎn)的取值。
S 型隸屬度函數(shù):S 型隸屬度函數(shù)的曲線特征是,指標(biāo)假定在一固定數(shù)值區(qū)間內(nèi)增長(zhǎng)時(shí),可以反映出評(píng)價(jià)指標(biāo)對(duì)樹種的生長(zhǎng)具有較大影響,一旦低于或高于此區(qū)間,則評(píng)價(jià)指標(biāo)影響較小,適用于數(shù)據(jù)離散程度較大的土壤有機(jī)質(zhì)、全氮、全磷、有效磷、速效鉀、堿解氮含量。式中,x 為測(cè)定的評(píng)價(jià)指標(biāo)的實(shí)際值,xa、xb為函數(shù)曲線轉(zhuǎn)折點(diǎn)取值。
綜上,根據(jù)公式(1)和(2)可將指標(biāo)數(shù)據(jù)轉(zhuǎn)化為0.1 ~1.0 的無綱量化數(shù)值,數(shù)值越接近1.0 表明指標(biāo)對(duì)樹木成長(zhǎng)的影響越大。依據(jù)全國(guó)第二次土壤普查數(shù)據(jù)[15]和查閱相關(guān)的文獻(xiàn)[16-17]。函數(shù)曲線轉(zhuǎn)折點(diǎn)的取值如表2 所示。
表2 隸屬度函數(shù)的轉(zhuǎn)折點(diǎn)取值
1.4.2 綜合評(píng)價(jià)方法
采用下列2 種評(píng)價(jià)方法對(duì)5 種林分類型進(jìn)行土壤肥力綜合評(píng)價(jià)。
(1)主成分分析法:運(yùn)用SPSS 22.0 對(duì)評(píng)價(jià)指標(biāo)實(shí)際值進(jìn)行降維處理,重新組合成一組獨(dú)立的綜合指標(biāo),采用加權(quán)求和法計(jì)算公因子方差,進(jìn)一步計(jì)算各個(gè)公因子方差占公因子方差總和的比例[18],將其作為指標(biāo)的權(quán)重αi。土壤肥力綜合得分C 計(jì)算公式為:
式中,αi為第i 個(gè)土壤肥力指標(biāo)占總體的權(quán)重,n為指標(biāo)數(shù)量,F(xiàn)(i)為某一林分類型各項(xiàng)指標(biāo)的隸屬度值。
(2)相關(guān)性分析法:主要研究?jī)蓚€(gè)或兩個(gè)以上處于同等地位的隨機(jī)變量間的相關(guān)關(guān)系的統(tǒng)計(jì)分析方法。在存在多變量的情況下,由于變量之間存在著直接或間接的影響,所以需要消除其他變量的影響。在土壤肥力評(píng)價(jià)體系中,在控制其他指標(biāo)的影響后,可根據(jù)各指標(biāo)間的相關(guān)程度,建立偏相關(guān)系數(shù)矩陣,計(jì)算其平均值,求出該平均值占所有指標(biāo)平均值之和的百分比。最終求得單項(xiàng)肥力指標(biāo)在表征土壤肥力狀況中的貢獻(xiàn)率,即權(quán)重系數(shù)Wi[19],最終基于Wi和各指標(biāo)的隸屬度函數(shù)值可求得綜合指數(shù)IFI[20]:
式中,Wi為第i 個(gè)土壤肥力指標(biāo)的權(quán)重系數(shù),n 為指標(biāo)數(shù)量,F(xiàn)(i)為某一林分類型各項(xiàng)指標(biāo)的隸屬度值。
數(shù)據(jù)處理采用Excel 2010,使用拋物線型和S型曲線,對(duì)實(shí)際觀測(cè)數(shù)據(jù)進(jìn)行統(tǒng)一化處理,使用SPSS 22.0 進(jìn)行數(shù)據(jù)分析,對(duì)評(píng)價(jià)指標(biāo)進(jìn)行描述性統(tǒng)計(jì)及相關(guān)性分析,利用相關(guān)性分析法以及主成分分析法分析計(jì)算得出不同林分類型林下土壤肥力分值。
根據(jù)表3 顯示的各項(xiàng)土壤肥力指標(biāo)的描述性統(tǒng)計(jì)結(jié)果顯示,數(shù)據(jù)均呈正態(tài)分布,變異系數(shù)由大到小為速效鉀>全氮>堿解氮>有效磷>全磷>有機(jī)質(zhì)>pH。其中pH、有機(jī)質(zhì)為低于10%的弱變異,數(shù)據(jù)離散程度較小,說明5 種林分土壤pH、有機(jī)質(zhì)兩項(xiàng)指標(biāo)數(shù)據(jù)較為接近;其他6 項(xiàng)指標(biāo)均屬中等變異,速效鉀為最高,說明5 種林分的土壤速效鉀含量差異較大。
表3 土壤指標(biāo)描述性統(tǒng)計(jì)
不同林分類型的土壤肥力指標(biāo)測(cè)定值如表4 所示,土壤pH 平均值為8.28,表明研究區(qū)土壤為堿性土壤。根據(jù)全國(guó)第二次土壤普查分類標(biāo)準(zhǔn)[16](表5),5 種林分類型中,5 種林分的土壤全磷、有效磷含量均為三級(jí),屬于中上水平;堿解氮含量平均為68.19 ~88.66 mg·kg-1,均為四級(jí)中下水平,整體含量偏低;落葉松林、山杏林的有機(jī)質(zhì)、全氮含量為四級(jí),油松林、喬灌混交林、針闊混交林為三級(jí);針闊混交林的速效鉀含量為二級(jí)較高水平,油松林、山杏林、喬灌混交林為三級(jí),落葉松林為四級(jí)。
表4 不同林分類型土壤肥力指標(biāo)測(cè)定值
表5 全國(guó)第二次土壤普查養(yǎng)分分級(jí)標(biāo)準(zhǔn)
由表4 可得,除pH 外,以落葉松+楊樹混交為主要樹種的針闊混交林的各項(xiàng)肥力指標(biāo)數(shù)值均高于其他4 種林分,土壤肥力狀況較其他林分更好,5 種林分中,山杏林肥力狀況為最差。且數(shù)據(jù)顯示,混交林的土壤肥力要優(yōu)于純林。5 種林分的土壤全磷、有效磷含量整體較豐富,其中針闊混交林與山杏林的有效磷含量差異顯著,差值達(dá)到4.69 mg·kg-1。速效鉀含量整體差異較大,平均值為156.16 ~97.65 mg·kg-1,其中油松林、山杏林和喬灌混交林之間無顯著差異,落葉松林與針闊混交林差異顯著。5 種林分有機(jī)質(zhì)、全氮和堿解氮含量平均值分 別 為18.06 ~19.84 g·kg-1、0.82 ~1.00 g·kg-1和68.19 ~88.66 mg·kg-1,整體無顯著差異。
根據(jù)雷達(dá)圖的幾何意義,坐標(biāo)點(diǎn)離原點(diǎn)越遠(yuǎn)表明指標(biāo)的狀態(tài)越好,各坐標(biāo)點(diǎn)圍成的面積越大,則表明其評(píng)價(jià)對(duì)象的整體狀態(tài)更佳。因此,雷達(dá)圖可更加直觀地反映出單項(xiàng)土壤肥力指標(biāo)在土壤中的狀態(tài)及土壤肥力的整體情況[3],5 種林分各項(xiàng)土壤肥力指標(biāo)隸屬度函數(shù)值如圖1 所示。5 種林分類型各項(xiàng)土壤肥力指標(biāo)的隸屬度值最小的均為pH,隸屬度值基本在0.5 左右,說明pH 對(duì)林下土壤質(zhì)量的作用分值較小。隸屬度值最高為全氮,隸屬度值基本在0.7 ~0.8 之間,作用分值較大。其他指標(biāo)隸屬度值基本在0.55 ~0.7 之間,處于中間狀態(tài)。其中,最高的pH 隸屬度值出現(xiàn)在落葉松林,最高的有機(jī)質(zhì)、全磷、全氮、有效磷、速效鉀和堿解氮隸屬度值均出現(xiàn)在針闊混交林,說明針闊混交林的整體土壤肥力狀況較其他4 種林分更好。
圖1 中各項(xiàng)指標(biāo)隸屬度值圍合而成的面積能夠直觀反映出各項(xiàng)指標(biāo)對(duì)土壤肥力的作用影響,總體上,面積越大,林下土壤肥力情況越好。但是,以上數(shù)據(jù)僅是依據(jù)土壤肥力指標(biāo)反映了土壤肥力的狀態(tài)且是建立在土壤肥力相等的前提下,實(shí)際上,由于指標(biāo)間對(duì)土壤肥力的影響具有差異性,有必要結(jié)合單一土壤肥力指標(biāo)的綜合平均貢獻(xiàn)率進(jìn)一步驗(yàn)證。
圖1 各項(xiàng)土壤肥力指標(biāo)隸屬度函數(shù)值雷達(dá)圖
對(duì)不同林分類型的土壤肥力指標(biāo)隸屬度值進(jìn)行相關(guān)分析,并通過相關(guān)系數(shù)逆矩陣求得各指標(biāo)的權(quán)重系數(shù),進(jìn)而可知各項(xiàng)土壤肥力指標(biāo)對(duì)土壤肥力的貢獻(xiàn)率,結(jié)果如圖 2 所示。不同林分類型各項(xiàng)指標(biāo)綜合平均貢獻(xiàn)率從大到小為:全氮(0.199)>全磷(0.163)>有效磷(0.108)>堿解氮(0.099)>有機(jī)質(zhì)(0.090)>速效鉀(0.088)>pH(0.078)。由此可以看出,全氮與全磷對(duì)改善土壤肥力有著重要作用。
圖2 各項(xiàng)土壤肥力指標(biāo)隸屬度函數(shù)值綜合平均貢獻(xiàn)率雷達(dá)圖
2.3.1 主成分分析法
運(yùn)用SPSS 22.0 對(duì)7 項(xiàng)土壤肥力指標(biāo)進(jìn)行檢驗(yàn),KMO 測(cè) 驗(yàn) 結(jié) 果 為0.611(>0.6),Bartlett 球 形檢驗(yàn)的顯著系數(shù)為0.00(<0.05),說明數(shù)據(jù)適宜進(jìn)行主成分分析。經(jīng)過主成分分析得到各主成分?jǐn)?shù)據(jù)值(表6),共提取出3 個(gè)主成分,各主成分的特征值分別為3.152、1.360 和1.054,前3 個(gè)主成分的累計(jì)方差貢獻(xiàn)率分別達(dá)39.562%、63.740%和79.510%。第1 主成分包括有機(jī)質(zhì)、全磷、有效磷和堿解氮,第2 主成分包括pH 和速效鉀,第3 主成分為全氮。基于公因子方差求出的權(quán)重系數(shù)表現(xiàn)為全氮(0.161)>有效磷(0.159)>全磷(0.153)>速效鉀(0.139)>堿解氮(0.135)>有機(jī)質(zhì)(0.127)>pH(0.123),全氮對(duì)土壤肥力的貢獻(xiàn)率最大,其次為有效磷和全磷。
按照各林分類型的隸屬度函數(shù)值(圖1)以及求出各指標(biāo)的權(quán)重(表6),通過對(duì)其加權(quán)計(jì)算,可得到的土壤肥力綜合評(píng)價(jià)值如表7 所示,綜合得分為0.586 ~0.737,大小順序?yàn)獒橀熁旖涣郑締坦嗷旖涣郑居退闪郑韭淙~松林>山杏林,針闊混交林分值最高,山杏林最小。
2.3.2 相關(guān)性分析法
基于各項(xiàng)指標(biāo)的隸屬度函數(shù)值(圖1)和各指標(biāo)對(duì)土壤肥力的貢獻(xiàn)率(圖2)進(jìn)行加權(quán)計(jì)算,對(duì)不同林分類型各項(xiàng)土壤肥力指標(biāo)進(jìn)行綜合評(píng)價(jià),結(jié)果如表7 所示。綜合指數(shù)為0.495 ~0.629,與主成分分析法結(jié)果一致。大小順序?yàn)獒橀熁旖涣郑締坦嗷旖涣郑居退闪郑韭淙~松林>山杏林,針闊混交林肥力狀況最好,山杏林最差。
綜上,兩種評(píng)價(jià)方法進(jìn)行土壤肥力綜合評(píng)價(jià)的最終結(jié)果一致,均為針闊混交林>喬灌混交林>油松林>落葉松林>山杏林,說明數(shù)據(jù)本身對(duì)兩種方法體現(xiàn)出統(tǒng)一性,更具有說服力,且兩種方法適用于土壤肥力的綜合評(píng)價(jià)。二者最終得出的評(píng)價(jià)分值實(shí)際較為接近,但指標(biāo)權(quán)重有一定差別,例如相關(guān)性分析法得出的全磷平均綜合貢獻(xiàn)率在所有指標(biāo)中排名第2,為0.163,但主成分分析法得出的權(quán)重為0.153,在指標(biāo)中位列第3,說明不同評(píng)價(jià)方法偏重的內(nèi)容有所不同。主成分分析法主要是通過降維將原來具有一定相關(guān)性的變量重新組合成一組新的相互無關(guān)的幾個(gè)綜合變量,在盡可能多地保留原始變量信息的情況下,從中導(dǎo)出少數(shù)幾個(gè)主成分[21];而相關(guān)性分析法更多的是考察多個(gè)變量間的相關(guān)性,通過指標(biāo)之間的敏感度大小來確定指標(biāo)權(quán)重[3]。所以,在方法的選擇上,還應(yīng)根據(jù)指標(biāo)體系的具體情況來點(diǎn)定。
表6 土壤肥力指標(biāo)的主成分分析
表7 不同林分類型的土壤肥力綜合評(píng)價(jià)
本研究結(jié)果表明,研究區(qū)內(nèi)蒙古清水河縣生態(tài)公益林區(qū)土壤pH 平均值為8.35,屬于堿性土壤,這與以往關(guān)于黃土丘陵區(qū)土壤pH 值研究結(jié)果一致[22]。其中,混交林的土壤肥力狀況要高于純林?;旖涣种嗅橀熁旖涣郑▎棠净旖唬┩寥婪柿?yōu)于喬灌混交林,純林中油松林>落葉松林>山杏林。
混交林土壤肥力表現(xiàn)較純林更好,基本與前人研究結(jié)果一致[23-24],分析原因是由于人為對(duì)林分結(jié)構(gòu)進(jìn)行改造,改變了樹木的空間分布格局[25],使其能夠更加合理地吸收土壤養(yǎng)分,使土壤通透性增加,土壤結(jié)構(gòu)發(fā)生改變,更有利于養(yǎng)分的積累[26];且混交林林內(nèi)枯落物較純林更為豐富,成分復(fù)雜,為林下土壤微生物的活動(dòng)提供有利條件[27],使得土壤微生物更為活躍,加速了枯落物的分解和養(yǎng)分釋放的速度,改良了土壤理化性質(zhì),進(jìn)而提高了土壤肥力。5 種林分中,山杏林土壤肥力表現(xiàn)最差,原因主要是山杏成長(zhǎng)速度比其余樹種快,在生長(zhǎng)過程中對(duì)營(yíng)養(yǎng)的需求量較大,且林下植被稀少,最終導(dǎo)致枯落物蓄積量較少,只是表層養(yǎng)分積累,土壤整體肥力狀況不佳。但也有研究表明,山杏林土壤肥力狀況較油松等針葉林更好[28],推測(cè)可能與當(dāng)?shù)負(fù)嵊芾泶胧┯嘘P(guān)。綜上所述,在植被恢復(fù)和撫育管理人工林的過程中,在遵從自然規(guī)律以及因地制宜的原則上,可以適當(dāng)增加混交林所占比例,引入適宜的混交樹種,并根據(jù)土壤肥力狀況進(jìn)行科學(xué)的土壤養(yǎng)分管理。
由主成分分析法計(jì)算出的單項(xiàng)指標(biāo)權(quán)重及雷達(dá)圖所呈現(xiàn)的貢獻(xiàn)率均說明全氮在當(dāng)?shù)赝寥婪柿χ胸暙I(xiàn)最大,其含量的大小與土壤養(yǎng)分環(huán)境密切相關(guān)[29],可在今后作為指示土壤肥力的第一因素。pH 是土壤化學(xué)性質(zhì)的綜合體現(xiàn)[30-31],對(duì)土壤養(yǎng)分的存在形態(tài)有直接影響[32]。本研究表明,pH 是土壤肥力的第一限制性因子,因此應(yīng)注意調(diào)節(jié)土壤pH,可在土壤中適當(dāng)添加以石膏或磷石膏為主的土壤改良劑用以降低土壤堿性[33]。