基于耕作指數(shù)評(píng)價(jià)耕作措施對(duì)東北風(fēng)沙土耕層土壤質(zhì)量的影響

苗淑杰，趙紅飛，喬云發(fā)*，陸欣春，李 琪

（1．南京信息工程大學(xué)應(yīng)用氣象學(xué)院，江蘇 南京 210044；2．中國(guó)科學(xué)院東北地理與農(nóng)業(yè)生態(tài)研究所，黑龍江 哈爾濱 150081）

土壤質(zhì)量是指土壤在生態(tài)系統(tǒng)的范圍內(nèi)，維持生產(chǎn)能力、保護(hù)環(huán)境質(zhì)量及促進(jìn)動(dòng)植物健康的綜合能力［1］。土壤質(zhì)量是一種指示土壤條件動(dòng)態(tài)變化最敏感的表示方法，反映土壤管理的變化［2］。改善土壤質(zhì)量對(duì)農(nóng)業(yè)生產(chǎn)有重要意義［3］。耕作措施影響土壤水分、養(yǎng)分的分布和遷移，是影響耕層土壤質(zhì)量變化最普通、最直接、最深刻的因素，合理耕作措施即提高土壤質(zhì)量,又改善生態(tài)環(huán)境［2，4-7］。東北風(fēng)沙土區(qū)，氣候干旱、土壤貧瘠、土壤保水保肥能力弱、耕層土壤質(zhì)量差，影響作物生長(zhǎng)［8］。因此，在東北風(fēng)沙土區(qū)農(nóng)田開(kāi)展耕作措施對(duì)耕層土壤質(zhì)量的影響研究具有重要意義。

土壤質(zhì)量是在土壤不同功能間尋求平衡和綜合而確定的土壤本身的內(nèi)在屬性，土壤質(zhì)量評(píng)價(jià)工作是根據(jù)已知的土壤外部性質(zhì)對(duì)土壤的這一內(nèi)在屬性進(jìn)行量化表達(dá)［9］。定量化分析耕作措施對(duì)土壤質(zhì)量的影響，根據(jù)量化評(píng)價(jià)的實(shí)際情況選擇適宜的耕作方式［7，10］。土壤質(zhì)量不能直接測(cè)定，目前國(guó)內(nèi)外對(duì)土壤質(zhì)量的評(píng)價(jià)主要通過(guò)選擇土壤理化性質(zhì)等土壤屬性作為評(píng)價(jià)指標(biāo)，建立評(píng)價(jià)模型，利用各種數(shù)學(xué)方法，根據(jù)量化的土壤屬性，計(jì)算出土壤質(zhì)量的“分?jǐn)?shù)”，最好的土壤得到最高的分?jǐn)?shù)［9，11］。土壤耕作指數(shù)是將特定土壤耕作下土壤質(zhì)量進(jìn)行定量描述的指數(shù)，通過(guò)計(jì)算不同耕作措施下的耕作指數(shù)（TI），評(píng)價(jià)耕作措施對(duì)土壤質(zhì)量的影響［7，12］。通過(guò)耕作指數(shù)對(duì)土壤生產(chǎn)力進(jìn)行評(píng)價(jià)，可為土壤生產(chǎn)力狀況和合理耕層構(gòu)建提供依據(jù)。徐尚起等［7］、Bockari-Gevao等［13］認(rèn)為利用土壤耕作指數(shù)定量評(píng)價(jià)土壤質(zhì)量，對(duì)于選擇合理耕作措施，提高耕層土壤質(zhì)量的可持續(xù)發(fā)展具有重要的意義。通過(guò)耕作與秸稈還田、有機(jī)肥施用集成模式是構(gòu)建合理耕層常用措施。合理耕作措施改善土壤理化學(xué)及生物學(xué)性狀，有利于農(nóng)田土壤質(zhì)量的提高，促進(jìn)農(nóng)業(yè)可持續(xù)發(fā)展［7，9］。土壤質(zhì)量的評(píng)價(jià)受研究尺度、土壤類(lèi)型、氣候等因素的綜合影響，選擇的評(píng)價(jià)指標(biāo)方法是相對(duì)的而不是絕對(duì)的；而對(duì)于農(nóng)田土壤，土壤質(zhì)量的評(píng)價(jià)與生產(chǎn)力評(píng)價(jià)類(lèi)似；在小區(qū)尺度上，需要進(jìn)行田間試驗(yàn)來(lái)理解土壤質(zhì)量和作物生產(chǎn)之間的關(guān)系［2，7］。

耕層土壤質(zhì)量評(píng)價(jià)是東北風(fēng)沙土區(qū)合理耕層構(gòu)建的重要內(nèi)容。由于土壤質(zhì)量評(píng)價(jià)的困難性，對(duì)耕層土壤質(zhì)量定量評(píng)價(jià)的研究相對(duì)較少，特別是以定量方法評(píng)價(jià)耕作措施對(duì)東北風(fēng)沙土耕層土壤質(zhì)量的影響鮮見(jiàn)報(bào)道。本文針對(duì)東北風(fēng)沙土區(qū)農(nóng)田，基于土壤生產(chǎn)性能指標(biāo)基礎(chǔ)上，利用聚類(lèi)分析建立最小數(shù)據(jù)集，并采用土壤質(zhì)量指數(shù)法驗(yàn)證最小數(shù)據(jù)集指標(biāo)的合理性，計(jì)算土壤耕作指數(shù)。應(yīng)用耕作指數(shù)定量評(píng)價(jià)耕作措施對(duì)東北風(fēng)沙區(qū)農(nóng)田合理耕層構(gòu)建效果，為選擇合理的耕作措施提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)于2017年4～10月在黑龍江省大慶市杜爾伯特蒙古族自治縣第一糧種場(chǎng)進(jìn)行。該區(qū)屬溫帶大陸性半干旱季風(fēng)氣候，夏季高溫多雨，冬季寒冷干燥，四季分明，雨熱同季，年均降水量400 mm，年平均氣溫3.1℃，≥10℃活動(dòng)積溫2 600～2 800℃，無(wú)霜期145～150 d，冬春季風(fēng)大，氣候干燥，蒸發(fā)力大，春旱嚴(yán)重，為典型春玉米農(nóng)作區(qū)。試驗(yàn)區(qū)土壤肥力屬于中下等，0～20 cm耕層土壤中砂粒、粉粒和粘粒含量分別為71.36%、10.65%和17.99%；20～40 cm土壤中砂粒、粉粒和粘粒含量分別為66.96%、12.33 %和20.71%；40～50 cm土壤中砂粒、粉粒和粘粒含量分別為68.01%、11.10 %和20.89%，土壤為砂壤土質(zhì)地。種植制度為一年一熟制，種植方式為壟作。試驗(yàn)前土壤基本理化性質(zhì):有機(jī)質(zhì)為15.2 g·kg-1，全氮為1.02 g·kg-1，全磷為 0.99 g·kg-1，全鉀為 24.1 g·kg-1，堿解氮53.8 mg·kg-1，有效磷為15.9 mg·kg-1，速效鉀為121.21 mg·kg-1，pH值為7.81，耕層厚度17.0 cm，犁底層厚度10.5 cm。

1.2 試驗(yàn)設(shè)計(jì)

試驗(yàn)設(shè)施化肥（NPK）、化肥+秸稈還田（簡(jiǎn)稱(chēng)，秸稈還田，NPK+S）和化肥+有機(jī)肥（簡(jiǎn)稱(chēng)，有機(jī)肥，NPK+OM）3種施肥模式，5種耕作方式（表1）共15個(gè)處理，每個(gè)處理3次重復(fù)，每個(gè)小區(qū)面積84 m2，隨機(jī)排列。5月5日播種，施玉米專(zhuān)用肥400kg·hm-2作基肥，在拔節(jié)期追施尿素1次，施肥量150kg·hm-2，秸稈全量還田，有機(jī)肥8 t·hm-2，10月2日測(cè)產(chǎn)，采集耕層土壤樣品。玉米品種為吉農(nóng)302，種植密度為5.5萬(wàn)株·hm-2。

表1 不同耕作措施田間操作描述

1.3 土壤理化分析

測(cè)定土壤理化性狀指標(biāo)16項(xiàng)，其中pH值、堿解氮、有效磷、速效鉀、全氮、全磷、全鉀、有機(jī)質(zhì)、容重、陽(yáng)離子交換量（CEC）、團(tuán)聚體平均重量直徑（MWD）、含水量、田間持水量和孔隙度參照《土壤農(nóng)業(yè)化學(xué)分析方法》進(jìn)行各項(xiàng)指標(biāo)的分析測(cè)定［14］；土壤滲透速率采用環(huán)刀法、耕層硬度采用田間原位硬度計(jì)測(cè)定，玉米產(chǎn)量測(cè)量采用多點(diǎn)采樣測(cè)量平均值。

1.4 最小數(shù)據(jù)集建立［7］

聚類(lèi)分析是對(duì)評(píng)價(jià)指標(biāo)的諸多特性進(jìn)行分類(lèi)。在SPSS 22中運(yùn)用Pearson距離對(duì)評(píng)價(jià)指標(biāo)進(jìn)行R型聚類(lèi)分析。在相應(yīng)聚合水平之間，將評(píng)價(jià)指標(biāo)劃分為反映耕層構(gòu)建不同方面特征的若干分組。根據(jù)相關(guān)分析，每組中相關(guān)性顯著的指標(biāo)可相互代替，選擇易度量、重現(xiàn)性好的指標(biāo)，剔除冗余指標(biāo)，選擇具有代表性且相互獨(dú)立的指標(biāo)進(jìn)入最小數(shù)據(jù)集。

1.5 耕作指數(shù)計(jì)算［15］

耕作指數(shù)（TI）是對(duì)耕層土壤結(jié)構(gòu)評(píng)價(jià)指標(biāo)的集成，耕作指數(shù)越大，土壤生產(chǎn)力越強(qiáng)，作物產(chǎn)量越高，則耕層土壤結(jié)構(gòu)越趨于合理。由于評(píng)價(jià)指標(biāo)對(duì)作物的作用方式不同，評(píng)價(jià)指標(biāo)得分隸屬函數(shù)類(lèi)型不同，本研究采取S型、反S型和拋物線型三類(lèi)函數(shù)。對(duì)數(shù)據(jù)集指標(biāo)做主成分分析，提取評(píng)價(jià)指標(biāo)的公因子方差，各指標(biāo)公因子方差占公因子方差之和的比例即為評(píng)價(jià)指標(biāo)權(quán)重值，計(jì)算公式如下:

式中Ai是各個(gè)評(píng)價(jià)指標(biāo)權(quán)重，X′i是各個(gè)評(píng)價(jià)指標(biāo)隸屬度值。

1.6 耕作指數(shù)評(píng)價(jià)精度驗(yàn)證［16］

利用Nash有效系數(shù)（Ef）和相對(duì)偏差系數(shù)（Er）評(píng)價(jià)基于聚類(lèi)分析最小數(shù)據(jù)集的精確程度。計(jì)算公式為:

式中I0和為基于全量數(shù)據(jù)集計(jì)算得出的耕作指數(shù)值和耕作指數(shù)平均值，Ical為基于MDS計(jì)算得出的耕作指數(shù)值。

1.7 數(shù)據(jù)處理與分析

采用SPSS 22軟件和Excel 2010對(duì)數(shù)據(jù)進(jìn)行分析處理。

2 結(jié)果與分析

2.1 東北風(fēng)沙土耕層土壤屬性描述性統(tǒng)計(jì)分析

耕層土壤質(zhì)量一般采用土壤理化指標(biāo)反映，耕層土壤屬性特征差異明顯（表2）。本研究選取16項(xiàng)耕層結(jié)構(gòu)因子作為評(píng)價(jià)指標(biāo)，有效磷（6.53 mg·kg-1）和滲透速率（5.68 mm·min-1）為中敏感度指標(biāo)，變異系數(shù)分別為37.29%和39.17%；堿解氮65.99～107.30 mg·kg-1、速效鉀100.33～ 205.93 mg·kg-1、全氮1.03～1.67 g·kg-1、團(tuán)聚體平均重量直徑0.41～0.74 mm和耕層土壤硬度113.20～196.73kg·cm-2的變異系數(shù)接近不敏感水平（＜10%）；pH值、有機(jī)質(zhì)、全磷、全鉀、陽(yáng)離子交換量、容重、含水量、田間持水量和孔隙度的平均值分別為7.67、18.40 g·kg-1、0.52 g·kg-1、24.09 g·kg-1、15.01 cmol·kg-1、0.82 g·cm-3、13.80%、27.55%和40.01%，它們的變異系數(shù)＜10%，屬于相對(duì)穩(wěn)定指標(biāo)，是維護(hù)耕層土壤結(jié)構(gòu)穩(wěn)定的主要因素。K-S檢驗(yàn)結(jié)果表明，各項(xiàng)指標(biāo)均服從正態(tài)分布。

表2 東北風(fēng)沙土農(nóng)田耕層質(zhì)量評(píng)價(jià)土壤參數(shù)統(tǒng)計(jì)特征

2.2 土壤質(zhì)量評(píng)價(jià)最小數(shù)據(jù)集

運(yùn)用SPSS 22中的聚類(lèi)分析法中R型聚類(lèi)，對(duì)全量數(shù)據(jù)集16個(gè)評(píng)價(jià)指標(biāo)進(jìn)行分類(lèi)。圖1表明，在聚合水平為10，評(píng)價(jià)指標(biāo)劃分為4類(lèi)，第一類(lèi)為有效磷、速效鉀、有機(jī)質(zhì)、全鉀、堿解氮、全氮、pH值、陽(yáng)離子交換量、全磷、田間持水量和孔隙度表征土壤基礎(chǔ)地力。由表3可見(jiàn)，土壤有機(jī)質(zhì)與堿解氮（0.792**）、全氮（0.647**）、速效鉀（0.812**）和全鉀（0.774**）均呈極顯著相關(guān)，土壤有機(jī)質(zhì)是反映土壤質(zhì)量的最重要的化學(xué)參數(shù)之一，除了作為營(yíng)養(yǎng)來(lái)源，它可以改善土壤的結(jié)構(gòu)和提高生物活性，是評(píng)價(jià)風(fēng)沙土土壤質(zhì)量重要評(píng)價(jià)指標(biāo)，選擇有機(jī)質(zhì)進(jìn)入最小數(shù)據(jù)集（MDS）。有效磷是作物從土壤中獲取的主要養(yǎng)分資源之一，有效磷與全磷（0.668**）、田間持水量（0.703**）均呈極顯著相關(guān)，有效磷變異性大（37.29%），選擇有效磷進(jìn)入MDS。全氮與陽(yáng)離子交換量（0.674*）和孔隙度（0.567*）都顯著相關(guān)，全氮是評(píng)價(jià)土壤肥力必備指標(biāo)，全氮進(jìn)入MDS。pH值是衡量土壤酸堿度的指標(biāo)，在土壤質(zhì)量評(píng)價(jià)使用率達(dá)90%［3］，故pH值進(jìn)入MDS。第二類(lèi)為團(tuán)聚體平均重量直徑，表征耕層土壤基本結(jié)構(gòu)性狀，是反映土壤物理結(jié)構(gòu)的重要指標(biāo)，則團(tuán)聚體平均重量直徑進(jìn)入MDS。第三類(lèi)為含水量和土壤水分滲透速率，表征土壤保水性能，土壤含水量和水分滲透速率相關(guān)系數(shù)0.569*，含水量可以直接反映土壤濕潤(rùn)狀況并且測(cè)定快速、方便，選擇土壤含水量進(jìn)入MDS。土壤容重和硬度為第四類(lèi)，表征土壤緊實(shí)度，對(duì)于特定的土壤類(lèi)型，容重可用于監(jiān)測(cè)土壤的硬度，通常情況下，土壤容重越大，土壤硬度相應(yīng)越大，土壤容重和硬度相關(guān)系數(shù)0.615*，土壤容重的變化除了本身影響水分和氧氣的供應(yīng)外，還可影響許多其它性質(zhì)和過(guò)程，選擇土壤容重進(jìn)入MDS。最終確定有機(jī)質(zhì)、有效磷、全氮、pH值、團(tuán)聚體平均重量直徑、含水量和容重7項(xiàng)指標(biāo)進(jìn)入MDS。提取評(píng)價(jià)指標(biāo)的公因子方差，得到進(jìn)入MDS指標(biāo)權(quán)重分別為有機(jī)質(zhì)0.088、有效磷0.148、全氮0.141、pH值0.050、團(tuán)聚體平均重量直徑0.117、含水量0.245和容重0.212。

圖1 東北風(fēng)沙土農(nóng)田耕層質(zhì)量評(píng)價(jià)聚類(lèi)分析樹(shù)圖

表3 東北風(fēng)沙土農(nóng)田耕層評(píng)價(jià)指標(biāo)Person相關(guān)系數(shù)矩陣

2.3 最小數(shù)據(jù)集合理性驗(yàn)證

最小數(shù)據(jù)集評(píng)價(jià)指標(biāo)體系合理性驗(yàn)證是耕層土壤質(zhì)量評(píng)價(jià)的重要環(huán)節(jié)［17］。計(jì)算基于不同數(shù)據(jù)集的耕作指數(shù)，全量數(shù)據(jù)集耕作指數(shù)（TI-TDS）變化幅度為0.35～0.87，均值為0.57，基于聚類(lèi)分析最小數(shù)據(jù)集的耕作指數(shù)（TI-MDS）在0.26～0.85之間變化，均值為0.53，TI-TDS和TI-MDS變化幅度、均值和變異系數(shù)接近。將TI-TDS與TI-MDS進(jìn)行線性擬合分析（圖2），相關(guān)系數(shù)r為0.928。TI-TDS與TI-MDS的Nash有效系數(shù)Ef為0.761，偏差系數(shù)Er為0.079，平均相對(duì)誤差分別為0.044和0.048。有效系數(shù)（Ef）越接近于1，表示基于MDS計(jì)算的耕層土壤質(zhì)量指數(shù)與TI-TSD基準(zhǔn)值越接近，精度越高［3，17］。相對(duì)偏差系數(shù)（Er）越接近于 0，表明基于MDS計(jì)算的耕層土壤質(zhì)量指數(shù)相對(duì)于TI-TDS基準(zhǔn)值偏差越小，結(jié)果越精確［17］。玉米產(chǎn)量分別與TI-TDS和TI-MDS進(jìn)行線性擬合（圖3），從擬合效果來(lái)看，玉米產(chǎn)量與TI-TDS、TI-MDS均呈極顯著正相關(guān)，相關(guān)系數(shù)r分別為0.722**和0.767**，達(dá)到極顯著水平，故MDS適合替代全量數(shù)據(jù)集對(duì)風(fēng)沙土耕地耕層土壤質(zhì)量進(jìn)行評(píng)價(jià)。

圖2 基于不同數(shù)據(jù)集東北風(fēng)沙土耕作指數(shù)的相關(guān)性

圖3 東北風(fēng)沙土不同數(shù)據(jù)集耕作指數(shù)與玉米產(chǎn)量的相關(guān)性

2.4 耕作措施對(duì)東北風(fēng)沙土耕層土壤質(zhì)量的影響

基于風(fēng)沙土耕地聚類(lèi)分析得出的最小數(shù)據(jù)集的耕作指數(shù)，耕作指數(shù)（TI-MDS）取值范圍在0～1之間，其值越高，表明作物產(chǎn)量越高，耕層結(jié)構(gòu)越合理［3］。根據(jù)耕作指數(shù)與作物產(chǎn)量的對(duì)應(yīng)關(guān)系，人為把土壤質(zhì)量劃定分為五級(jí)，耕作指數(shù)0.80～1.0為Ⅰ級(jí)、0.60～0.80為Ⅱ級(jí)、0.40～0.60為Ⅲ級(jí)、0.20～0.40為Ⅳ級(jí)、0～0.20為V級(jí)。由圖4表明，施有機(jī)肥耕層土壤耕作指數(shù)為0.70～0.85，平均為0.76±0.03，達(dá)到Ⅱ級(jí)水平；秸稈還田耕層土壤耕作指數(shù)為0.26～0.51，平均為0.41±0.04，達(dá)到Ⅲ級(jí)水平；施用化肥耕層土壤耕作指數(shù)為0.26～0.52，平均為0.41±0.04，達(dá)到Ⅲ級(jí)水平，這表明風(fēng)沙土區(qū)化肥配施有機(jī)肥是構(gòu)建合理耕層的有效培肥方式，秸稈還田對(duì)合理耕層構(gòu)建效果不顯著。從耕作方式角度來(lái)看，不考慮施肥因素，把3種施肥相同耕作方式平均耕作指數(shù)，翻耕＞旋耕＞常規(guī)＞超深翻＞深翻，翻耕達(dá)到Ⅱ級(jí)，其他耕作為Ⅲ 級(jí)，可見(jiàn)，翻耕可以改善耕層結(jié)構(gòu)。從耕作措施分析，翻耕施用有機(jī)肥耕作指數(shù)最大，為0.85，達(dá)到I級(jí)，深翻打破犁底層秸稈還田和超深翻施化肥耕作措施土壤質(zhì)量嚴(yán)重下降，耕作指數(shù)為0.26，耕層土壤為Ⅳ級(jí)。除超深翻以外，其他4種耕作方式秸稈還田沒(méi)有增加耕作指數(shù)，深翻打破犁底層秸稈還田耕作指數(shù)下降40.0%，由Ⅳ級(jí)地降為V級(jí)，說(shuō)明秸稈直接還田不適用于風(fēng)沙土區(qū)玉米生產(chǎn)。

圖4 耕作措施對(duì)東北風(fēng)沙土耕作指數(shù)的影響

3 討論

3.1 土壤耕作指數(shù)評(píng)價(jià)指標(biāo)選擇

耕層土壤質(zhì)量評(píng)價(jià)一般是基于農(nóng)作物的生產(chǎn)力、土壤理化及力學(xué)性質(zhì)［15］。準(zhǔn)確提取適宜評(píng)價(jià)指標(biāo)是耕層土壤質(zhì)量評(píng)價(jià)的重要環(huán)節(jié)，而在評(píng)價(jià)指標(biāo)篩選中，除了利用原始變量作為評(píng)價(jià)指標(biāo)外，從大量土壤參數(shù)中篩選出相對(duì)獨(dú)立、影響土壤質(zhì)量的敏感性指標(biāo)建立最小數(shù)據(jù)集［15］。最小數(shù)據(jù)集作為一種評(píng)價(jià)指標(biāo)篩選方法，在國(guó)內(nèi)外土壤質(zhì)量評(píng)價(jià)及監(jiān)測(cè)中得到了廣泛的應(yīng)用［18-19］。不同的研究目的對(duì)土壤質(zhì)量評(píng)價(jià)時(shí)選取的指標(biāo)存在差異［7］。從全量數(shù)據(jù)集指標(biāo)使用頻率排序來(lái)看，使用頻率最高的前3位是土壤容重、黏粒、pH值，是土壤質(zhì)量評(píng)價(jià)中必須考慮的指標(biāo)，使用頻率達(dá)90%；全氮使用頻率為70%；團(tuán)聚體平均重量直徑、粉粒、砂粒、有機(jī)質(zhì)和孔隙度次之，均為60%；有效磷使用頻率為50%［15］。本研究全量數(shù)據(jù)集16項(xiàng)評(píng)價(jià)指標(biāo)中有土壤容重、pH值、全氮、團(tuán)聚體平均重量直徑、有機(jī)質(zhì)、孔隙度和有效磷7項(xiàng)進(jìn)入全量數(shù)據(jù)集評(píng)價(jià)指標(biāo)使用頻率前10位。最小數(shù)據(jù)集中土壤容重、pH值和團(tuán)聚體平均重量直徑三者使用率均為40%；有機(jī)質(zhì)、全氮和鈉吸附比次之，使用頻率為30%；使用頻率為20%的有粘粒、土壤有效含水量、有效磷和有機(jī)碳［15］。最小數(shù)據(jù)集評(píng)價(jià)指標(biāo)中容重、pH值、團(tuán)聚體平均重量直徑、有機(jī)質(zhì)、全氮、含水量和有效磷7項(xiàng)評(píng)價(jià)指標(biāo)全部進(jìn)入最小數(shù)據(jù)集評(píng)價(jià)指標(biāo)使用頻率前10位，與前人研究結(jié)果相似，表明本研究中全量數(shù)據(jù)集和最小數(shù)據(jù)集評(píng)價(jià)指標(biāo)體系均有較好代表性，適宜于耕層土壤質(zhì)量評(píng)價(jià)。

3.2 土壤耕作指數(shù)評(píng)價(jià)土壤質(zhì)量

應(yīng)用耕作指數(shù)評(píng)價(jià)土壤質(zhì)量在國(guó)內(nèi)外得到廣泛應(yīng)用，對(duì)于土壤質(zhì)量的評(píng)價(jià)方法多種多樣，但目前尚無(wú)統(tǒng)一的被廣泛使用的土壤質(zhì)量評(píng)價(jià)方法，原因在于絕大多數(shù)方法無(wú)法同時(shí)做到準(zhǔn)確和簡(jiǎn)單實(shí)用［20］。近年來(lái)，越來(lái)越多的學(xué)者將聚類(lèi)分析與主成分分析方法相結(jié)合運(yùn)用于土壤質(zhì)量評(píng)價(jià)，以各原始指標(biāo)為變量進(jìn)行聚類(lèi)分析篩選出主要評(píng)價(jià)指標(biāo)，再以各主成分得分作為評(píng)價(jià)新指標(biāo)權(quán)重［17，21-22］，這一方法不僅簡(jiǎn)化了原始變量，而且消除多指標(biāo)間的相互影響，使評(píng)價(jià)結(jié)果更為準(zhǔn)確、可靠［17］。本研究利用聚類(lèi)分析與主成分分析相結(jié)合方法計(jì)算耕作指數(shù)，結(jié)果表明TI-TDS與TI-MDS呈極顯著正相關(guān)，相關(guān)系數(shù)r為0.928（圖2），TI-TDS與TIMDS的Nash有效系數(shù)Ef為0.761，偏差系數(shù)Er為0.079，TI-MDS取值范圍為0.26～0.85，與玉米產(chǎn)量顯著正相關(guān)（r=0.767），可作為定量評(píng)價(jià)土壤質(zhì)量的指標(biāo)。有機(jī)質(zhì)、有效磷、全氮和pH值屬于土壤基礎(chǔ)理化性質(zhì)，團(tuán)聚體平均重量直徑、含水量和容重測(cè)定方法簡(jiǎn)單便捷，利用公式即可計(jì)算TIMDS，可作為評(píng)價(jià)不同耕作措施的簡(jiǎn)單有效的方法。可以通過(guò)計(jì)算TI評(píng)價(jià)與試驗(yàn)地氣候土壤情況相似的東北風(fēng)沙土玉米田土壤質(zhì)量，是否能在其它地區(qū)及種植制度下運(yùn)用仍需進(jìn)一步的試驗(yàn)驗(yàn)證。

4 結(jié)論

依據(jù)聚類(lèi)構(gòu)建的最小數(shù)據(jù)集計(jì)算土壤耕作指數(shù)，能夠準(zhǔn)確地反映耕層土壤質(zhì)量情況，耕作指數(shù)TI（0.26～0.85）與玉米產(chǎn)量呈顯著正相關(guān)（r=0.767），利用此法計(jì)算TI-MDS簡(jiǎn)單有效，可以用做評(píng)價(jià)東北風(fēng)沙土玉米田土壤質(zhì)量的指標(biāo)，能預(yù)測(cè)不同耕作措施下的玉米產(chǎn)量。不同耕作措施下，風(fēng)沙土區(qū)化肥配施有機(jī)肥是構(gòu)建合理耕層的有效培肥方式，秸稈還田對(duì)風(fēng)沙土合理耕層構(gòu)建效果不明顯。不考慮施肥因素，從耕作方式角度來(lái)看土壤耕作指數(shù)，翻耕＞旋耕＞常規(guī)＞超深翻＞深翻。把耕作方式與有機(jī)物料培肥融為一體，翻耕+有機(jī)肥耕作指數(shù)最大，為0.85，達(dá)到I級(jí)，深翻+秸稈還田和超深翻+化肥耕作措施土壤質(zhì)量嚴(yán)重下降，耕作指數(shù)均為0.260，土壤為Ⅳ級(jí)。除超深翻耕作外，其它4種耕作方式秸稈還田沒(méi)有增加耕作指數(shù)，深翻+秸稈還田耕作指數(shù)下降40%，由Ⅳ級(jí)地降為V級(jí)，說(shuō)明秸稈還田不適用于風(fēng)沙土區(qū)玉米生產(chǎn)。