袁晶晶，同延安※，盧紹輝，袁國軍

（1. 西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院，楊凌 712100；2. 農(nóng)業(yè)部西北植物營養(yǎng)與農(nóng)業(yè)環(huán)境重點實驗室，楊凌 712100；3. 河南省林業(yè)科學(xué)研究院，鄭州 450003）

0 引 言

土壤肥力是土壤質(zhì)量的核心，也是農(nóng)業(yè)可持續(xù)發(fā)展的基礎(chǔ)。合理進(jìn)行土壤培肥是維持土壤質(zhì)量和保證土壤資源可持續(xù)利用的最重要措施之一[1]。生物炭特有的結(jié)構(gòu)和性質(zhì)（如孔隙豐富、比表面積大），使其作為一種土壤改良劑受到很多關(guān)注[2]，并廣泛應(yīng)用于農(nóng)業(yè)土壤、環(huán)境生態(tài)等領(lǐng)域[3-5]。

華北平原作為中國主要的紅棗產(chǎn)區(qū)之一，由于不合理的管理措施，施肥模式尤其是對化肥和農(nóng)藥的過度依賴，導(dǎo)致土壤酸化、肥料利用率低及生產(chǎn)效益降低等諸多問題，不僅造成了資源浪費，也帶來了一系列的環(huán)境問題[6]，所以培肥地力仍是實現(xiàn)該地區(qū)產(chǎn)業(yè)持續(xù)發(fā)展的關(guān)鍵措施。不少研究表明，生物炭與化學(xué)肥料配合施用不僅能減少土壤養(yǎng)分的淋失，延緩養(yǎng)分的釋放，從而提高作物產(chǎn)量，提高肥料利用率[7-12]。但也有研究顯示，生物炭對土壤肥力和碳庫存潛力無影響或有抑制作用[13-14]。這些差異結(jié)果與生物炭種類、熱解溫度、施用量、土壤類型以及其它不確定環(huán)境因素有關(guān)。因此，客觀的、全面的綜合評價生物炭和其它化學(xué)肥料配合施用對土壤肥力的績效，并找出其對土壤質(zhì)量影響的驅(qū)動因子，尋找最佳互作模式，在提高經(jīng)濟效益的同時，達(dá)到改良土壤質(zhì)量、提高肥料利用率的作用。但目前，利用統(tǒng)計學(xué)方法綜合評價生物炭與氮肥不同配比條件下對棗園土壤質(zhì)量的相關(guān)研究仍較為鮮見。本研究采用因子分析與聚類分析法，選取具有代表性的土壤理化及生物學(xué)14個指標(biāo)進(jìn)行綜合分析研究，明確生物炭與氮肥不同配比條件下的應(yīng)用效果，探尋最佳配比施肥量，以期為棗園土壤資源高效利用、培肥改良技術(shù)提供可靠的科學(xué)依據(jù)。

1 材料與方法

1.1 試驗地概況

田間試驗于2013年4月～2015年10月在河南省濮陽市王助鄉(xiāng)潘莊進(jìn)行。潘莊（114°52′ E, 35°20′ N），屬溫帶大陸性氣候，四季分明，光熱資源適中。年平均氣溫13.4℃，年平均降水量為502.3～601.3 mm，年平均日照時數(shù)為2 454.5 h，無霜期一般為205 d。試驗前耕層土壤（0～20 cm）基本理化指標(biāo)如下：pH值8.47、容重1.42 g/cm3、有機質(zhì)質(zhì)量分?jǐn)?shù)10.61 g/kg、全氮0.68 g/kg、速效磷11.67 mg/kg、速效鉀95.61 mg/kg。

1.2 試驗設(shè)計

試驗采用隨機區(qū)組設(shè)計，每個小區(qū)面積30 m2，每個處理重復(fù) 5次。目前中國生物質(zhì)炭在大田的施用量介于3～40 t/hm2之間[15]，結(jié)合當(dāng)?shù)貙嶋H情況，試驗采用4 × 3完全方案設(shè)計，設(shè)生物炭用量4個水平（即C0、C1、C2和C3）、氮肥用量3個水平（N1、N2和N3），加上1個完全空白處理CK（不施生物炭和氮肥），共計13個處理，具體施肥用量詳見表1。種植品種為15 a生扁核酸紅棗，栽植密度為2 m×3 m，約1 650棵/hm2，設(shè)置保護(hù)行。

表1 田間試驗各處理生物炭和氮肥用量Table 1 Biochar and nitrogen fertilizer amount under different treatments in field experiment

研究所用生物炭原料為花生殼秸稈，購自于河南商丘三利新能源有限公司，生產(chǎn)設(shè)備采用連續(xù)豎式生物質(zhì)炭化爐，炭化溫度為350～500℃。其理化性質(zhì)如下：pH值10.35、有機碳質(zhì)量分?jǐn)?shù)467.20 g/kg、全氮5.90 g/kg、全鉀11.5 g/kg、CEC 21.7 mmol/kg、容重0.45 g/cm3。

過磷酸鈣（P2O5質(zhì)量分?jǐn)?shù)16%）300 kg/hm2，硫酸鉀（K2O質(zhì)量分?jǐn)?shù)45%）300 kg/hm2均作為基肥在春季一次性施入，尿素（質(zhì)量分?jǐn)?shù) 70%）作為基肥施入，剩下的30 %在7月中旬作為追肥施入。在樹冠下兩側(cè)0.5 m左右挖深20～30 cm的條狀溝，使生物炭與土混勻后覆土填平，此后的2 a生物炭不再施入。在作物生長期間根據(jù)天氣及作物生長狀況（不同生育期）適量灌水，約 2 250 m3/hm2，以滿足作物正常生長發(fā)育所需。

1.3 樣品采集及測定

土壤樣品于2015年10月紅棗采收后，每個小區(qū)用土鉆在耕層( 0～20 cm) 分別按照“S”形取5個點，剔除雜物?；旌暇鶆蚝笥盟姆址ǚ殖?份，一份鮮樣過2 mm篩4℃保存供土壤微生物量碳、氮、土壤酶活性以及微生物數(shù)量的測試，另一份風(fēng)干供理化性質(zhì)測試。

土壤微生物量碳、氮采用氯仿熏蒸浸提方法測定[16]；土壤酶活性測定[17-18]：脲酶測定采用苯酚鈉比色法，堿性磷酸酶測定采用磷酸苯二鈉比色法，蔗糖酶采用Na2S2O3滴定法。土壤微生物數(shù)量的測定[19]：細(xì)菌采用牛肉膏蛋白胨培養(yǎng)基，放線菌采用改良高氏一號培養(yǎng)基。

土壤基本性狀測定方法[20]：土壤有機質(zhì)采用重鉻酸鉀-濃硫酸外加熱法，全氮采用凱氏定氮法；全磷、全鉀采用NaOH熔融-鉬銻抗比色法、火焰光度計法測定；土壤速效氮1 mol/L KCL浸提-流動分析儀法；速效磷采用0.5 mol/L NaHCO3浸提-鉬銻抗比色法；速效鉀采用1 mol/L NH4OAC浸提-火焰光度計法測定。

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

試驗數(shù)據(jù)采用 Excel 2016、SPSS22.0（IBM corporation,Aromonk, New York）統(tǒng)計軟件進(jìn)行方差分析與原始數(shù)據(jù)的標(biāo)準(zhǔn)化，選擇LSD最小顯著檢驗方法（least significant difference）進(jìn)行顯著性多重比較。同時利用因子分析（factor analysis）和聚類分析（hierarchical cluster analysis）分別對土壤性質(zhì)進(jìn)行因子分析和聚類分析。

1.5 因子分析和聚類分析的基本原理

因子分析是指通過顯在變量測評潛在變量，在多個變量中找出少數(shù)具有代表性的綜合因子，減少變量數(shù)目，達(dá)到降維目的[21-22]。因子分析法具有命名清晰性高和應(yīng)用上側(cè)重成因清晰性的綜合評價的特點，所提取的公因子比主成分分析提取的主成分更具有解釋性。

聚類分析是指將物理或抽象對象的集合分組成為由類似的對象組成的多個類的分析過程。以歐氏距離作為衡量各處理肥力差異大小，采用最短距離法將各處理按土壤肥力水平的親疏相似程度進(jìn)行系統(tǒng)聚類[23]。

2 結(jié)果與分析

2.1 生物炭與氮肥配施措施下土壤肥力評價指標(biāo)的選取

為了更加全面客觀地評價土壤培肥效果，需要對土壤肥力指標(biāo)進(jìn)行篩選，使評價指標(biāo)更加科學(xué)化、合理化[24]。通常選取土壤理化和生物指標(biāo)作為評價指標(biāo)[25-26]，且遵循代表性、主導(dǎo)性、穩(wěn)定性、差異性、可比性等原則[27]。本研究根據(jù)以上原則且結(jié)合試驗區(qū)實際指標(biāo)的測定，選取14項指標(biāo)進(jìn)行土壤肥力質(zhì)量評價, 具體測定值見表2。

2.2 土壤肥力質(zhì)量的因子分析

由于本試驗中14個土壤指標(biāo)具有不同的量綱，且它們的數(shù)量級差別較大，為確保分析結(jié)果的客觀性和科學(xué)性，需要對研究所選取的各土壤指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理。同時，對標(biāo)準(zhǔn)化后的數(shù)據(jù)進(jìn)行KMO和Bartlett檢驗，KMO的取值為0.85，Sig.值為0.000 1＜0.001，符合因子分析的前提條件。依據(jù)因子分析方法的原理[28]，運用SPSS 22.0統(tǒng)計軟件計算出各指標(biāo)的相關(guān)系數(shù)矩陣（見表3），分析得出的特征值與貢獻(xiàn)率如表 4所示。特征值在一定程度上可看作公共因子影響原變量力度大小的指標(biāo)，一般來說，特征值＞1能較好地說明該公因子對原變量的解釋力度。本試驗可提取出 3個公因子，其累計貢獻(xiàn)率分別達(dá)到82.49%，即說明前3個公因子涵蓋了原始數(shù)據(jù)信息的82.49%，能較明確地解釋該地區(qū)土壤肥力水平。因此，將這 3個公因子作為綜合變量來評價本試驗的土壤肥力狀況是可行的[22]。

表2 土壤肥力各指標(biāo)平均值Table 2 Mean values of soil fertility indices

表3 土壤各肥力指標(biāo)的相關(guān)系數(shù)Table 3 Correlation coefficient of the soil properties

對本試驗選用的 3個公因子進(jìn)行分析，通過方差最大化旋轉(zhuǎn)后的因子載荷矩陣見表5。由表5可知，第1公因子上，具有高載荷因子的指標(biāo)有全氮、全鉀、速效鉀、微生物量碳，其中全氮、全鉀反映了土壤養(yǎng)分的總儲量，微生物量碳作為土壤肥力的評價指標(biāo)，可反映土壤的有效養(yǎng)分狀況及生物活性[29]?？梢?，第 1公因子實質(zhì)上反映了土壤環(huán)境和土壤養(yǎng)分的貯存和轉(zhuǎn)化，可作為保肥供肥因子。第 2公因子上脲酶、微生物量氮、蔗糖酶、有機質(zhì)、速效氮和全磷為主要影響因子。堿性磷酸酶為第3公因子上的主要影響因子（表5）。3個公因子涵蓋了土壤的化學(xué)和生物指標(biāo)，使評價更具科學(xué)性、合理性[30-31]。

表5 旋轉(zhuǎn)成分矩陣Table 5 Rotated component matrix

2.3 不同處理土壤質(zhì)量的得分及排名

為了更直觀地比較生物炭與氮肥不同用量處理下土壤肥力質(zhì)量狀況，通過公因子得分系數(shù)矩陣（表6），計算出各處理的因子得分和綜合得分（表 7）。在F1公因子上C3N3、C3N1和C2N1處理的得分位于前3名，說明其在全氮、全鉀、速效鉀和微生物數(shù)量碳上具有相對較高的優(yōu)勢，這 3個處理土壤養(yǎng)分保肥供肥能力強；在公因子F2上，C2N3、C2N1和C2N23個處理的得分位于前3名。說明生物炭與氮肥配施增加了土壤脲酶、蔗糖酶的活性以及微生物量氮的含量，使土壤供給作物養(yǎng)分能力增加。

通過把各因子的特征值貢獻(xiàn)率作為權(quán)數(shù)進(jìn)行加權(quán)求和，就可得到綜合評價得分。F=F1×49.618%+F2×25.420%+F3×7.450%，結(jié)果見表7。棗區(qū)中以C3N1處理的綜合得分最高，其次為C2N3和C2N2處理。該試驗區(qū)綜合得分最低的均為CK，其次為不施生物炭的3個氮肥處理。

表6 成分得分系數(shù)矩陣Table 6 Scores of component matrix

表7 不同處理各因子得分及綜合得分Table 7 Scores and general scores of principal components under different treatments

作物產(chǎn)量是土壤內(nèi)在屬性外在的間接綜合表現(xiàn)，產(chǎn)量的高低在一定程度上可以直觀反映土壤質(zhì)量的好壞[32]，因此也常被作為檢驗評價結(jié)果客觀性及準(zhǔn)確性的依據(jù)。由圖1可知，C3N1處理的扁核酸紅棗產(chǎn)量最高，為22 716.4 kg/hm2，其次為C2N2和C2N3處理，不同處理土壤質(zhì)量的綜合得分與產(chǎn)量的變化趨勢比較吻合。

圖1 生物炭與氮肥配施對扁核酸紅棗產(chǎn)量的影響Fig.1 Effects of biochar and nitrogen fertilizer amendments on jujube (Bianhesuan) yields

2.4 聚類分析

從圖 2可以看出各配施處理土壤質(zhì)量狀況的大致分類，按土壤肥力水平的親疏相似程度進(jìn)行系統(tǒng)聚類。為體現(xiàn)生物炭與氮肥不同配比施肥對土壤肥力影響的差異，將距離閾值定為 8，此時類間距較大，13 個處理大致可以分為5類，即CK、C0N1、C0N2和C0N34個處理為1個等級，C1N1、C1N2和C1N3處理為1個等級，C2N2、C2N3、C3N1和C3N3為1 個等級，C2N1和C3N2分別各為一類。其中CK為代表的這一等級的土壤肥力程度最低，C3N1、C2N2、C2N3和C3N3這一等級屬于高肥力水平。說明在本研究試驗條件下，生物炭和氮肥配施可不同程度地提高棗園土壤肥力水平，且生物炭和氮肥用量是影響土壤肥力水平高低的關(guān)鍵因素。

圖2 不同施肥處理下土壤質(zhì)量評價聚類圖Fig.2 Cluster graph of soil quality assessment in different fertilizer treatments

3 討 論

土壤質(zhì)量是土壤物理、化學(xué)性質(zhì)和生物學(xué)特性的綜合反映，其評價結(jié)果可直觀地反映土壤質(zhì)量總體狀況[33]。由于評價目的和側(cè)重評價的尺度不一樣，評價指標(biāo)和評價方法的選擇亦有所差異[34]，但到目前為止，土壤質(zhì)量的評價方法國際上尚沒有統(tǒng)一的標(biāo)準(zhǔn)，也沒有固定的方法[35]。一般來說，評價土壤質(zhì)量需要土壤物理、化學(xué)和生物指標(biāo)。近年來，國內(nèi)外關(guān)于生物炭的應(yīng)用開展了大量的研究工作[9,36-38]。本研究3a田間試驗表明，生物炭與氮肥配施顯著提高了土壤養(yǎng)分含量、微生物量、土壤酶活性及微生物數(shù)量，在一定程度上改善了土壤的理化和生物學(xué)性質(zhì)，提高了土壤質(zhì)量，進(jìn)而利于作物高產(chǎn)（表2）。同時相關(guān)性分析表明，土壤養(yǎng)分與土壤微生物量碳氮及微生物數(shù)量之間存著不同程度的相關(guān)性（表3），這與宋大利等[39]研究結(jié)果具有一定的相似性，說明土壤理化、生物學(xué)性狀之間是彼此聯(lián)系，相互作用，共同影響土壤質(zhì)量水平。本文選取代表土壤質(zhì)量的土壤理化性質(zhì)和生物學(xué)性質(zhì)作為評價指標(biāo)，采用因子分析對生物炭與氮肥配施對棗園土壤質(zhì)量進(jìn)行綜合評價，將14個原始指標(biāo)降維、提取出3個公因子，累計貢獻(xiàn)率達(dá)82.49%。通過聚類分析對生物炭與氮肥不同配比處理下的土壤進(jìn)行分級，將13個處理分為5類不同肥力等級。結(jié)合土壤質(zhì)量綜合評價得分來看，10 t/hm2的生物炭，配施300 kg/hm2的氮肥培肥模式在改善華北平原棗園土壤結(jié)構(gòu)和提高作物產(chǎn)量方面優(yōu)勢更突出。

因子分析和聚類分析結(jié)果可以較好地反映實際土壤質(zhì)量及變化趨勢，對于生物炭和氮肥配施在農(nóng)田果園的科學(xué)管理和合理施肥具有重要意義。目前有關(guān)生物炭與氮肥配施土壤質(zhì)量評價方面的研究尚少，還有待進(jìn)一步研究驗證。

4 結(jié) 論

1）通過對生物炭與氮肥不同配比條件下的 14個土壤生化指標(biāo)進(jìn)行因子分析，提取出的 3個公因子，反映了原信息總量的82.49%。第1公因子以全氮、全鉀、速效鉀等貢獻(xiàn)最大，第 2公因子以脲酶、微生物量氮、蔗糖酶貢獻(xiàn)最大，堿性磷酸酶為第 3公因子上的主要影響因子。3個公因子涵蓋了土壤的生化指標(biāo)，使評價更具科學(xué)性。

2）土壤肥力指標(biāo)相關(guān)性分析表明，土壤養(yǎng)分、微生物量、土壤酶活性和微生物數(shù)量之間存在多種顯著或極顯著正相關(guān)關(guān)系。合理的生物炭與氮肥配施有利于改善棗園土壤理化性質(zhì)和微生物學(xué)特性，進(jìn)而提高土壤肥力和紅棗產(chǎn)量。

3）聚類分析將13 個不同施肥處理分為5類，其分析結(jié)果與各處理因子綜合得分評價較為一致，且生物炭與氮肥對紅棗產(chǎn)量的影響與土壤質(zhì)量因子的綜合得分也比較吻合，其中 C3N1處理的扁核酸紅棗產(chǎn)量最高，為22 716.4 kg/hm2。綜合本文試驗結(jié)果，10 t/hm2的生物炭，配施300 kg/hm2的氮肥為此地區(qū)最佳培肥模式。這一結(jié)果對生物炭和氮肥配施在農(nóng)田果園的科學(xué)管理和合理施肥具有指導(dǎo)意義。

[1] 鄭昭佩，劉作新. 土壤質(zhì)量及其評價[J]. 應(yīng)用生態(tài)學(xué)報，2003，14(1)：131–134.Zheng Zhaopei, Liu Zuoxin. Soil quality and its evaluation[J].Chinese Journal of Applied Ecology, 2003, 14(1): 131－134.(in Chinese with English abstract)

[2] 王丹丹，鄭紀(jì)勇，顏永毫，等. 生物炭對寧南山區(qū)土壤持水性能影響的定位研究[J]. 水土保持學(xué)報，2013，27(2)：101－104.Wang Dandan, Zheng Jiyong, Yan Yonghao, et al. Effect of biochar application on soil water holding capacity in the southern region of Ningxia[J]. Journal of Soil and Water Conservation, 2013, 27(2): 101－104. (in Chinese with English abstract)

[3] Glaser B, Haumaier L, Guffenberger G, et al. Black carbon indensity fractions of anthropogenic soils of the Brazilian Amazon region[J]. Organic Geochemistry, 2000, 31(7): 669－678.

[4] Glaser B, Haumaier L, Guggenberger G, et al. The “Terra Preta” phenomenon: A model for sustainable agriculture in the humid tropic[J]. Naturwissenschaften, 2001, 88(1): 37－41.

[5] Lehmann J, Gaunt J, Rondon M. Biochar sequestration in terrestrial ecosystems-A review[J]. Mitigation and Adaptation Strategies for Global Change, 2006, 11(2): 403－427.

[6] 黨維勤，鄭妍，卜曉鋒，等. 談黃土丘陵溝壑區(qū)紅棗產(chǎn)業(yè)的發(fā)展[J]. 中國水土保持，2007，(5)：52－54.Dang Weiqin, Zheng Yan, Bu Xiaofeng, et al. Discussion on the development of chinese date industry in gullied rolling loess area[J]. Soil and Water Conservation in China, 2007,(5): 52–54. (in Chinese with English abstract)

[7] Bass A M, Bird M I, Kay G, et al. Soil properties, greenhouse gas emissions and crop yield under compost, biochar and co-composted biochar in two tropical agronomic systems[J].Science of the Total Environment. 2016, 550: 459－470.

[8] Haider G, Steffens D, Moser G, et al. Biochar reduced nitrate leaching and improved soil moisture content without yield improvements in a four-year field study[J]. Agriculture,Ecosystem & Environment. 2017, 237: 80－94.

[9] Gul S, Whalen, J K. Biochemical cycling of nitrogen and phosphorus in biochar-amended soils[J]. Soil Biology &Biochemistry. 2016, 103: 1－15.

[10] Li B, Fan C H, Zhang H, et al. Combined effects of nitrogen fertilization and biochar on the net global warming potential,greenhouse gas intensity and net ecosystem economic budget in intensive vegetable agriculture in southeastern China[J].Atmospheric Environment. 2015, 100: 10－19.

[11] Herath H M S K, Camps-Arbestain M, Hedley M. Effect of biochar on soil physical properties in two contrasting soils:An Alfisol and an Andisol[J]. Geoderma, 2013, (209–210):188－197.

[12] Ippolito J A, Ducey T F, Cantrell K B, et al. Designer, acidic biochar influences calcareous soil characteristics[J]. Chemosphere,2016, 142: 184－191.

[13] Bruun E W, Ambus P, Egsgaard H, et al. Effects of slow and fast pyrolysis biochar on soil C and N turnover dynamic[J].Soil Biology& Biochemistry, 2012, 46(1): 73－79.

[14] Tammeorg P, Simojoki A, M?kel? P, et al. Short-term effects of biochar on soil properties and wheat yield formation with meat bone meal and inorganic fertiliser on a boreal loamy sand[J]. Agriculture Ecosystems & Environment, 2014, 191:108－116.

[15] 趙紅，孫濱峰，逯非，等. Meta分析生物質(zhì)炭對中國主糧作物痕量溫室氣體排放的影響[J]. 農(nóng)業(yè)工程學(xué)報，2017，33(19)：10－16.Zhao Hong, Sun Binfeng, Lu Fei, et al. Meta analysis on impacts of biochar on trace greenhouse gases emissions from staple in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017,33(19): 10－16. (in Chinese with English abstract)

[16] Jenkinson D S, Brookles P C, Powlson D S. Measuring soil microbial biomass[J]. Soil Biology and Biochemistry, 2004,36(1): 5－7.

[17] 關(guān)松蔭. 土壤酶及其研究方法[M]. 北京：農(nóng)業(yè)出版社，1986.Guan Songyin. Soil Enzyme and Its Research Methods[M].Beijing: Agriculture Press, 1986. (in Chinese with English abstract)

[18] 周禮愷. 土壤酶學(xué)[M]. 北京：科學(xué)出版社，1987.Zhou Likai. Soil Enzymology[M].Beijing: Science Press,1987. (in Chinese with English abstract)

[19] 袁麗紅. 微生物學(xué)實驗[M]. 北京：化學(xué)工業(yè)出版社，2010.Yuan Lihong. Laboratory Experiments in Microbiology[M].Beijing: Chemical Industry Press, 2010. (in Chinese with English abstract)

[20] 鮑士旦. 土壤農(nóng)業(yè)化學(xué)分析[M]. 北京：中國農(nóng)業(yè)出版社，2005.Bao Shidan. Soil and Agricultural Chemistry Analysis[M].Beijing: China Agriculture Press, 2005. (in Chinese with English abstract)

[21] 王芳. 主成分分析與因子分析的異同比較及應(yīng)用[J]. 統(tǒng)計教育. 2003，(5)：14－17.Wang Fang. Comparison and application of principal component analysis & factor analysis[J]. Statistic Education,2003, (5): 14-17. (in Chinese with English abstract)

[22] 高惠璇. 應(yīng)用多元統(tǒng)計分析[M]. 北京：北京大學(xué)出版社，2005.Gao Huixuan. Applied Multivariate Statistical Analysis[M].Beijing: Beijing University Press, 2005. (in Chinese with English abstract)

[23] 方開泰. 實用多元統(tǒng)計分析[M]. 上海：華東師范大學(xué)出版社，1989.Fang Kaitai. Practical Multivariate Statistical Analysis[M].Shanghai: East China Normal University Press, 1989. (in Chinese with English abstract)

[24] 黃勇，楊忠芳. 土壤質(zhì)量評價國外研究進(jìn)展[J]. 地質(zhì)通報，2009，28(1)：130－136.Huang Yong, Yang Zhongfang. Recent research progress of overseas soil quality evaluation[J]. Geological Bulletin of China, 2009, 28(1): 130－136. (in Chinese with English abstract)

[25] 劉曉冰，邢寶山，Stephen J H. 土壤質(zhì)量及其評價指標(biāo)[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，2002，18(2)：109-112.Liu Xiaobing, Xing Baoshan, Stephen J H. Soil quality and its assessment indicators[J]. System Science and Comprehensive Study in Agriculture, 2002, 18(2): 109－112. (in Chinese with English abstract)

[26] 路鵬，蘇以榮，牛錚，等. 土壤質(zhì)量評價指標(biāo)及其時空變異[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報，2007，15(4)：190－194.Lu Peng, Su Yirong, Niu Zheng, et al. Soil quality assessment indicators and their spatial temporal variability[J]. Chinese Journal of Eco-Agriculture, 2007, 15(4): 190－194. (in Chinese with English abstract)

[27] 周勇，田有國，任意，等. 定量化土地評價指標(biāo)體系及評價方法探討[J]. 生態(tài)環(huán)境，2003，12(1)：37－41.Zhou Yong, Tian Youguo, Ren Yi, et al. Research on index system and appraisal method of quantitative land assessment[J]. Ecology and Environment, 2003, 12(1): 37－41. (in Chinese with English abstract)

[28] Martin Y, Victor G J, David G R. Developing a minimum data set for characterizing soil dynamic in shifting cultivation system[J]. Soil & Tillage Research, 2006, 86(1): 84－98.

[29] Eivazi F, Bayan M R. Select soil enzyme activities in the historic sanborn field as affected by long-term cropping systems[J]. Communications in Soil Science and Plant Analysis, 2003, 34(15–16): 2259－2275.

[30] 郭篤發(fā)，王秋兵. 主成分分析法對土壤養(yǎng)分與小麥產(chǎn)量關(guān)系的研究[J]. 土壤學(xué)報，2005，42(3)：523－527.Guo Dufa, Wang Qiubing. Principal component analysis of relationship between various nutrients in albiudic cambosols profile and wheat yield[J]. Acta Pedologica Sinica, 2005,42(3): 523－527. (in Chinese with English abstract)

[31] Fox G A, Metla R. Soil property analysis using principal component sanalysis, soil line, and regression models[J]. Soil Science Society America Journal, 2005, 69: 1782－1788.

[32] 吳玉紅，田霄鴻，南雄雄，等. 基于因子和聚類分析的保護(hù)性耕作土壤質(zhì)量評價研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報，2010，18(2)：223－228.Wu Yuhong, Tian Xiaohong, Nan Xiongxiong, et al. Evaluation of soil quality under conservation tillage via factor and cluster analyses[J]. Chinese Journal of Eco-Agriculture, 2010, 18(2):223－228. (in Chinese with English abstract)

[33] 劉廣明，呂真真，楊勁松，等. 基于主成分分析及GIS的環(huán)渤海區(qū)域土壤質(zhì)量評價[J]. 排水機械工程學(xué)報，2015，33(1)：67－72.Liu Guangming, Lü Zhenzhen, Yang Jinsong, et al. Quality evaluation of soil in land Bohai Sea coast by principal component analysis and GIS[J]. Journal of Drainage and Irrigation Machinery Engineering, 2015, 33(1): 67－72. (in Chinese with English abstract)

[34] 貢璐，張海峰，呂光輝，等. 塔里木河上游典型綠洲不同連作年限棉田土壤質(zhì)量評價[J]. 生態(tài)學(xué)報，2011，31(14)：4136－4143.Gong Lu, Zhang Haifeng, Lü Guanghui, et al. Soil quality assessment of continuous cropping cotton fields for different years in a typical oasis in the upper reaches of the Tarim River[J]. Acta Ecologica Sinica, 2011, 31(14): 4136－4143.(in Chinese with English abstract)

[35] 崔東，肖治國，孫國軍，等. 伊犁河谷不同土地利用方式下土壤質(zhì)量評價[J]. 西北師范大學(xué)學(xué)報：自然科學(xué)版，2017，53(2)：112－117.Cui Dong, Xiao Zhiguo, Sun Guojun, et al. Soil quality evaluation of different land use in Yili river valley[J]. Journal of Northwest Normal University: Natural Science, 2017,53(2): 112－117. (in Chinese with English abstract)

[36] Ajayi A E, Holthusen D, Horn R. Changes in microstructural behavior and hydraulic functions of biochar amended soils[J].Soil and Tillage Research, 2016, 155: 166－175.

[37] Castellini M L, Giglio M, Niedda A D, et al. Impact of biochar addition on the physical and hydraulic properties of a clay soil[J]. Soil and Tillage Research, 2015, 154: 1－13.

[38] Hussain M, Farooq M, Nawaz A, et al. Biochar for crop production: Potential benefits and risks[J]. Journal of Soils and Sediments. 2016, 17(3): 658–716.

[39] 宋大利，習(xí)向銀，黃紹敏，等. 秸稈生物炭配施氮肥對潮土土壤碳氮含量及作物產(chǎn)量的影響[J]. 植物營養(yǎng)與肥料學(xué)報，2017，23(2)：369－379.Song Dali, Xi Xiangyin, Huang Shaomin, et al. Effect of combined application of straw biochar and nitrogen on soil carbon and nitrogen contents and crop yields in a fluvo-aquic soil[J]. Journal of Plant Nutrition Fertilizer, 2017, 23(2):369－379. (in Chinese with English abstract)