楊旭初,葉會(huì)財(cái)，李大明，余喜初，柳開樓,胡秋萍，胡志華，谷子寒，胡惠文，周利軍，黃慶海，林小兵

(江西省紅壤研究所/國(guó)家紅壤改良工程技術(shù)研究中心,江西 南昌 331717)

玉米作為中國(guó)3大主產(chǎn)糧食之一，據(jù)《2015中國(guó)統(tǒng)計(jì)年鑒》顯示，中國(guó)玉米播種面積占農(nóng)作物播種總面積的22.9%，玉米產(chǎn)量占農(nóng)作物總產(chǎn)量的36.1%[1]。隨著畜牧業(yè)的發(fā)展，玉米用途主要由糧食向飼料原料轉(zhuǎn)變，據(jù)廖永松等[2]估計(jì)，到2020年，80%以上的玉米將用來(lái)做飼料，我國(guó)飼料玉米仍有缺口，因此玉米生產(chǎn)同樣關(guān)乎國(guó)家糧食安全。

施肥對(duì)土壤生產(chǎn)力、養(yǎng)分回收率、作物產(chǎn)量等方面具有主要貢獻(xiàn)作用。國(guó)內(nèi)外長(zhǎng)期定位試驗(yàn)表明，均衡施肥或有機(jī)無(wú)機(jī)肥料配施是維持作物高產(chǎn)、穩(wěn)產(chǎn)的有效措施，但是長(zhǎng)期施用化肥(尤其是過量)，土壤質(zhì)量及生產(chǎn)力出現(xiàn)下降[3]。據(jù)段英華等[4]研究表明，湖南紅壤長(zhǎng)期施用化肥，土壤pH值平均每年下降0.06個(gè)單位，在土壤pH值4.5～6.3范圍內(nèi)，pH值每降低1個(gè)單位，玉米氮肥回收率下降10.9%。

土壤肥力評(píng)價(jià)至今無(wú)統(tǒng)一的方法、指標(biāo)和模型，目前主要有專家打分法、Fuzzy綜合評(píng)判法、全量數(shù)據(jù)集法、最小數(shù)據(jù)集法等[5-6]。除專家打分是主觀評(píng)價(jià)外，其余都是采用模糊數(shù)學(xué)和相關(guān)分析的理論，由于不同研究者對(duì)不同指標(biāo)重要性的認(rèn)識(shí)不同，因此土壤肥力指標(biāo)的選擇及權(quán)重確定各有不同[7]。例如包耀賢等[8]選用土壤有機(jī)質(zhì)、速效氮、速效磷、速效鉀、粘粒、團(tuán)聚度、團(tuán)聚體穩(wěn)定率8個(gè)指標(biāo)作為最小數(shù)據(jù)集評(píng)價(jià)紅壤性水稻土綜合肥力。溫延臣等[9]選用土壤物理、化學(xué)、生物學(xué)共14個(gè)指標(biāo)評(píng)價(jià)了潮土小麥-玉米輪作區(qū)土壤綜合肥力。

我國(guó)紅壤面積達(dá)218萬(wàn)km2，資源豐富，約占全國(guó)總土地面積的21.8%，主要分布在江西、湖南等長(zhǎng)江以南的低山丘陵區(qū)。目前南方紅壤旱地的土壤肥力評(píng)價(jià)方法較少，且缺少對(duì)其評(píng)價(jià)結(jié)果的客觀性驗(yàn)證。因此，本研究選自1986年建立的江西省進(jìn)賢縣紅壤旱地玉米定位試驗(yàn)，基于模糊數(shù)學(xué)和主成分分析的方法評(píng)價(jià)長(zhǎng)期不同施肥模式下土壤綜合肥力，并結(jié)合產(chǎn)量及產(chǎn)量穩(wěn)定性進(jìn)行驗(yàn)證，以期為玉米施肥提供科學(xué)指導(dǎo)。

1 材料與方法

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

玉米長(zhǎng)期定位試驗(yàn)位于江西省進(jìn)賢縣江西省紅壤研究所內(nèi)(116°20′24″N、28°15′30″E)，為典型低丘地形(海拔25 m，坡度5°)，土壤類型為第四紀(jì)紅粘土發(fā)育的紅壤。中亞熱帶季風(fēng)氣候，年平均氣溫18.7 ℃，年降水量1 807.5 mm ，年日照時(shí)數(shù) 1 657.8 h，>10℃有效積溫3 511.1℃。試驗(yàn)前耕層土壤理化性狀：pH值6.0，有機(jī)質(zhì)16.2 g/kg，全氮0.9 g/kg，全磷14.2 g/kg，全鉀15.8 g/kg，堿解氮60.3 mg/kg，有效磷12.9 mg/kg，速效鉀102.5 mg/kg。種植制度為春玉米-秋玉米-冬閑，1986～2015年春秋季玉米品種為掖單13號(hào)。

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

長(zhǎng)期定位試驗(yàn)共設(shè)10個(gè)處理，分別為：不施肥(CK)、單施豬糞(M)、單施氮肥(N)、單施磷肥(P)、單施鉀肥(K)、施氮磷肥(NP)、施氮鉀肥(NK)、氮磷鉀平衡施肥(NPK)、施2倍氮磷鉀化肥(NPK2)、施氮磷鉀化肥+豬糞(NPKM)。3次重復(fù)，小區(qū)面積22.2 m2，種植密度79 200株/hm2，春、秋季玉米肥料用量相同(表1)。氮肥為尿素(N 46%)，磷肥為鈣鎂磷肥(P2O512%)，鉀肥為氯化鉀(K2O 60%)。磷肥、豬糞作基肥施用，鉀肥在播種后第7 d全部做追肥施用，氮肥在播種后第7和第14 d各追肥50%。

表1 不同處理施肥量 (kg/hm2)

1.3 土壤肥力綜合指數(shù)計(jì)算

本研究構(gòu)建土壤肥力綜合指數(shù)(soil nutrient index，SNI)作為土壤養(yǎng)分綜合評(píng)價(jià)依據(jù)，其值越高，土壤肥力越好，具體計(jì)算過程如下：

(1)

(2)

式中，Wj為第j個(gè)指標(biāo)的權(quán)重。

(3)

式中，Qj表示第j個(gè)指標(biāo)的隸屬度值。

要計(jì)算SNI，一是土壤指標(biāo)的選擇；二是指標(biāo)權(quán)重的確定，三是隸屬函數(shù)的建立。

1.3.1 土壤指標(biāo)的選擇

參考文獻(xiàn)[10]選取土壤pH值、有機(jī)質(zhì)含量、有效磷含量、速效鉀含量、陽(yáng)離子交換量進(jìn)行紅壤旱地土壤肥力評(píng)價(jià)。土壤pH值的測(cè)定采用電極法(水土比為2.5∶1)，土壤有機(jī)質(zhì)的測(cè)定采用水合重鉻酸鉀氧化-比色法，土壤有效磷的測(cè)定采用碳酸氫鈉浸提-鉬銻抗比色法，土壤速效鉀的測(cè)定采用醋酸銨浸提-火焰光度法，土壤陽(yáng)離子交換量的測(cè)定采用醋酸銨交換-火焰光度法[11]。

1.3.2 指標(biāo)權(quán)重的確定

利用主成分分析對(duì)上述土壤指標(biāo)進(jìn)行權(quán)重確定。

1.3.3 隸屬函數(shù)建立和隸屬度的計(jì)算

隸屬度函數(shù)實(shí)際上是評(píng)價(jià)指標(biāo)與作物生長(zhǎng)曲線之間的數(shù)學(xué)表達(dá)式，它可以將不同量綱的指標(biāo)轉(zhuǎn)化數(shù)值為0～1之間的無(wú)量綱值，單一指標(biāo)隸屬度越大，反映該指標(biāo)的肥力水平越高。土壤pH值采用梯型隸屬函數(shù)，按公式(4)計(jì)算隸屬度；土壤有機(jī)質(zhì)含量、有效磷含量、速效鉀含量、陽(yáng)離子交換量采用S型隸屬函數(shù)，按公式(5)計(jì)算隸屬度。

(4)

(5)

公式(4)和(5)中X為上述土壤指標(biāo)實(shí)際含量，X1、X2、X3、X4分別為各土壤指標(biāo)的下臨界值、上臨界值、最優(yōu)值下限、最優(yōu)值上限。其值參考文獻(xiàn)[12]及江西紅壤旱地實(shí)際情況來(lái)確定。

表2 土壤指標(biāo)的隸屬函數(shù)類型和拐點(diǎn)值

1.4 土壤綜合肥力驗(yàn)證

玉米產(chǎn)量受當(dāng)季氣候影響波動(dòng)較大，因此選取試驗(yàn)30年的平均年產(chǎn)量及其變異系數(shù)(CV)、穩(wěn)定系數(shù)(SYI)3個(gè)指標(biāo)來(lái)驗(yàn)證土壤肥力[13]。產(chǎn)量統(tǒng)計(jì)方法：于1986～2015年每季玉米收獲時(shí)，整個(gè)小區(qū)單收單曬，并按14%籽粒含水率折算產(chǎn)量，年產(chǎn)量為春、秋季產(chǎn)量之和；年產(chǎn)量變異系數(shù)、穩(wěn)定系數(shù)按公式(6)、(7)計(jì)算，標(biāo)準(zhǔn)差越大，穩(wěn)定系數(shù)越低，變異系數(shù)越大。

(6)

(7)

1.5 數(shù)據(jù)統(tǒng)計(jì)

數(shù)據(jù)采用SPSS 17.0進(jìn)行方差分析， Duncan法進(jìn)行顯著性檢驗(yàn)。

2 結(jié)果與分析

2.1 長(zhǎng)期施肥對(duì)紅壤旱地土壤養(yǎng)分及特性的影響

由表3可知，長(zhǎng)期不同施肥處理對(duì)土壤pH值及有機(jī)質(zhì)含量、有效磷含量、速效鉀含量、陽(yáng)離子交換量均有顯著影響(P<0.05)。土壤pH值、有效磷含量、陽(yáng)離子交換量以單施豬糞(M)最高，偏施氮肥(N)最低；有機(jī)質(zhì)含量以化肥配施豬糞(NPKM)處理最高，CK最低；速效鉀含量以施2倍氮磷鉀化肥(NPK2)最高，偏施氮磷肥(NP)最低。

表3長(zhǎng)期施肥對(duì)土壤養(yǎng)分含量及特性的影響

處理pH值有機(jī)質(zhì)(g/kg)有效磷(mg/kg)速效鉀(mg/kg)陽(yáng)離子交換量(cmol/kg)CK5 20c13 50c10 20c80 00ef9 71eM6 40a17 83a183 50a141 67d12 78aN4 40e14 67c9 80c78 33ef9 47eP5 34c14 73bc23 00bc103 33e10 33cdK5 26c13 73bc20 70bc268 33b10 17cdNP5 18c13 90bc22 90bc65 00f10 60cNK4 80d15 70c14 90bc206 67c9 89deNPK4 85d15 27bc16 70c185 00c9 59eNPK24 84d15 43b29 80b321 67a9 94deNPKM6 03b20 00a177 30a246 67b11 73b

注：同列不同小寫字母表示在0.05水平上差異顯著，下同。

2.2 長(zhǎng)期施肥紅壤旱地土壤肥力評(píng)價(jià)

從圖1可知，5個(gè)土壤指標(biāo)的平均隸屬度以氮磷鉀配施豬糞(NPKM)最高，為0.94；其次為單施豬糞(M)，為0.92；單施氮肥(N)最低，為0.41。從5個(gè)指標(biāo)來(lái)看，土壤有效磷平均隸屬度最高，為0.83；土壤速效鉀次之，為0.75；土壤pH值最低，為0.52。

將土壤pH值、有機(jī)質(zhì)含量、有效磷含量、速效鉀含量、陽(yáng)離子交換量數(shù)據(jù)進(jìn)行球形假設(shè)檢驗(yàn)，KMO值為0.73，Bartlett值為127.98，sig=0.00<0.01，且Person相關(guān)性分析(表4)表明除土壤速效鉀外，土壤pH值、有機(jī)質(zhì)含量、有效磷含量、陽(yáng)離子交換量之間均呈極顯著正相關(guān)(P<0.01)，說明5個(gè)指標(biāo)之間可以進(jìn)行主成分分析。前2個(gè)主因子的方差貢獻(xiàn)率為92.0%，說明前2個(gè)主因子可以代表原始數(shù)據(jù)92.0%的數(shù)據(jù)。土壤pH值和土壤養(yǎng)分指標(biāo)的權(quán)重結(jié)果見表5。

圖1 土壤指標(biāo)平均隸屬度雷達(dá)圖

表4 土壤養(yǎng)分及特性的相關(guān)性

注：**表示在0.01水平上顯著相關(guān)。

從圖2可知，試驗(yàn)30年后，玉米土壤肥力綜合指數(shù)(SNI)從高到低依次為NPKM>M>K>P>NPK2>NPK>NP>NK>CK>N。化肥配施豬糞(NPKM)土壤肥力綜合指數(shù)最高，為0.94；單施豬糞(M)次之，為0.92；單施氮肥(N)最差，為0.42。

表5 變量在主成分上的載荷矩陣及變量權(quán)重

圖2 各土壤指標(biāo)綜合肥力貢獻(xiàn)率雷達(dá)圖

2.3 紅壤旱地土壤肥力驗(yàn)證

從表6可知，除單施豬糞(M)處理外，其余各處理平均年產(chǎn)量、產(chǎn)量穩(wěn)定系數(shù)、產(chǎn)量變異系數(shù)與土壤肥力綜合指數(shù)相比，排名規(guī)律基本相同。長(zhǎng)年單施豬糞(M)，土壤肥力綜合指數(shù)排第2，而平均周年產(chǎn)量、產(chǎn)量穩(wěn)定系數(shù)、產(chǎn)量變異系數(shù)分別排第3、第6、第6。

表6 玉米平均周年產(chǎn)量描述性統(tǒng)計(jì)

3 討論與結(jié)論

土壤基礎(chǔ)肥力是作物可持續(xù)生產(chǎn)的物質(zhì)基礎(chǔ)，長(zhǎng)期不同施肥模式通過影響土壤肥力從而影響土壤生產(chǎn)力。南方紅壤具有酸、瘦、板、粘的特點(diǎn)，其中土壤酸化嚴(yán)重[14]。本文研究表明，與試驗(yàn)前相比，試驗(yàn)30年后，單施豬糞和氮磷鉀配施豬糞處理土壤pH值提高0.40、0.03個(gè)單位，其余處理下降0.74～1.60個(gè)單位，對(duì)照處理下降0.80個(gè)單位。說明酸沉降等自然因素是紅壤旱地酸化的主要原因之一，施用有機(jī)肥豬糞和氮磷鉀配施豬糞可以維持和提高土壤pH值，是阻控紅壤酸化的有效措施，不均衡施肥顯著加劇紅壤酸化。這與前人[15-16]的研究結(jié)果一致。

南方紅壤富含鋁鐵氧化物，對(duì)有效磷固定能力很強(qiáng)，使磷素由可溶態(tài)向難溶態(tài)轉(zhuǎn)變，大大降低了土壤磷素的有效性，造成紅壤普遍缺磷[17]。本研究表明，不同處理間土壤有效磷含量為9.80～177.30 mg/kg，長(zhǎng)期施用豬糞和氮磷鉀配施豬糞土壤有效磷含量較對(duì)照分別提高16.99、16.38倍，活化了土壤磷素，提高了土壤有效磷含量?？赡艿脑蚴鞘┯秘i糞降低了土壤交換性Al3+濃度，減輕了鋁鐵氧化物對(duì)磷素的固定，從而緩解了鋁毒。

長(zhǎng)期定位試驗(yàn)系統(tǒng)科學(xué)地反映了土壤肥力演變及肥效變化規(guī)律，是研究土壤的經(jīng)典試驗(yàn)[18]。土壤肥力評(píng)價(jià)關(guān)鍵在于土壤指標(biāo)的選擇、指標(biāo)權(quán)重和隸屬函數(shù)閥值的確定，本評(píng)價(jià)體系表明，采用主成分分析確定指標(biāo)權(quán)重的方法是可行的。本文研究表明，紅壤旱地土壤肥力綜合指數(shù)排名以氮磷鉀配施豬糞最高，施用豬糞次之，偏施氮肥最低。對(duì)比產(chǎn)量、產(chǎn)量穩(wěn)定系數(shù)、變異系數(shù)排名順序，除單施豬糞處理外，其余各處理土壤肥力驗(yàn)證結(jié)果均較好。施用豬糞處理土壤肥力綜合指數(shù)和平均年產(chǎn)量排名靠前，產(chǎn)量穩(wěn)定系數(shù)和變異系數(shù)排名靠后?？赡艿脑蚴秦i糞中有機(jī)質(zhì)含量高，養(yǎng)分釋放過程較慢，其主要功能在于培肥，玉米產(chǎn)量逐年緩慢增加，因此，產(chǎn)量穩(wěn)定系數(shù)較低、變異系數(shù)較高。

[1] 郭澤林, 趙旭. 山東省糧經(jīng)作物播種比例演變及結(jié)構(gòu)優(yōu)化研究-基于1986～2015統(tǒng)計(jì)數(shù)據(jù)[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃, 2017, 38(7):164-171.

[2] 廖永松,黃季焜.21世紀(jì)全國(guó)及九大流域片糧食需求預(yù)測(cè)分析[J].南水北調(diào)與水利科技,2004,2(1):29-32.

[3] 朱兆良,金繼運(yùn).保障我國(guó)糧食安全的肥料問題[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2013,19(2):259-273.

[4] 段英華,徐明崗,王伯仁,等.紅壤長(zhǎng)期不同施肥對(duì)玉米氮肥回收率的影響[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2010,16(5):1108-1113.

[5] 高洪軍,彭暢,張秀芝,等.長(zhǎng)期不同施肥對(duì)東北黑土區(qū)玉米產(chǎn)量穩(wěn)定性的影響[J].中國(guó)農(nóng)業(yè)科學(xué),2015，48(23):4790-4799.

[6] 葉回春,張世文,黃元仿,等.北京延慶盆地農(nóng)田表層土壤肥力評(píng)價(jià)及其空間變異[J].中國(guó)農(nóng)業(yè)科學(xué),2013,46(15):3151-3160.

[7] 葉回春,張世文,黃元仿,等.粗糙集理論在土壤肥力評(píng)價(jià)指標(biāo)權(quán)重確定中的應(yīng)用[J].中國(guó)農(nóng)業(yè)科學(xué),2014,47(4):710-717.

[8] 包耀賢,黃慶海,徐明崗,等.長(zhǎng)期不同施肥下紅壤性水稻土綜合肥力評(píng)價(jià)及其效應(yīng)[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2013,19(1):74-81.

[9] 溫延臣,李燕青,袁亮,等.長(zhǎng)期不同施肥制度土壤肥力特征綜合評(píng)價(jià)方法[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(7):91-99.

[10] 曹志洪.中國(guó)土壤質(zhì)量[M].北京：科學(xué)出版社,2008.

[11] 鮑士旦.土壤農(nóng)化分析(3版)[M].北京：中國(guó)農(nóng)業(yè)出版社,2000.

[12] 何園球,孫波.紅壤質(zhì)量演變與調(diào)控[M].北京：科學(xué)出版社,2008.

[13] 戴茨華,王勁松,徐紅,等.旱地紅壤長(zhǎng)期定位施肥養(yǎng)分下移特征及作物產(chǎn)量效應(yīng)研究[J].中國(guó)農(nóng)學(xué)通報(bào),2012,28(9):70-74.

[14] 顏雄,彭新華,張楊珠,等.長(zhǎng)期施肥對(duì)紅壤旱地玉米生物量及養(yǎng)分吸收的影響[J].水土保持學(xué)報(bào),2013,27(2):120-125.

[15] 趙其國(guó),黃國(guó)勤,馬艷芹.中國(guó)南方紅壤生態(tài)系統(tǒng)面臨的問題及對(duì)策[J].生態(tài)學(xué)報(bào),2013,33(24):7615-7622.

[16] 張祥,王典,姜存?zhèn)},等.生物炭對(duì)我國(guó)南方紅壤和黃棕壤理化性質(zhì)的影響[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2013,21(8):979-984.

[17] 李昌新.長(zhǎng)期施肥對(duì)紅壤旱地玉米生產(chǎn)力和土壤肥力的影響及其機(jī)制研究[D].南京：南京農(nóng)業(yè)大學(xué),2009.

[18] 馬力,楊林章,沈明星,等.基于長(zhǎng)期定位試驗(yàn)的典型稻麥輪作區(qū)作物產(chǎn)量穩(wěn)定性研究[J].農(nóng)業(yè)工程學(xué)報(bào),2011,27(4):117-124.