孫翠平，張英鵬，羅加法，仲子文，孫 明，李 彥，劉兆輝,5，井永蘋，薄錄吉

(1.山東省農(nóng)業(yè)科學(xué)院 農(nóng)業(yè)資源與環(huán)境研究所，山東 濟(jì)南 250100；2.農(nóng)業(yè)部黃淮海平原農(nóng)業(yè)環(huán)境重點(diǎn)實(shí)驗(yàn)室，山東 濟(jì)南 250100；3.農(nóng)業(yè)部廢棄物基質(zhì)化利用重點(diǎn)實(shí)驗(yàn)室，山東 濟(jì)南 250100；4.山東省農(nóng)業(yè)面源污染防控重點(diǎn)實(shí)驗(yàn)室，山東 濟(jì)南 250100；5.山東省農(nóng)業(yè)科學(xué)院，山東 濟(jì)南 250100；6.新西蘭農(nóng)業(yè)科學(xué)院 魯亞庫(kù)拉研究中心，漢密爾頓 3240)

磷是作物生長(zhǎng)發(fā)育所需的三大營(yíng)養(yǎng)元素之一，其土壤有效磷可以直接被農(nóng)作物吸收利用，對(duì)作物生長(zhǎng)及產(chǎn)量具有重要影響；而土壤有效磷的變化與土壤磷盈虧具有密切關(guān)系[1-3]。因此，土壤有效磷含量的變化以及磷平衡一直受到研究者的關(guān)注[4]。試驗(yàn)證明，通過合理施肥可增加土壤磷素肥力，從而增加土壤有效磷含量，提高作物產(chǎn)量；但作物產(chǎn)量和土壤有效磷之間不呈完全的線性關(guān)系，兩者之間存在一個(gè)臨界值(磷農(nóng)學(xué)閾值)，而且超過農(nóng)學(xué)閾值一定限度時(shí)，還會(huì)造成磷環(huán)境污染[5]。磷農(nóng)學(xué)閾值是不同土壤類型、不同作物磷肥用量和最佳經(jīng)濟(jì)產(chǎn)量的重要指導(dǎo)依據(jù)。因此，研究施肥與磷素盈虧及與土壤有效磷的關(guān)系，明確不同作物有效磷農(nóng)學(xué)閾值的大小，對(duì)保障糧食安全和農(nóng)田生態(tài)可持續(xù)發(fā)展具有重要意義。

展曉瑩等[6]通過分析21個(gè)監(jiān)測(cè)點(diǎn)(包括5種土壤類型)單位磷平衡的有效磷變化量(有效磷效率)發(fā)現(xiàn)，其中81%的土壤監(jiān)測(cè)點(diǎn)有效磷與磷素盈虧之間呈顯著正相關(guān)關(guān)系。還有研究表明，土壤有效磷對(duì)磷素盈虧的響應(yīng)特征受土壤類型、施肥模式、肥料種類、施肥量及輪作模式等多種因素的影響[7-9]。有研究者在哈爾濱黑土中通過不同施肥模式研究得出土壤每虧缺和盈余100 kg/hm2時(shí)Olsen-P分別下降1.36～3.35 mg/kg和上升4.80～16.04 mg/kg，且認(rèn)為此變化范圍與土壤較低的pH值有很大的關(guān)系[10]。楊軍等[11]研究發(fā)現(xiàn)，PK、NP、NPK處理中土壤每盈余100 kg/hm2磷土壤有效磷分別增加3.59，1.19，1.75 mg/kg。Singh等[12]在不同初始土壤有效磷含量下施用相同磷量發(fā)現(xiàn)，隨著種植年限的增加有效磷的變化趨勢(shì)明顯不同；可見，有效磷對(duì)磷盈虧的響應(yīng)關(guān)系還與有效磷本底值有關(guān)。土壤磷素盈余會(huì)增加土壤有效磷含量，進(jìn)而提高作物產(chǎn)量，但當(dāng)土壤有效磷含量達(dá)到農(nóng)學(xué)閾值時(shí)，作物產(chǎn)量不再提高，過多磷素的添加反而會(huì)造成環(huán)境的污染[13]。不同模式、不同作物和不同模擬方法下的土壤有效磷農(nóng)學(xué)閾值有顯著的差異。沈浦[14]通過研究表明，單季旱作、雙季旱作和水旱輪作模式下，小麥和玉米土壤有效磷的農(nóng)學(xué)閾值分別位于7.5～23.5，5.7～15.2 mg/kg。還有研究者利用線性、線性+平臺(tái)和米切爾西方程模擬農(nóng)學(xué)閾值，其中米切爾西方程模擬得到的小麥和玉米農(nóng)學(xué)閾值最大，分別為16.6，15.9 mg/kg[15]，農(nóng)學(xué)閾值大小還與使用模擬方程有關(guān)。

潮土在山東省總面積達(dá)466.58 萬hm2，占全省土壤面積的38.53%，是山東省最大的土壤種類；耕地面積也居全省耕地面積之首，達(dá)48.12%，主要集中分布在魯西北黃河沖積平原地區(qū)[16]。而潮土磷素具有全磷含量較高、有效磷含量較低的特點(diǎn)[8]，需要通過合理施肥達(dá)到培肥地力、提高作物產(chǎn)量的目的。因此，在山東潮土區(qū)不同施肥模式下研究長(zhǎng)期定位試驗(yàn)對(duì)磷素盈虧、有效磷對(duì)磷盈虧響應(yīng)及土壤有效磷農(nóng)學(xué)閾值大小的影響具有重要意義。

1 材料和方法

1.1 試驗(yàn)地概況

長(zhǎng)期施肥定位試驗(yàn)在山東省濟(jì)南市山東省農(nóng)業(yè)科學(xué)院試驗(yàn)農(nóng)場(chǎng)院內(nèi)(北緯36°40′；東經(jīng)117°00′)，海拔為27.5 m，屬暖溫帶大陸性季風(fēng)氣候，年均溫度為12～14 ℃，年降水量580～700 mm，年日照時(shí)數(shù)1 870.9 h，全年無霜期216.4 d左右。供試土壤為成土母質(zhì)近代河流沉積物，土壤中的黏土礦物主要以水云母、綠泥石成分為主。該定位試驗(yàn)于1982年小麥季開始，試驗(yàn)開始時(shí)耕層土壤(0～20 cm)基本性質(zhì)為有機(jī)質(zhì)5.71 g/kg，全氮0.47 g/kg，全磷1.28 g/kg，堿解氮15.19 mg/kg，速效磷5.90 mg/kg，速效鉀75.3 mg/kg，pH 值8.20。

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

采用小麥-玉米一年兩熟輪作制，試驗(yàn)設(shè)單施化肥和有機(jī)無機(jī)肥配施兩大主處理，其中有機(jī)無機(jī)肥配施處理的無機(jī)肥施用量與單施化肥處理相同；單施化肥副處理包括：不施肥(CK)、氮(N)、氮磷(NP)、氮鉀(NK)、磷鉀(PK)、氮磷鉀(NPK)、減量氮磷鉀(N15PK)、增量氮磷鉀(N25PK)，有機(jī)無機(jī)配施副處理包括：有機(jī)肥(CK+M)、氮和有機(jī)肥(N+M)、氮磷和有機(jī)肥(NP+M)、氮鉀和有機(jī)肥(NK+M)、磷鉀和有機(jī)肥(PK+M)、氮磷鉀和有機(jī)肥(NPK+M)、減量氮磷鉀和有機(jī)肥(N15PK+M)、增量氮磷鉀和有機(jī)肥(N25PK+M)；每個(gè)小區(qū)1 m2，處理重復(fù)3次，試驗(yàn)采取隨機(jī)區(qū)組設(shè)計(jì)。供試小麥、玉米品種及各處理肥料用量同參考文獻(xiàn)[17]，合每年化肥帶入磷量為131.0 kg/hm2，有機(jī)肥帶入磷量為33.5 kg/hm2。

1.3 土壤采集與分析

土壤樣品分析及測(cè)定方法見參考文獻(xiàn)[17-18]。由于管理措施不變，且鑒于微區(qū)試驗(yàn)可持續(xù)性和土壤養(yǎng)分變化不顯著性，1982-1990年間隔1 a、1991-2010年間隔2～3 a、2011-至今間隔5 a左右采集1次土樣。

1.4 計(jì)算與統(tǒng)計(jì)方法

當(dāng)季土壤表觀磷盈虧=每年施入土壤磷素總量-每年作物(籽粒+秸稈)吸磷量 ①

土壤累積磷盈虧=∑[當(dāng)季土壤表觀磷盈虧] ②

農(nóng)學(xué)閾值的計(jì)算方法[19]：Yr=Yi/Ym×100 ③，式中：Yr(%)為籽粒的相對(duì)產(chǎn)量；Yi為每年各處理的籽粒產(chǎn)量；Ym為每年各處理的最大籽粒產(chǎn)量。

作物相對(duì)產(chǎn)量對(duì)土壤有效磷的響應(yīng)關(guān)系通過Mitscherlich方程模擬，公式如下：Y=A×(1-e-bx) ④，其中，Y是預(yù)測(cè)的相對(duì)產(chǎn)量；A是最大的相對(duì)產(chǎn)量；b是產(chǎn)量對(duì)土壤有效磷的響應(yīng)系數(shù)；由方程模擬出的相對(duì)產(chǎn)量為最大值的90%時(shí)，土壤有效磷含量為農(nóng)學(xué)閾值。

公式①～②單位均為kg/hm2，利用Origin 8.6進(jìn)行統(tǒng)計(jì)分析并作圖，差異顯著性檢驗(yàn)采用單因素Duncan法。

2 結(jié)果與分析

2.1 長(zhǎng)期施肥下土壤磷素盈虧

2.1.1 土壤當(dāng)季表觀磷素盈虧 不同施肥模式下33 a當(dāng)季土壤表觀磷盈虧如圖1(單施化肥A，有機(jī)無機(jī)肥配施B，下同)。隨年限的增加，當(dāng)季土壤表觀磷盈虧總體呈現(xiàn)略微下降的趨勢(shì)。單施化肥處理中，不施磷肥的3個(gè)處理(CK、N、NK)當(dāng)季土壤表觀磷盈虧一直呈現(xiàn)虧缺狀態(tài)，虧缺量為23.27～29.20 kg/hm2，顯著低于其他處理；由于N和NK處理作物品質(zhì)和產(chǎn)量原因帶走的磷多于CK處理，表現(xiàn)出更大磷虧缺。而施磷肥的5個(gè)處理(NP、PK、NPK、N15PK、N25PK)當(dāng)季土壤表觀磷呈現(xiàn)盈余狀態(tài)，平均位于68.34～85.99 kg/hm2，且以PK處理磷盈余最大，多數(shù)年限顯著高于氮磷鉀施肥。有機(jī)無機(jī)肥配施處理中，CK+M、N+M和NK+M處理當(dāng)季土壤表觀磷仍然處于虧缺狀態(tài)，虧缺量為12.10～22.27 kg/hm2，顯著低于其他處理，說明作物帶走的磷量大于有機(jī)肥帶入的磷量。NP+M、PK+M、NPK+M、N15PK+M和N25PK+M處理中，當(dāng)季土壤表觀磷處于盈余狀態(tài)，盈余量位于88.84～107.64 kg/hm2；其中，以PK+M和NP+M處理最高。

圖中不同字母表示各處理間差異顯著(P<0.05)。圖2同。Different lowercase letters in the figure indicate significant differences among different treatments (P<0.05).The same as Fig.2.

2.1.2 土壤累積磷素盈虧 不同施肥模式下33 a土壤累積磷盈虧如圖2。對(duì)于單施化肥處理，不施磷肥3個(gè)處理(CK、N、NK)土壤累積磷一直處于虧損狀態(tài)，3個(gè)處理間無顯著性差異，而盈虧量顯著低于其他5個(gè)處理，且隨種植時(shí)間延長(zhǎng)虧缺量累積增大。施用磷肥的5個(gè)處理(NP、PK、NPK、N15PK和N25PK)土壤累積磷一直處于盈余狀態(tài)，且盈余量隨種植年限的增加而增加，其中以PK處理最高，種植10 a后盈余量顯著高于其他處理；2015年5個(gè)處理土壤累積磷盈余值分別為2 567.25，2 837.63，2 255.20，2 399.51，2 378.98 kg/hm2，平均每年盈余77.80，86.00，68.34，72.71，72.10 kg/hm2。有機(jī)無機(jī)肥配施處理中，由于有機(jī)肥的施入，CK+M、N+M和NK+M處理土壤累積磷虧損值較單施化肥處理小；其中CK+M處理土壤磷虧損量最少，種植25 a后虧缺量顯著低于N+M和NK+M處理。其余5個(gè)處理累積磷盈余量更大，其中以NP+M和PK+M處理盈余量最大，種植20 a后顯著高于其他3個(gè)處理；2015年5個(gè)處理(NP+M、PK+M、NPK+M、N15PK+M和N25PK+M)土壤累積磷盈值分別為3 513.98，3 552.18，3 125.41，2 931.57，3 087.31 kg/hm2。由圖2還可得出，施氮量的大小對(duì)土壤磷平衡值的影響不大。

圖2 不同施肥處理土壤累積磷盈虧Fig.2 Accumulated P balance of different fertilization treatments

2.2 長(zhǎng)期施肥下土壤有效磷變化對(duì)土壤磷素盈虧的響應(yīng)

在兩大主處理下，對(duì)潮土土壤有效磷變量與土壤累積磷盈虧量進(jìn)行相關(guān)分析，如圖3所示，土壤有效磷增量和土壤磷素盈虧量之間呈現(xiàn)極顯著正相關(guān)關(guān)系(P<0.01)。對(duì)單施化肥和有機(jī)無機(jī)肥配施處理，土壤每盈余100 kg/hm2磷，Olsen-P含量分別上升0.76，1.46 mg/kg，可見有機(jī)無機(jī)肥配施單位土壤盈余磷更能增加土壤有效磷含量。

**.相關(guān)關(guān)系達(dá)到0.01水平。圖4同。**. Significant correlation at P<0.01. The same as Fig.4

2.3 土壤有效磷農(nóng)學(xué)閾值研究

米切里西方程在確定土壤磷素農(nóng)學(xué)閾值的研究方法中應(yīng)用較為廣泛，針對(duì)山東潮土不同施肥模式下有效磷含量與作物產(chǎn)量的長(zhǎng)期定位試驗(yàn)數(shù)據(jù)，利用Mitscherlich方程模擬作物相對(duì)產(chǎn)量與土壤有效磷二者的關(guān)系。結(jié)果為：在有機(jī)無機(jī)肥配施處理中，長(zhǎng)期施肥下小麥和玉米作物相對(duì)產(chǎn)量和土壤有效磷之間響應(yīng)關(guān)系不顯著(圖4)；單施化肥處理中，根據(jù)模擬方程得出小麥、玉米有效磷農(nóng)學(xué)閾值分別為12.60，8.90 mg/kg，潮土區(qū)冬小麥有效磷農(nóng)學(xué)閾值高于夏玉米。

圖4 單施化肥下小麥及玉米相對(duì)產(chǎn)量對(duì)土壤有效磷的響應(yīng)Fig.4 Response of relative wheat and maize yield to soil Olsen-P under single inorganic fertilizer treatments

2.4 磷肥應(yīng)用的指導(dǎo)

如圖5所示，不同土壤性質(zhì)、輪作制度和肥料種類及施用量下，土壤有效磷本底值和增加量不同，土壤有效磷農(nóng)學(xué)閾值也不同?；谕寥烙行Я妆镜字岛土灼胶獾捻憫?yīng)關(guān)系以及土壤有效磷農(nóng)學(xué)閾值的大小來推算肥料用量，對(duì)磷肥應(yīng)用具有重要的指導(dǎo)作用。

通過對(duì)土壤有效磷農(nóng)學(xué)閾值的研究得出小麥、玉米有效磷農(nóng)學(xué)閾值分別為12.60，8.90 mg/kg，而本試驗(yàn)土壤本底有效磷為5.90 mg/kg，與研究得出的小麥和玉米農(nóng)學(xué)閾值分別差6.70，3.00 mg/kg。以NP、PK、NPK 3個(gè)處理為例，從以上研究可知，到2015年時(shí)土壤累積磷平均每年盈余77.80，86.00，68.34 kg/hm2；且單施化肥處理每盈余100 kg/hm2土壤磷平均增加0.76 mg/kg土壤有效磷。因此，達(dá)到小麥農(nóng)學(xué)閾值NP、PK、和NPK處理分別需要12，11，13 a，達(dá)到玉米農(nóng)學(xué)閾值3個(gè)處理分別需要6，5，6 a。

圖5 磷肥應(yīng)用指導(dǎo)模式圖Fig.5 Picture of phosphate fertilizer application guidance model

3 討論與結(jié)論

3.1 長(zhǎng)期施肥下土壤磷素盈虧

由于不同區(qū)域生態(tài)因子的差異以及施肥量、種植結(jié)構(gòu)、氣候特征等因素，不同施肥模式下土壤磷盈虧及有效磷的響應(yīng)特征不同。首先，由于施肥量的關(guān)系有機(jī)無機(jī)肥配施處理磷盈余量明顯大于單施化肥處理。不施化學(xué)磷肥處理磷盈虧均處于缺乏狀態(tài)；可見只施用現(xiàn)在用量的有機(jī)肥也不能維持磷素平衡，需補(bǔ)充磷肥或增加有機(jī)肥用量來保持土壤磷素供應(yīng)。因此，合理施用磷肥或有機(jī)肥是土壤磷素供應(yīng)的保障。與平衡施肥處理相比，偏施磷肥處理(NP、PK、NP+M和PK+M)土壤磷素累積更多，這是由于在相同施磷量下，作物產(chǎn)量小，從土壤中帶走磷量也小，這與很多研究結(jié)論一致[11]。磷盈虧隨種植年限呈波浪狀變化，種植27 a磷盈余均明顯下降，是天氣、有機(jī)肥質(zhì)量和作物品質(zhì)與產(chǎn)量等綜合作用的結(jié)果。施用磷肥處理中，土壤中磷盈余量較大，隨種植年限增加導(dǎo)致磷素長(zhǎng)期累積于土壤中，勢(shì)必會(huì)造成磷素的淋溶及流失，綜合考慮下有必要適當(dāng)減少磷肥施用量，或利用有機(jī)肥替代部分磷肥。

3.2 長(zhǎng)期施肥下土壤有效磷對(duì)土壤磷素盈虧響應(yīng)

土壤有效磷是土壤磷素養(yǎng)分的重要指標(biāo)，最能反映土壤的供磷水平。通過土壤有效磷含量對(duì)磷素盈虧的響應(yīng)程度可用來指導(dǎo)施肥用量和施肥種類，對(duì)指導(dǎo)作物生產(chǎn)、施肥以及評(píng)價(jià)農(nóng)業(yè)環(huán)境磷風(fēng)險(xiǎn)具有重要意義[20]。

本研究中，兩大主處理每盈余100 kg/hm2磷素土壤有效磷分別增加0.76，1.46 mg/kg。中國(guó)7個(gè)樣點(diǎn)調(diào)查每盈余100 kg/hm2磷，Olsen-P含量上升范圍為1.44～5.74 mg/kg[21]，相較之本研究中有效磷轉(zhuǎn)化效率較低。從圖3也可以看出，單施化肥處理累積磷盈虧量<500 kg/hm2和有機(jī)無機(jī)肥配施處理磷盈虧量<2 000 kg/hm2時(shí)，Olsen-P增加較快，隨著磷累積量進(jìn)一步增加，Olsen-P有所降低；這可能主要與長(zhǎng)期耕作下單一的輪作制度和土壤性質(zhì)、結(jié)構(gòu)變差有關(guān)，生物多樣性降低，導(dǎo)致有效磷的轉(zhuǎn)化率降低。土壤有效磷對(duì)磷盈虧的響應(yīng)關(guān)系在很大程度上受土壤性質(zhì)和氣候條件的影響[9]。其差異還與試驗(yàn)規(guī)模、種植制度和試驗(yàn)誤差有很大的關(guān)系。本試驗(yàn)中單施化肥處理有效磷轉(zhuǎn)化低于有機(jī)無機(jī)肥配施處理，與很多研究結(jié)論一致[21]。林誠(chéng)等[22]研究得出，等磷素盈虧量下，有機(jī)無機(jī)肥配施的有效磷響應(yīng)系數(shù)高于化肥。

3.3 土壤有效磷農(nóng)學(xué)閾值研究

生產(chǎn)實(shí)踐中，作物產(chǎn)量對(duì)土壤肥力的響應(yīng)有一定的閾值(臨界值)，其中作物產(chǎn)量對(duì)土壤有效磷響應(yīng)閾值更為明顯，是評(píng)價(jià)施肥合理性的重要指標(biāo)。由米切里西模擬方程得出，單施化肥處理中，在獲得最大相對(duì)產(chǎn)量的90%時(shí)，小麥和玉米的土壤有效磷農(nóng)學(xué)閾值分別為12.60，8.90 mg/kg。郭斗斗等[19]同樣利用米切里西方程模擬河南潮土區(qū)小麥-玉米輪作下的農(nóng)學(xué)閾值，以獲得最大相對(duì)產(chǎn)量的95%為依據(jù)得出小麥玉米Olsen-P農(nóng)學(xué)閾值分別為13.1，7.5 mg/kg，與本研究結(jié)論較相似。劉彥伶等[7]對(duì)云南貴州黃壤性水稻土的玉米農(nóng)學(xué)閾值的研究發(fā)現(xiàn)，玉米磷農(nóng)學(xué)閾值平均值為15.8 mg/kg。李渝等[23]研究發(fā)現(xiàn)，貴州黃壤旱地的玉米農(nóng)學(xué)閾值平均為22.4 mg/kg?？梢姡魑镛r(nóng)學(xué)閾值受作物類型、土壤類型以及氣候條件等諸多因素的影響。

3.4 磷肥應(yīng)用指導(dǎo)

本試驗(yàn)中，有機(jī)無機(jī)肥配施處理當(dāng)季和累積磷盈余量明顯大于單施化肥處理，不施磷肥處理需施用磷肥維持土壤磷素供應(yīng)，其他處理需控制施磷量。

兩大主處理每盈余100 kg/hm2磷素時(shí)土壤有效磷分別增加0.76，1.46 mg/kg，有效磷轉(zhuǎn)化率較低，可通過調(diào)理土壤結(jié)構(gòu)性質(zhì)或有機(jī)肥替代部分磷肥來提高Olsen-P轉(zhuǎn)化率。

由米切里西模擬方程得出，單施化肥模式下小麥和玉米的土壤有效磷農(nóng)學(xué)閾值分別為12.60，8.90 mg/kg，現(xiàn)試驗(yàn)至33 a，理論上和實(shí)際中有效磷均已超過農(nóng)學(xué)閾值，需適當(dāng)減少磷肥施用。當(dāng)土壤有效磷含量下降到農(nóng)學(xué)閾值時(shí)，根據(jù)作物帶走的磷量，此時(shí)保證施入的磷量與作物帶走的磷量相當(dāng)，以此來確定磷肥施用量。在農(nóng)業(yè)生產(chǎn)應(yīng)用中，還需要及時(shí)觀測(cè)作物的生長(zhǎng)及土壤有效磷含量的變化來及時(shí)調(diào)整磷肥施用量。

綜合長(zhǎng)期施肥情況下土壤磷盈虧量變化、農(nóng)學(xué)閾值和土壤本底磷來看，每年施磷量可控制在33.5～131.0 kg/hm2。