項(xiàng)大力,楊學(xué)云,孫本華,張樹(shù)蘭,古巧珍

(西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院,陜西楊凌712100)

磷素的面源污染已成為影響許多國(guó)家和地區(qū)的地表水質(zhì),導(dǎo)致水環(huán)境惡化和水體富營(yíng)養(yǎng)化的主要因素。資料表明,太湖流域農(nóng)田面源磷對(duì)水體磷的貢獻(xiàn)率高達(dá)19%[1],美國(guó)農(nóng)田面源排放的磷約占各種污染源總排放磷量的50%[2],水體污染問(wèn)題已成為環(huán)境研究的熱點(diǎn)[3]。農(nóng)田土壤中磷向地表和地下水體遷移的主要途徑包括地表徑流、土壤侵蝕和亞地表徑流(淋失或滲漏)。由于土壤有很強(qiáng)的磷固定能力,而磷肥主要施在耕層,含磷量很低的下層土壤是一個(gè)吸持磷素的巨大的容量庫(kù),所以一般認(rèn)為磷沿土壤剖面垂直向下淋失的可能性不大或淋失并不重要[4-6]。但隨著磷肥和有機(jī)肥的長(zhǎng)期大量施用,土壤磷素在耕層土壤大量累積甚至在一些局部地區(qū)接近飽和,最終可能發(fā)生較強(qiáng)的淋溶,特別是在大量施用有機(jī)肥時(shí)尤其如此[7-8]。

磷素淋失的形態(tài)一般為可溶性全磷(Total dissolved phosphorus,簡(jiǎn)寫(xiě)為T(mén)DP)和顆粒磷(Particulate phosphorus,簡(jiǎn)寫(xiě)為PP),可溶性磷包括鉬酸鹽反應(yīng)磷(Molybdate reactive phosphorus,簡(jiǎn)寫(xiě)為MRP)和可溶性有機(jī)磷(Dissolved organic phosphorus,簡(jiǎn)寫(xiě)為DOP)。有報(bào)道認(rèn)為,磷素淋失主要以可溶性形態(tài)為主[9-11],而另有報(bào)道則認(rèn)為以顆粒態(tài)磷為主[12-13]。磷素主要以何種形態(tài)遷移,取決于氣候因子、土壤因子、水文地理因子、農(nóng)作管理因子(如灌溉水平)等[11,14]。目前國(guó)內(nèi)外在施肥、土壤類型等對(duì)磷素淋失的影響方面做了很多研究,但對(duì)不同灌溉水平的影響報(bào)道極少。

大型滲漏池與各種各樣的土柱模擬試驗(yàn)相比更接近于田間實(shí)際情況;比較大田試驗(yàn)而言,更便于控制實(shí)驗(yàn)條件,也使得收集土壤溶液更加容易和便捷,因此是研究養(yǎng)分(溶質(zhì))運(yùn)移的一種非常有效的手段。本研究利用設(shè)在“國(guó)家黃土肥力與肥料效益監(jiān)測(cè)基地”的滲漏池群,研究土大量施用有機(jī)肥導(dǎo)致的土壤磷素快速累積條件下磷素能否淋出作物根區(qū)以及不同灌溉水平對(duì)土磷素淋失影響的程度和以何種形態(tài)為主,以期為該地區(qū)磷肥和灌溉水資源合理管理,有效防止農(nóng)業(yè)非點(diǎn)源污染提供依據(jù)。

1 材料與方法

1.1 試驗(yàn)地概況

試驗(yàn)設(shè)在“國(guó)家黃土肥力與肥料效益監(jiān)測(cè)基地”的滲漏池群上。該基地位于黃土高原南部的陜西省“楊凌國(guó)家農(nóng)業(yè)高新技術(shù)產(chǎn)業(yè)示范區(qū)”渭河三級(jí)階地,海拔524.7 m,年均氣溫13℃,年平均降水量550～600 mm。滲漏池建于1989年,為鋼筋混凝土結(jié)構(gòu),內(nèi)徑100 cm,深度120 cm,為擾動(dòng)回填土,按自然土壤發(fā)生層次和容重分層裝土,底部為石英砂顆粒,其上覆蓋一層石棉布,裝有PVC管(直徑為約2.5 cm)導(dǎo)流并用容器承接收集滲濾液,到本試驗(yàn)開(kāi)始已歷經(jīng)15年種植,可以認(rèn)為已基本接近自然土壤狀態(tài)。供試土壤屬褐土類土亞類,紅油土屬,厚層紅油土種,為粉砂粘壤土,黃土母質(zhì)。試驗(yàn)為冬小麥-夏玉米一年兩熟制。試驗(yàn)開(kāi)始時(shí)表層0—20 cm土壤有機(jī)質(zhì)10.92 g/kg,全氮 0.832 g/kg,全磷 0.61 g/kg,Olsen-P 4.2 mg/kg。自然土壤在77—155 cm處有粘化層,剖面性質(zhì)參見(jiàn)文獻(xiàn)[15]。

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

本試驗(yàn)始于2004年10月,設(shè)灌溉量為600(當(dāng)?shù)卮筇锷a(chǎn)常用灌水定額)、900、1200m3/hm23個(gè)處理水平;3個(gè)處理施肥相同,施肥量為N 180 kg/hm2(尿素,含 N 46%)、P2O5144 kg/hm2(過(guò)磷酸鈣,含P2O516%)、135 t/hm2有機(jī)肥(牛糞)。化肥每季施,有機(jī)肥每年在秋季小麥種植前一次施入;小麥季灌溉2～3次,一般在12月底冬灌一次,翌年3月春灌一次,或根據(jù)天氣在5月補(bǔ)充灌溉一次。玉米季灌溉3～4次,根據(jù)墑情而定。由于該深度只有8個(gè)滲漏池,除了 1200 m3/hm2重復(fù)2次外,其余均重復(fù) 3次。共種植7季作物,其中冬小麥4季。

1.3 滲濾液采集與分析

根據(jù)降水和灌溉不定期采集滲濾液,量其體積,并取400 mL水樣于當(dāng)天分析或4℃貯存。2004年10月至2008年5月底共采樣分析34次。

取一定量滲濾液,經(jīng)0.28 μ m微孔濾膜過(guò)濾,用鉬銻抗直接顯色分析[9]樣品磷含量,所得結(jié)果為鉬酸鹽反應(yīng)磷(MRP),一般多為無(wú)機(jī)磷。全磷(TP)測(cè)定用一定體積未經(jīng)過(guò)濾的滲濾液,采用改進(jìn)的Brookes等的方法[16],即:加入分析純高氯酸(70%～72%)2 mL和0.25 mL飽和氯化鎂溶液,高溫消化至近干,用0.6 mol/L稀鹽酸5mL加熱溶解,然后用鉬銻抗法顯色分析。取一定體積經(jīng)過(guò)0.28 μ m微孔濾膜過(guò)濾的滲濾液,采用與全磷相同的方法消化測(cè)定,所得結(jié)果為可溶性全磷(TDP)?？扇苄匀着c鉬酸鹽反應(yīng)磷的差值為可溶性有機(jī)磷(DOP)。全磷與可溶性全磷的差值為顆粒磷(PP)。

試驗(yàn)數(shù)據(jù)的統(tǒng)計(jì)分析用SPSS軟件的One-Way ANOVA進(jìn)行,1200m3/hm2處理由于只有2次重復(fù),按缺省值對(duì)待。

2 結(jié)果分析

2.1 不同灌溉水平對(duì)滲濾液量和累積淋失全磷的影響

不同灌溉條件下,各處理滲濾液和全磷累積量在不同采樣時(shí)間的變化(圖1)看出,累積滲濾液量隨灌溉水平的提高而升高,并且在試驗(yàn)初期差別較小,隨試驗(yàn)時(shí)間的延長(zhǎng)差異逐漸變大。至最后一次采樣(2008年5月底),灌溉量為900 m3/hm2滲濾液累積總量達(dá)193.8 mm,較600 m3/hm2處理(121.0 mm)處理增加了60.2%;而1200 m3/hm2處理則高達(dá)266.6 mm,分別比600和900 m3/hm2處理提高了120.3%和37.6%。累積淋失的全磷與滲濾液量有相似的趨勢(shì),也隨著灌溉量的增加而增加,但增加幅度要小。1200 m3/hm2處理在前23次取樣(前5季作物)和900 m3/hm2處理相差很小;900 m3/hm2處理累積淋出120 cm土層的滲濾液量自第23次以后才明顯的高于600 m3/hm2處理,而累積淋出的磷量則在第9次取樣時(shí)就明顯的高于后者。

圖1 不同灌溉處理累積滲濾液量與累積全磷量的變化Fig.1 Changes of cumulative leachate and TP of the irrigation treatments

2.2 不同灌溉水平對(duì)滲濾液中各形態(tài)磷濃度和累積量的影響

不同灌溉水平土體中鉬酸鹽反應(yīng)磷(MRP)、可溶性全磷(TDP)、顆粒磷(PP)和可溶性有機(jī)磷(DOP)等都有不同程度的淋移(表 1)。灌溉量為 1200 m3/hm2處理滲濾液中各個(gè)形態(tài)磷的濃度低于600和900 m3/hm2處理;900 m3/hm2處理的變幅略高于600 m3/hm2處理;各處理MRP的變幅總體上高于DOP和PP,后兩者較接近(表1)。

不同灌溉處理各形態(tài)磷累積量在不同采樣次數(shù)(時(shí)間)時(shí)的變化趨勢(shì)(圖2)看出,各形態(tài)磷的淋失量在試驗(yàn)開(kāi)始階段(第一季作物生長(zhǎng)期)處理間差異很小,后期則隨灌溉水平的提高而升高,即:灌溉量1200 m3/hm2處理＞900 m3/hm2處理＞600 m3/hm2處理。3個(gè)處理的MRP(圖2A)和PP(圖2C)在前20次沒(méi)有差異,而后1200和900 m3/hm2處理明顯高于600m3/hm2處理;而累積淋失的DOP在前8次(第一季作物生長(zhǎng)期)處理間無(wú)明顯區(qū)別,之后則隨灌溉量的增加而明顯升高(圖2D)。兩個(gè)高灌溉水平處理的TDP和DOP隨試驗(yàn)時(shí)間的延長(zhǎng)其累積淋失量差異較大,但MRP和PP的差異較小。

表1 滲濾液中各形態(tài)磷的濃度(mg/L)Table 1 Phosphorus concentrations of the fractions in the leachates

2.3 滲濾液中磷形態(tài)的比率及累積淋失量

為了更清楚地說(shuō)明土壤磷素淋失的主要形態(tài),表2列出了不同灌溉水平滲濾液中各形態(tài)磷的累積淋失量以及其與累積淋失的全磷或可溶性磷量的比率。淋失到120 cm以下土體的磷大部分為可溶態(tài)(TDP),約占66.4%～74.9%,且隨灌溉量增加而增加。在TDP中MRP和DOP約各占一半;MRP、DOP和PP分別占全部淋失磷總量的大概三分之一。灌溉量為900和1200 m3/hm2處理的滲濾液中MRP較600 m3/hm2處理高,而600 m3/hm2處理中PP占全磷量的比率高于其他兩個(gè)處理?？扇苄粤自谌字械姆蓊~有隨灌溉量增加而增加的趨勢(shì)。

圖2 不同灌溉處理滲濾液中各個(gè)形態(tài)磷累積量的變化Fig.2 Changes of cumulative amounts of the phosphorus fractions in the leachates under the irrigation treatments

表2 滲濾液中各形態(tài)磷累積淋失量及其在全磷(可溶性磷)中的比率Table 2 The amounts of the P forms and their ratios to TP or TDP in the lysimeter leachates

表2還看到,試驗(yàn)期間各形態(tài)磷的累積淋失量均隨灌溉量增加而增加,900和1200 m3/hm2處理的TDP和MRP顯著高于600 m3/hm2處理(P＜0.05),但900 m3/hm2處理的DOP和600 m3/hm2處理沒(méi)有顯著差異。3個(gè)處理的PP均沒(méi)有顯著差別。試驗(yàn)期間3個(gè)處理累積淋失的全磷量分別為83.4、164.0和209.4 g/hm2,900 m3/hm2處理較600 m3/hm2處理顯著增加(P＜0.05),而1200 m3/hm2處理則達(dá)到極顯著水平(P＜0.01)。

圖3看出,600 m3/hm2處理DOP、PP和MRP等3個(gè)形態(tài)磷累積淋失量很接近,DOP略高于PP和MRP(圖3A)。900 m3/hm2處理在第1季作物以后(第5次取樣)DOP高于PP和MRP;從第27次取樣后MRP高于PP,甚至在最后幾次取樣時(shí)間的累積量高于DOP(圖3B)。1200 m3/hm2處理也是在第1季作物種植后DOP開(kāi)始明顯高于其他形態(tài);自第24次取樣后MRP高于PP(圖3C)。

2.4 土壤剖面Olsen P的分布

圖4表明,經(jīng)過(guò)6季作物的種植施肥,2007年玉米收獲后土壤耕層(0—20 cm)Olsen-P的分布已高于70mg/kg,平均每年大約增加20 mg/kg;20—40 cm土層也增加了約6 mg/kg,但在0—100 cm土壤剖面中的各個(gè)層次,灌水處理間均沒(méi)有顯著差異。

3 討論

圖4 不同灌溉處理Olsen-P在土壤剖面的分布Fig.4 The distribution of Olsen-P in the soil profile under the irrigation treatments

水分運(yùn)動(dòng)和土壤磷素狀況是決定土壤磷向土壤深層移動(dòng)的最基本的兩個(gè)條件。本研究土壤滲濾液量和灌溉量相一致,但移出土體的磷總量卻并不完全與滲濾液的趨勢(shì)一致,尤其是900和1200 m3/hm2兩個(gè)處理在第5季作物之前(第23次取樣)全磷累積淋失量基本接近,之后兩處理間差異增大。其中原因可能與較厚的土層可以貯存較多的灌水,造成了試驗(yàn)前期3個(gè)處理淋出土體的滲濾液量沒(méi)有差異,因而磷總量也沒(méi)有差異(圖1);另外也可能是由于土壤磷含量在此之前相對(duì)較低,處于臨界點(diǎn)以下,盡管滲濾液量有差異,但并未明顯影響到磷移動(dòng)的總量。因?yàn)榈?季作物收獲后耕層土壤Olsen-P平均為68.6 mg/kg,第6季后為72.5 mg/kg,據(jù)此估計(jì)第5季作物收獲后大體上在70 mg/kg左右,這較有關(guān)報(bào)道的“轉(zhuǎn)折點(diǎn)”高出約10 mg/kg[9,17-18],這也可能是造成PP和MRP在前5季作物期間相對(duì)平緩的原因。在第5季作物之后,也同樣觀察到MRP和PP累積量的一個(gè)“躍遷”式增加,表明磷淋失臨界點(diǎn)的形成。

不同灌溉水平磷素淋失的形態(tài)都以可溶性為主,這與Heckrath等[9]在英國(guó)Broadbalk長(zhǎng)期肥料試驗(yàn)地排水管(埋深65 cm)中排水的結(jié)果一致;也與楊學(xué)云等[11]的原狀土柱模擬試驗(yàn)和張翠榮等[19]用盆栽試驗(yàn)所得結(jié)果也較一致。所不同的是,在模擬試驗(yàn)中,可溶性MRP隨取樣時(shí)間延長(zhǎng)升高,平均MRP占了淋失總量的三分之二[11],其原因可能是模擬試驗(yàn)土壤長(zhǎng)期處于高含水狀態(tài),水分與土壤作用時(shí)間長(zhǎng),還有相對(duì)處于較強(qiáng)的還原條件下,所以估計(jì)的結(jié)果可能會(huì)偏高。本試驗(yàn)以顆粒態(tài)磷淋失的比重也較模擬試驗(yàn)大,這可能是滲漏池是回填土壤,和原位土壤相比仍然存在一定的差異所致,也可能與模擬試驗(yàn)過(guò)高估計(jì)了可溶性磷降低了顆粒磷有關(guān)。

不同灌水處理都觀察到了磷素被淋到120 cm土體,盡管其淋失總量不到施入土壤總磷量的0.02%(表2),沒(méi)有任何農(nóng)學(xué)意義;在短期內(nèi)也不可能通過(guò)土壤剖面的化學(xué)分析(如土壤全磷和Olsen-P等)而感知其移動(dòng)和下層累積(圖4)。但是其平均濃度卻達(dá)到了0.03 mg/L以上,最高甚至達(dá)到了0.6 mg/L。而引發(fā)湖水富營(yíng)養(yǎng)化的磷(P)的臨界濃度只有0.02～0.035 mg/L[20-21]。實(shí)際上,滲漏池試驗(yàn)有可能低估了淋失量,因?yàn)楸驹囼?yàn)灌溉次數(shù)雖然較大田稍多,但一方面由于滲漏池本身的“邊際效應(yīng)”,使其土壤蒸發(fā)增加;另一方面由于土體在池的底部出現(xiàn)斷層和空氣接觸,而非自然土壤的連續(xù)體,從而失去了自然土體中下層土壤的毛管作用,水分有可能被空氣“托”在底部不能下移。因而從滲漏池淋出土體的滲濾液量就有可能低于自然條件下的淋失量。

盡管1200 m3/hm2處理的總淋失量高于其他兩個(gè)處理,但是其滲濾液中全磷和各形態(tài)磷的平均濃度都較其他兩個(gè)處理低,只有總滲濾液量高于600和900 m3/hm2處理,這可能是因?yàn)?200 m3/hm2處理的部分土壤水是通過(guò)優(yōu)先流途徑下移的。比較而言,水分和土體接觸時(shí)間較短,因此溶解的土壤(肥料)磷較低。

土壤磷的淋失或向下的淋移是一個(gè)緩慢的過(guò)程,肥料和水資源的管理對(duì)磷的淋失有極其重要的影響,因?yàn)樗墓芾碓谝欢ㄇ闆r下是不可控的,如強(qiáng)降雨情況。因而加強(qiáng)磷肥特別是有機(jī)肥的管理,是減少磷素淋失的重要手段,尤其是在地下水位較淺地區(qū),對(duì)于需水較多的作物,在肥水兩方面都要加強(qiáng)協(xié)調(diào)管理。

4 結(jié)論

[1]張維理,武淑霞,冀宏杰,等.中國(guó)農(nóng)業(yè)面源污染形勢(shì)估計(jì)及控制對(duì)策[J].中國(guó)農(nóng)業(yè)科學(xué),2004,37(7):1008-1017.Zhang W L,Wu S X,Ji H J et al.Estimation of agricultural non point source pollution in china and the alleviating strategies[J].Sci.Agric.Sin.,2004,37(7):1008-1017.

[2]Sharpley A N,Gburek W J,FolmarG J et al.Sources of phosphorus exported from an agricultural watershed in Pennsylvania[J].Agric.WaterManag.,1999,41:77-89.

[3]夏天翔,李文朝.撫仙湖北岸有機(jī)與常規(guī)種植菜地土壤氮、磷流失及累積特征[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2008,16(3):560-564.Xia T X,Li W C.Nitrogen and phosphorus loss and accumulation in organic and conventional vegetable fields in northern bank of Fuxian Lake[J].Chin.J.Eco-Agric.,2008,16(3):560-564.

[4]章明奎,王麗平.旱耕地土壤磷垂直遷移機(jī)理的研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2007,26(1):282-285.Zhang M K,Wang L P.Study on mechanisms of phosphorus downward transfer in arable soils[J].J.Agro-Environ.Sci.,2007,26(1):282-285.

[5]呂家瓏,Fortune S,Brookes P C.土壤磷淋溶狀況及其Olsen磷“突變點(diǎn)”研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2003,22(2):142-146.Lü J L,Fortune S,Brookes P C.Research on phosphorus leaching from soil and its Olsen-P “shreshold volume” [J].J.Agro-Environ.Sci.,2003,22(2):142-146.

[6]楊學(xué)云,孫本華,古巧珍,等.長(zhǎng)期施肥磷素盈虧及其對(duì)土壤磷素狀況的影響[J].西北農(nóng)業(yè)學(xué)報(bào),2007,16(5):118-123.Yang X Y,Sun BH,Gu Q Z et al.Phosphorus balances and its effects on soil phosphorus status under a 12-year long-term fertilization[J].Acta Agric.Bor.-Occid.Sin.,2007,16(5):118-123.

[7]Sims J T,Simard R R,Joern B C.Phosphorus loss in agricultural drainage:Historical perspective and current research[J].J.Environ.Qual.,1998,27:277-293.

[8]JagerP C,Claassens A S.Long-term phosphorus desorption kinetics of an acid sand clay soil fromMpumalanga,South Africa[J].Commun.Soil Sci.Plant Anal.,2005,36:309-319.

[9]Heckrath G,Brookes P C,Poulton P R et al.Phosphorus leaching from containing different phosphorus concentrations in the Broadbalk experiment[J].J.Environ.Qual.,1995,24:904-910.

[10]Stamm C,Fluhler H,Gachter R et al.Preferental transport of phosphorus in drainaged grassland soils[J].Environ.Qual.,1998,27:515-522.

[11]楊學(xué)云,古巧珍,馬路軍,等.土磷素淋移的形態(tài)研究[J].土壤學(xué)報(bào),2005,42(5):792-798.Yang X Y,Gu Q Z,Ma L J et al.Forms of phosphorus leaching in Loessial soil[J].Acta Pedol.Sin.,2005,42(5):792-798.

[12]Haygarth P M,Chapman P J,Jarvis S C et al.Phosphorus budgets for two contrasting grassland farming systems in the UK[J].Soil Use Manag.,1998,14:160-167.

[13]Simard R R,Beauchemm S,Haygarth P M.Potential for preferential pathways of phosphorus transport[J].Environ.Qual.,2000,29:97-105.

[14]呂家瓏.農(nóng)田土壤磷素淋溶及其預(yù)測(cè)[J].生態(tài)學(xué)報(bào),2003,23(12):2689-2701.Lü J L.Phosphorus leaching from agricultural soils and its prediction[J].Acta Ecol.Sin.,2003,23(12):2689-2701.

[15]楊學(xué)云,同延安,張樹(shù)蘭,等.長(zhǎng)期定位施肥對(duì)土剖面硝態(tài)氮分布與累積的影響[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2002,8(增刊):42-46.Yang X Y,Tong Y A,Zhang S L et al.Effect of long-termfertilization on nitrate nitrogen distribution and accumulation in the loess soil profile[J].Plant Nutr.Fert.Sci.,2002,8(Suppl.):42-46.

[16]Brookes P C,Powlson D S,JenkinsonD S.Measurement of microbial biomass phosphorus in soil[J].Soil Biol.Biochem.,1982,14:319-321.

[17]Tunney H,Carton O T,Brookes P C et al.Phosphorus loss from soil to water[J].CAB Inter.,1997,253-271.

[18]Heckrath N,Brookes P C.Development of an indicator for risk of phosphorus leaching[J].Environ.Qual.,2000,29:105-110.

[19]張翠榮,呂家瓏,潘楊,等.冬麥不同施磷水平下土壤磷素淋溶試驗(yàn)[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25(3):105-109.Zhang C R,Lü J L,P Y et al.Simulate of winter wheat phosphorus leaching in potted plant[J].Agric.Res.Arid Areas,2007,25(3):105-109.

[20]OECD.Eutrophication of water[R].Monitoring assessment and control.1982.156.

[21]Vollenweider R A.Scientific fundamentals of the eutrophication of lakes and flowing waters,with particular reference to nitrogen and phosphorus as factors in eutrophication[R].Tech.Rep.DAS/CSI,1968.