制種玉米連作恒量施磷對灌漠土與潮土中磷素利用的影響

趙蕓晨, 王治江, 孫曉娟, 付余業

(1.河西學院 農業與生物技術學院, 甘肅 張掖734000; 2.甘肅省河西走廊特色資源利用重點實驗室, 甘肅 張掖 734000; 3.甘肅張掖市中天農業科技有限公司, 甘肅 張掖 734000)

磷素在農業生產與環境安全中具有重要影響，磷肥當季利用率低，在土壤中移動性小，固定強烈[1]；磷素施用后的產量效應逐漸降低，農田土壤磷素環境風險逐漸增大[2-3]，磷素已經成為農業生產的重要限制因素，研究表明隨著土壤磷素大量積累，土壤磷素逐漸由虧轉盈[4-5]，70%～90%的磷成為難以被作物吸收利用的固定形態[6]。石灰性土壤12 a磷肥施用，無機態磷大量增加[6]。棕壤26 a化學磷肥施用，除十鈣磷以外的其他磷素均增加[7]。

河西走廊張掖是全國玉米制種大區，玉米種子產量占全國用種量70%以上，為當地經濟提供了巨大貢獻。高利益驅動下，使玉米制種地塊連作嚴重，且85%以上農田采用恒量施肥模式，勢必會對玉米種子生產造成相應的影響，因此研究制種玉米連作恒量施磷模式下土壤磷素動態變化、土壤—植物磷素平衡，以及植物磷素吸收利用效率問題對該區玉米種子生產及當地生態環境變化都具有重要意義。作者前期研究表明10 a制種玉米連作生產，土壤總磷含量增加，速效磷先增后降[8]，但恒量施肥模式，長期連作，不同土類土壤磷素遷移、分級的研究報道還未見，外源磷素進入后在土體后發生的運移、轉化、積累、吸收利用動態，及其變化對生態環境影響的研究也未見報道。基于此，本研究擬以當地2種代表性石灰性灌漠土與石灰性潮土為研究對象，玉米制種連作8 a，研究連作條件2種不同土類土壤磷素動態變化特征及磷素養分虧缺、吸收利用規律，以期闡明制種玉米連作土壤磷素動態變化情況，為合理施磷，提高磷肥利用率，高產優質玉米種子生產提供合理的理論依據。

1 材料與方法

1.1 試驗地概況

試驗地位于甘肅省張掖市烏江鎮(海拔1 400 m，東經106°33′—107°47′，北緯37°04′—38°10′，土壤類型灌漠土)與張掖市臨澤縣鴨暖鄉(海拔1 600 m，東經99°51′—99°57′，北緯38°57′—39°02′，土壤類型潮土)，2種土壤基本理化性狀：灌漠土，pH值8.4，有機質17.5 g/kg，全氮0.91 g/kg，全磷0.661 g/kg，堿解氮63 mg/kg，速效磷7.5 mg/kg，速效鉀160 mg/kg。潮土pH 8.7，有機質13.5 g/kg，全氮0.83 g/kg，全磷0.635 g/kg，堿解氮56 mg/kg，速效磷7.1 mg/kg，速效鉀167 mg/kg。

1.2 供試品種

供試玉米品種為浚單22(河南秋樂種業集團提供)，播前種子包衣，包衣劑代號Q23，種子成熟后收獲全株，籽粒與秸稈分別收獲，秸稈不還田，測定籽粒、秸稈產量及籽粒、秸稈全磷含量。

1.3 試驗設計

試驗采用隨機區組設計，連續定位恒量施肥模式，于2006年開始，小區面積30 m×20 m，3次重復，小區四周設置保護行，制種玉米種植模式為“行比加滿天星”，播前覆膜。種植期間只施尿素與磷二銨，其中二銨用量525 kg/hm2(磷二胺，P2O548%)，氮肥用量675 kg/hm2(尿素，N 46%)，播前以底肥形式施入全部磷肥與1/3氮肥，拔節期及孕穗期追施剩余2/3尿素，通過藥劑防治紅蜘蛛及其他常見病蟲害，常規水肥管理，灌漠土小區代號A，潮土區代號B。

1.4 分析項目及方法

土壤與植物樣品：分別在連作的第1，3，5，7 a地上部分分秸稈與籽粒分別取樣(因2014年鴨暖鄉試驗地受蟲害影響沒有取樣，于連作第8 a取樣)為植物樣品，測定全磷含量，以末施磷區為空白對照CK。全株采收后利用5點法，每小區取6個樣點混合成一個土壤樣品，風干后測定全磷、速效磷，并分析磷的活化系數，樣品名稱分別為A1,A3,A5,A7和B1,B3,B5,B8；在0—20,20—40,40—60 cm分層采樣，樣品名稱仍為A1,A3,A5,A7和B1,B3,B5,B8，分析無機與有機分級磷。

測定項目與方法：土壤全磷量采用碳酸鈉熔融—鉬銻抗比色法[9]。速效磷采用NaHCO3浸提—鉬銻抗比色法。土壤無機磷分級采用顧益初和蔣柏藩法[10]，用0.25 mol/L的NaHCO3溶液浸提土壤二鈣磷，用0.5 mol/L的NH4OAC(pH=4.2)溶液浸提土壤八鈣磷，用0.5 mol/L的NH4F溶液浸提土壤鋁磷，用0.1 mol/L的NaOH與0.1 mol/L的Na2CO3混合溶液浸提鐵磷，用0.3 mol/L檸檬酸鈉-Na2S2O4溶液浸提閉蓄態磷，用0.5 mol/L的H2SO4溶液浸提十鈣磷，然后在分光光度計上采用硫酸鉬銻比色法測定。土壤有機磷分級采用Bowman-Cole法[11]，用0.5 mol/L的NaHCO3浸提活性有機磷，1 mol/L的H2SO4與1 mol/L的NaOH浸提中等活性有機磷，0.5 mol/L的NaOH浸提中穩性有機磷，0.5 mol/L的NaOH浸提高穩性有機磷，硫酸鉬銻抗比色法測定。玉米籽粒與秸稈中全磷含量采用H2O2—H2SO4消煮—鉬銻抗比色法測定[9]。

1.5 數據處理與統計

采用Excel 2007，SPSS 10.0軟件對數據進行統計分析和處理，采用單因素方差分析對數據進行顯著性檢驗，采用Person法對所有數據進行相關性分析。

表觀磷平衡(%)(SPB)[12]=(輸入磷-

輸出磷)/輸出磷×100%

(1)

磷素活化系數(%)(PAC)[13]=(土壤速效磷

含量/土壤全磷含量)×100%

(2)

作物磷素移出量(PTA)[14]=(籽粒全磷量×

籽粒產量)+(秸稈全磷量×秸稈產量)

(3)

磷肥表觀利用率(%)(APUE)[15]=

(施肥區作物養分含量-不施肥作物養分含量)×100/施肥量

(4)

2 結果與分析

2.1 制種玉米連作恒量施肥模式土壤磷素構成特征

連作8 a，2種連作土壤全磷含量分別為719.82與705.24 mg/kg，較處理前(0.661與0.635 g/kg)增加，增加量分別為58.01與69.91 mg/kg，年平均增加量為8.28與8.75 mg/kg，增加率為8.76%～11.01%。灌漠土速效磷含量先增后降，潮土中速效磷含量持續增加，2種土類速效磷較處理前分別增加1.66～4.64與1.01～4.04 mg/kg。磷的活化系數(PAC)為1.274%～1.794%，均小于2%(表1)。

表1 制種玉米連作總磷與速效磷的變化

注：同行不同小寫字母表明差異顯著性(p<0.05); T-P代表總磷; Av-P代表速效性磷; PA代表磷活化系數。下同。

連作8 a，土壤磷素組成以無機磷為主，占全磷總量65.3%～70.2%，無機磷主要以鈣磷(Ca-P)為主，其次為鋁磷(Al-P)和鐵磷(Fe-P)，而閉蓄態磷(O-P)只占全磷很少部分。鈣磷中十鈣磷(Ca10-P)含量最高，其次為八鈣磷(Ca8-P)， 最后為二鈣磷(Ca2-P)(表2)。總有機磷占全磷含量的6.3%～11.4%，主要以中度活性有機磷(MLO-P)為主，其次為中穩性(MRO-P)和高穩性有機磷(HPO-P)，活性有機磷(LO-P)含量最低(表3)。隨連作年限增加、無機化肥連續施用，土壤不同磷素：Al-P,Fe-P,Ca8-P,Ca10-P,LO-P,MLO-P均增加(表2—3)。制種玉米連作8 a，灌漠土中T-P,Ca10-P,Ca2-P,Av-P,LO-P含量較潮土高，而Fe-P,Ca8-P,MLO-P,HPO-P與MRO-P含量較潮土低，O-P相差不大(表2—3)。

表2 研究區不同連作年限土壤無機磷在不同土層含量特征 mg/kg

注：表中不同小寫字母表明差異顯著性(p<0.05); Ca2-P代表二鈣磷; Ca8-P代表八鈣磷; Al-P代表鋁磷，Fe-P代表鐵磷; O-P代表閉蓄態磷; Ca10-P代表十鈣磷; T-IP代表總無機磷。下同。

表3 研究區制種玉米連作有機磷在不同土層的分部特征 mg/kg

注：LO-P代表活性有機磷; MLO-P代表中度活性有機磷； MRO-P代表中穩性有機磷; HRO-P代表高穩性有機磷; T-OP代表總有機磷。下同。

2.2 制種玉米連作無機磷動態變化特征

恒量施磷模式下，隨著施肥年限增加，無機磷總量增加，增加比率為14%～18%。無機磷分級鈣磷所占比重最大，占無機磷總量81.15%～88.16%，其次為鋁磷，占無機磷總量10.73%～16.20%，閉蓄態磷含量最小，只占0.36%～0.64%。鈣磷中Ca2-P，Ca8-P，Ca10-P分別占無機磷總量5.32～9.95%，15.50～23.52%，55.66～64.19%，Ca10-P所占比重最大。2種不同土類隨磷肥施用年限增加：除Ca2-P顯著減少外， Ca8-P，Ca10-P，Al-P，Fe-P，O-P均增加；Fe-P在連作第3 a降至最低；灌漠土Ca10-P在連作第5 a，潮土Ca10-P在連作第8 a達最大，Ca8-P動態變化。灌漠土Ca2-P下降速率較潮土大，Ca10-P增加量較潮土高，Ca8-P與Al-P增加量較潮土低。不同剖面土層由上至下，各分級無機磷均呈現遞減趨勢(表4)，但各形態無機磷在不同土層所占總無機磷比例不同，Ca-P在20—40 cm土層比例最高，其次為0—20 cm土層，40—60 cm土層最低，O-P在40—60 cm土層比例最高，Al-P在20—40 cm土層比例最低，Fe-P在0—20 cm 比例最高。

Ca-P中，Ca10-P在40—60 cm土層比例最高，Ca2-P在0—20 cm土層比例最高，Ca8-P在20—40 cm土層最高(表4)。無機磷增加主要以Ca8-P與Ca10-P增加，呈現Ca10-P,O-P底聚，Ca2-P,Fe-P,Al-P表聚在0—20 cm，Ca8-P中層(20—40 cm)聚集現象(表4)。

表4 不同無機磷在不同土層占總無機磷的比例 %

2.3 制種玉米連作有機磷組分動態變化特征

如表5所示，有機磷各組分以MLO-P含量最高，占總有機磷52.10%～79.49%，其次為HRO-P，占9.31%～38.01%，MRO-P和LO-P含量較低，分別占6.64%～10.31%和3.21%～8.18%，隨連作年限增加，2種不同土類，土壤總有機磷各組分均增加，其中灌漠土中增加幅度最大的是MLO-P，其次為MRO-P和LO-P，HRO-P有下降趨勢，LO-P,MLO-P與MRO-P持續增加。

而潮土中有機磷各組分含量均隨連作年限增加而持續增加。在不同剖面土層中，2種不同土類各不同有機磷組分均呈現表層含量最高—表聚現象，但不同組分在不同土層中所占比例不同，灌漠土LO-P在20—40 cm土層最高,HRO-P在0—20 cm土層最高，MLO-P，MRO-P在40—60 cm土層最高，而潮土LO-P，MRO-P在40—60 cm土層最高，MLO-P在20—40 cm土層最高，HRO-P在0—20 cm土層最高。

表5 各有機磷在不同土層占總有機磷的分布比例 %

2.4 制種玉米連作土壤磷素動態平衡變化

依據不同連作年限制種玉米籽粒與秸稈產量，籽粒與秸稈中磷素含量，磷肥施入量以及土壤中磷素含量，得出土壤磷素收支表觀利用率等特征(表6)。由表6可以看出，當前恒量施肥模式，外源磷肥進入土壤產生173.10～175.65 kg/(hm2·a)的磷肥盈余量。灌漠土在連作第5 a，潮土在連作第3 a植物磷素帶走量最大，分別達到79.05與78.15 kg/hm2，隨后下降；而灌漠土在連作第3 a，潮土在連作第5 a磷素吸收量達到最大，分別為147.60與145.2 kg/hm2，然后下降，磷素的表觀利用率只有2.08%～4.89%，磷素表觀平衡系數均超過200%以上。

表6 不同土類作物磷素的帶走、吸收及可及性

注：籽粒/秸稈磷帶走量=籽粒/秸稈產量×籽粒/秸稈全磷量，總磷=籽粒磷+秸稈磷，剩余磷=施入磷肥量-總磷，磷肥表觀利用率=(施肥區作物養分含量-不施肥作物養分含量)×100/施肥量。

長期單施化肥使土壤磷素的可及性在一定程度上增加，但其增加量有限，分別在灌漠土上連作第3 a，潮土上連作第5 a達到最高，隨后下降，作物磷素吸收量也呈現相同規律，先增后降，制種玉米種子生產，總體磷肥利用率極低，大量磷肥在土體中盈余。

2.5 制種玉米連作土壤磷素間相關性分析

土壤各形態磷的相關關系如表7—8所示。由表7可知，在恒量施磷模式，制種玉米連作8 a，土壤全磷與土壤速效磷、總無機磷，八鈣磷、十鈣磷、鋁磷均呈現顯著正相關關系，相關系數分別為0.958,0.798,0.875,0.896,0.766，土壤中磷素總量的增加與上述不同形態磷素增加關系密切，土壤無機磷總量與土壤中八鈣磷、十鈣磷關系密切。土壤速效磷與土壤中全磷、二鈣磷、鋁磷、八鈣磷4項指標關系密切，說明土壤中速效磷的增加與上述磷密切相關。

如表8所示，速效磷與活性有機磷和中度活性有機磷也顯著相關。有機磷總量與中活性有機磷、高穩性有機磷關系密切，中活性有機磷與全磷間關系密切。

表7 各無機磷素間的相關系數

注：*表示顯著相關； **表示極顯著相關。下同。

表8 各有機磷素間的相關關系

3 討 論

3.1 衡量施肥模式對連作制種玉米各形態磷庫的影響

外源磷進入土壤后被有機、無機組分吸持，并最終以不同形態賦存于土壤，表現出不同的移動性和生物活性[16]。在本研究中，恒量施磷肥，制種玉米連作8 a，除Ca2-P以外的其他分級磷，全磷、速效磷均不同程度增加，說明外源無機磷肥施用可使土壤磷素增加，這與裴瑞娜及劉建玲長期施用無機磷均能顯著增加土壤中Olsen-P含量的研究結論相符[17,6]，但不同分級磷在不同土層所占比例不同，說明種子連作生產過程，不同分級磷在土體中轉化、累積、移動與被利用的程度不同。其中活性較高的Ca2-P含量顯著減少，表明其在連作過程被持續利用，是玉米種子生產主要的有效磷源，特別生產力較高的灌漠土消耗比例較潮土高。本研究中，Ca8-P與Ca10-P的量呈現積累—消耗—再積累模式，其中Ca8-P層聚在20—40 cm，Ca10-P在40—60 cm底聚，說明Ca8-P潛在被利用程度較Ca10-P高，Ca10-P被利用程度低。Al-P，Fe-P，O-P穩步增加，Fe-P呈現表聚現象，有被潛在利用的可能性。Al-P隨連作年限由表聚變底聚，可利用性由大變小。前人研究表明，Al-P是石灰性土壤中一種相當有效磷源[18-20]，Ca8-P，Al-P與Fe-P是作物第2有效磷源，無機磷對玉米有效性表現為Ca2-P>Al-P>Ca8-P>Fe-P>Ca10-P[20]。O-P與Ca10-P被認為是土壤中非活性磷的賦存形態[21]。本研究中，從相關系數的大小可以看出，各種形態磷對總磷的貢獻為：Ca10-P>Ca8-P>Ca2-P>Al-P>Fe-P >O-P，各形態磷對速效磷的貢獻為Ca8-P>Ca2-P>Al-P>Fe-P>Ca10-P>O-P，這說明本研究石灰性土壤中可有效增加土壤磷含量，但其底聚趨勢使其潛在有效被利用程度低，Ca8-P，Al-P，Fe-P分布在植物根系較豐富的區域，潛在被利用的程度較大，為潛在有效磷源，但其活性不同，在土壤中移動性也不同，在不同土層的積累量之間存在差異。由無機磷、速效磷的增加比率可以看出，河西走廊石灰性土壤恒量施磷，速效磷轉化比率小，有部分磷轉化為潛在植物磷源，還有大部分磷轉化成難利用性無機磷，如Ca10-P，逐漸向土壤下層移動，造成土壤中磷的累積。O-P也是作物潛在磷源[22]，但在本研究中，閉蓄態磷在土壤剖面中含量最低，且呈現明顯底聚現象，潛在被轉化利用的可能性較小。由上可以看出，不同磷素增加，但主要以難溶性磷增加為主，且磷土壤中的移動性都較小，造成土壤磷空間分布不均，土體剖面磷素動態分布不平衡現象，從而影響作物根系對養分的吸收，使所施磷肥的空間有效性降低。

本研究中，有機磷各形態增加，且均呈現表層含量較高現象，可以作為潛在磷源。中活性與活性有機磷可以作為植物有效磷供應的指標之一。從相關系數可以看出，長期單施化肥能促進植物難以利用的高穩性有機磷向活性、中活性有機磷轉化，施肥處理有利于增加耕層活性、中活性有機磷的含量。

3.2 衡量施磷模式制種玉米連作磷素動態利用與磷肥施用關系

作物對磷素的吸收利用不僅取決于土壤磷的含量和形態，還取決于磷素在土壤中的運移，無機磷是作物吸收利用的主要磷源，其在土壤中的含量和形態影響作物對磷素的吸收與利用，而有機磷對作物磷素的吸收也具有一定的影響。土壤全磷的變化與許多因素有關，許多研究表明長期過量施肥會導致土壤全磷含量升高[23]，本研究中，磷的增加量與土壤全磷變化線性相關，每年在土體中盈余較多磷肥，外源輸入磷量遠大于產品磷輸出量，說明全磷增加與外源磷肥進入直接相關，外源磷肥施入導致磷素在土壤中的積累，這與孫寧科等[24]24 a田間定位研究結果相似。本研究不同分級磷不同土層增加程度不同，說明肥料磷進行了部分轉化與移動。施入土壤的無機磷肥，主要向植物難利用的Ca10-P,MRO-P,HRO-P轉化，這可能與二銨的效應有關，二銨施入土壤后，會引起土壤磷酸鹽及土壤原有鐵鋁礦物發生分解，釋放鐵鋁，減少對肥料磷的固定和限制磷肥遷移有關[25]，不過金亮等[25]認為短時間內肥料磷不會轉化形成結晶態的難溶性含磷礦物，仍以水溶態存在，可以得出水溶態磷肥短期肥料效應與長期肥料效應不同，這造成不同連作年限外源磷肥的轉化效應在土體中進行了不同的表征。另外，本研究中2種土壤類型均為典型石灰性土壤，理論上存在碳酸鈣對磷素固定現象，使鈣磷有效性降低。有研究表明，當土壤磷素負荷超過一定的臨界點，磷的釋放能力將迅速增加，產生磷素在土體中的垂直遷移[26-28,14]，且Eghball等[29]認為外源化肥在土壤中的遷移距離為1.1 m，一旦超過，將進入其他循環中，磷素向下累積，超過農作物的需求水平，下層土壤對遷移的磷素具有很強的吸附能力[29]。本研究中，雖然土壤速效磷含量< 57 mg/kg，各分級磷在土體中垂直遷移性較小，但外源磷肥殘余量較大，有部分成為結合態磷，有部分成為潛在磷源，且活性很低的Ca10-P與O-P在土體下層呈現底聚趨勢，磷素活化系數<2%，磷素表觀利用率<5%，被利用的肥料磷只占到施用磷肥的15%左右，且磷肥一次性以基肥形式施入。Djodjic等研究表明，即使土壤固磷能力末達到飽和，若土壤中存在大孔隙，水溶性磷與吸附磷也會向下移[30],且磷肥后效較高[31]，超過磷素利用的臨界點，將使土壤磷極度不平衡，造成植物對磷的吸收發生紊亂。說明本試驗制種玉米連作恒量施磷水平條件下，全磷很難轉化為速效磷，其有效性不高[29]，在土壤中形成殘留，并向下積累，超過農作物的需求水平，土壤中磷素負荷正在逐漸加重，形成土壤磷空間分布不均狀態，使磷肥的空間有效性降低，從而造成植物磷素吸收紊亂，極有可能進入其它循環。恒量施磷是制種區長期制種玉米生產過程中一種保護產地穩定的政策，但其磷肥施用量遠遠超過農田植物獲取經濟產量實際所需要磷素水平，具有極大的潛在環境威脅性，其模式不利于生產區土壤磷素的合理循環與利用，增加了磷素離開土體，向水體或生態環境中遷移的可能，必將使生產區面臨巨大面源污染，導致相應的生態環境隱患。同時，磷素過量施用又會引起作物早熟，誘發土壤缺鋅，缺鉬，有害重金屬如鎘元素的富集[8]，使土壤堿性加重，理化性狀惡化。在今后的制種生產過程中，本區制種田磷素施用應減量或基本停止施用，進一步研究如何有效轉化土壤中殘余磷素。

3.3 衡量施肥模式2種不同土類磷素動態變化特征

石灰性潮土與灌漠土上定位施用磷肥，無論是分級磷含量還是分級磷在不同土層變化趨勢間均存在一定的差異性，從產量水平可以看出，灌漠土生產力較潮土高，從磷素的分級積累看，Ca8-P在潮土中移動性較灌漠土中強，灌漠土Ca10-P達到最高年限較潮土中早，潮土中除HRO-P外其他有機磷均向下移動性增加，2種不同土類磷素的移動性特征及磷素含量高低可能與土壤本身的肥力特征、水文運動特性及與地下水相連有關。研究表明，土壤水分條件與土壤中各個形態磷的相對含量有關，隨著水分的增加和還原作用的加強，會對土壤磷的吸附發生改變[32-33]。灌漠土發自于漠土綱，水分條件較差，在后期灌耕過程中，肥力逐漸增加，而潮土多源發于較低濕洼且地下水位較高的地方，水分動態較灌漠土活躍，因此其磷素運移特性較灌漠土更明顯，磷素的累積量也較灌漠土中明顯，但水文動態特性增加了磷素向土體下部遷移且進入水體及其他生態環境的危險性，潮土中磷素的施用更應該關注磷素的累積問題。

4 結 論

制種玉米連作恒量施磷8 a土體含磷量增加。2種不同農業土壤，除Ca2-P下降外，其他總磷及各分級磷均不同程度增加。Ca10-P及MLO-P含量最高，增加量最大，O-P及LO-P相對含量最低。總磷轉化率低，磷素在土體中移動性較小，磷肥表觀利用率小于5%，磷素活化系數<2%,外源輸入磷肥以174.16 kg/(hm2·a)量殘余在土壤中，大部分磷以低溶解態形式聚積在土壤表層。但隨連作年限增加，表層土壤對磷素特別對鈣磷及有機磷的固持及轉化率下降，賦存容量下降，呈現底聚趨勢，對生態環境健康存在極大威脅，應減量或停止施磷。潮土磷肥施用需更加精細化。

致謝：本研究在實施過程中，受秋樂種業生產部經理魯銀國大力支持，在土樣采集與管理過程受中天種業的大力支持，在此表示衷心的感謝！

[1] 夏立忠,Roy Anderson.長期施用牛糞條件下草原土壤磷的等溫吸附與解吸動力學[J].土壤,2000,32(3):160-164.

[2] 劉建玲,李仁崗,廖文華.白菜—辣椒輪作中磷肥的產量效應及土壤磷積累研究[J].中國農業科學,2005,38(8),1616-1620.

[3] Shen Jianbo, Yuan Lixing, Zhang Junliang, et al. Phosphorus dynamics: From soil to plant[J]. Plant Physiology, 2011,156:997-1005.

[4] 魯如坤,時正元,施建平.我國南方6省農田養分平衡現狀評價和動態變化研究[J].中國農業科學,2000,33(2):63-67.

[5] 李書田,金繼運.中國不同區域農田養分輸入、輸出與平衡[J].中國農業科學,2011,44(20):4207-4229.

[6] 劉建玲,張福鎖.小麥—玉米輪作長期肥料定位試驗中土壤磷庫的變化(Ⅱ):土壤Olsen-P 及各種形態無機磷的動態變化[J].應用生態學報,2000,11(3):360-364.

[7] 韓曉日,馬玲玲,王曄青,等.長期定位施肥對壤無機磷形態及剖面分布的影響[J].水土保持學報,2007,21(4)51-55,144.

[8] 趙蕓晨,秦加海,肖占文,等.長期定點施肥對制種玉米土壤理化性狀及重金屬含量的影響[J].水土保持學報,2012,26(6):204-208.

[9] 魯如坤.土壤農業化學分析方法[M].北京:中國農業科技出版社,2000.

[10] 蔣柏藩,顧益初.石灰性土壤無機磷分級體系的研究[J].中國農業科學,1989,22(3):58-66.

[11] Bowman R A, Cole C V. An exploratory method for fractionation of organic phosphorus from grassland soils[J]. Soil Science, 1978,125(2):95-101.

[12] 黃欣欣,廖文華,劉建玲,等.長期秸稈還田對潮土土壤各形態磷的影響[J].土壤學報,2016,53(3):779-789.

[13] 陳新,梁成華,張恩平,等.長期定位施肥對蔬菜保護地土壤磷素空間分布的影響[J].中國農學通報,2005,21(12):209-213.

[14] 王月立,張翠翠,馬強,等.不同施肥處理對潮棕壤剖面磷素累積與分布的影響[J].土壤學報,2013,50(4):135-142.

[15] 宋永林,李小平.長期施肥對作物氮磷利用及土壤速效氮磷供應能力的影響[J].磷肥與復肥,2008,23(2):71-78.

[16] 鄭順安,鄭向群,張鐵亮,等.石灰性紫色土上磷肥老化過程中土壤Olsen-P的變化特征及其模擬[J].環境化學,2012,31(9):1454-1455.

[17] 裴瑞娜,楊生茂,徐明崗,等.長期施肥條件下黑壚土有效磷對磷盈虧的響應[J].中國農業科學,2010,43(19):4008-4015.

[18] 史瑞和,邱嘉璋,陳幫本.江蘇省幾種主要土壤磷素供應狀況和磷肥效果[J].土壤學報,1962(4):374-379.

[19] 顧益初,蔣柏藩,魯如坤.風化對土壤粒級中磷素形態轉化及其有效性的影響[J].土壤學報,1984,21(2):134-143.

[20] 馮固,楊茂秋,白燈莎.32P示蹤研究磷素在石灰性土壤中形態及有效性變化[J].土壤學報,1996,33(3):301-307.

[21] 李中陽,李菊梅,徐明崗,等.長期施用化肥有機肥下我國典型土壤無機磷的變化特征[J].土壤通報,2010,41(6):1434-1439.

[22] 蔣柏藩,沈仁芳.土壤無機磷分級的研究[J].土壤學進展,1990,18(1):1-8.

[23] 王艷玲,何園球,吳洪生,等.長期施肥下紅壤磷素積累的環境風險分析[J].土壤學報,2010,47(5):880-887.

[24] 孫寧科,李隆,索東讓,等.河西農田磷鉀養分平衡及肥料利用率長期定位研究[J].土壤,2013,45(6):1009-1016.

[25] 金亮,周健民,王火焰,等.石灰性土壤肥際磷酸二銨的轉化與肥料磷的遷移[J].磷肥與復肥,2008,23(5):14-18.

[26] Heckrath G, Brookes P C, Poulton P R, et al., Phosphorus leaching from soils containing different phosphorus concentration in the Broad balk experiment[J]. Journal of Environmental Quality, 1995,24(5):904-910.

[27] Hesketh N, Brookes P C. Development of an indicator for risk of phosphorus leaching[J]. Journal of Environmental Quality, 2000,29(1):105-110.

[28] Eghball B, Binford G D, Baltensperger D D. Phosphorus movement and adsorption in a soil receiving long-term manure and fertilizer application[J]. Journal of Environmental Quality, 1996,25(6):1339-1343.

[29] 李學敏,張勁苗.河北潮土磷素狀態的研究[J].土壤通報,1994,25(6):259-260.

[30] Djodjic F, Borling K, Bergstrom L, Phosphorus leaching in relation to soil type and soil phosphorus content[J]. Journal of Environmental Quality, 2004,33(2):678-684.

[31] 王立春.充分發揮磷肥后效作用是實現玉米節本增效的重要舉措[J].玉米科學,2004(S2):91-94.

[32] De Mello J W V, Barron V, Torrent J. Phosphorus and iron mobilization in flooded soils from Brazil[J]. Soil Science, 1998,163(2):122-132.

[33] Miller A J, Schuur E A G, Chadwick O A. Redox control of phosphorus pools in Hawaiian montane forest soils[J]. Geoderma, 2001,102(3):219-237.