晏 軍,朱建強(qiáng),馬 艷,蔣舜堯,吳啟俠
(1.長(zhǎng)江大學(xué)農(nóng)學(xué)院,湖北 荊州 434025;2. 濕地生態(tài)與農(nóng)業(yè)利用教育部工程研究中心,湖北 荊州 434025)
水稻種植和水產(chǎn)養(yǎng)殖都是農(nóng)業(yè)中的重要產(chǎn)業(yè),對(duì)維護(hù)我國(guó)糧食安全和促進(jìn)農(nóng)民增收具有舉足輕重的作用。我國(guó)南方平原湖區(qū)水土資源優(yōu)越,適宜水稻種植和水產(chǎn)養(yǎng)殖,但是近年來(lái)也面臨諸多矛盾。一方面,我國(guó)水資源緊缺,水資源時(shí)空分布極不均勻,而占全國(guó)總用水量45%~50%的水稻種植系統(tǒng),水浪費(fèi)現(xiàn)象卻相當(dāng)嚴(yán)重[1]。而稻田過(guò)量施用化肥可導(dǎo)致氮磷等化學(xué)元素的大量流失和稻田周邊水體的污染[2,3]。另一方面,我國(guó)池塘養(yǎng)殖已越來(lái)越多地采用高密度、高投餌的精養(yǎng)技術(shù)模式,70%~80%投喂飼料以溶解和顆粒物的形式排入水體中,而被魚(yú)類(lèi)食用飼料中的營(yíng)養(yǎng)物質(zhì)最終約有51%的氮和64%的磷會(huì)成為廢棄物[4]。在池塘換水期間被集中排入外部環(huán)境后,會(huì)進(jìn)一步加劇了周?chē)雍秃吹人虻母粻I(yíng)養(yǎng)化程度,也造成一定的水資源浪費(fèi)。因此,合理利用池塘養(yǎng)殖廢水來(lái)灌溉稻田,進(jìn)行農(nóng)業(yè)利用是一條有效的途徑。
許多學(xué)者就養(yǎng)殖廢水灌溉作物做了大量的研究,陳柏湘[5]、吳啟俠等[6]分析了不同水管理模式下稻田對(duì)養(yǎng)殖廢水的凈化效果,杜會(huì)英等[7-9]研究了牛場(chǎng)肥水灌溉對(duì)冬小麥-夏玉米產(chǎn)量,養(yǎng)分吸收利用及土壤氮素平衡,黃紅英等[10]研究表明豬糞沼液施用可替代化肥50%~75%,可獲得與化肥處理相當(dāng)?shù)漠a(chǎn)量,且在一定程度上提高水稻氮素利用率,章明奎等[11]調(diào)查研究表明長(zhǎng)期畜禽養(yǎng)殖污水灌溉可促進(jìn)氮、磷在土壤剖面中的垂直遷移,增加對(duì)地下水的污染風(fēng)險(xiǎn)。這些研究多集中畜禽養(yǎng)殖肥水灌溉對(duì)作物產(chǎn)量、土壤理化性質(zhì)、修復(fù)及其處理工藝上,且研究多集中在北方旱地作物上,針對(duì)我國(guó)南方平原湖區(qū),池塘養(yǎng)殖廢水灌溉對(duì)水稻產(chǎn)量、養(yǎng)分吸收及土壤養(yǎng)分殘留等方面研究較少。本試驗(yàn)探究了減量施肥條件下,池塘養(yǎng)殖廢水灌溉對(duì)水稻氮磷吸收量、不同生育期稻田養(yǎng)分供應(yīng)及收獲后稻田土壤剖面氮磷殘留量進(jìn)行了研究分析。旨在實(shí)現(xiàn)池塘養(yǎng)殖廢水水分及其養(yǎng)分的再利用,節(jié)水節(jié)肥,達(dá)到水資源綜合利用,降低水田的灌溉施肥成本,為我國(guó)南方濱湖地區(qū)池塘養(yǎng)殖和水稻生產(chǎn)的有機(jī)結(jié)合提供數(shù)據(jù)理論依據(jù)。
本試驗(yàn)于2015年在湖北省荊州市荊州區(qū)太湖農(nóng)場(chǎng)梅槐分場(chǎng)進(jìn)行,該地區(qū)屬于亞熱帶季風(fēng)氣候區(qū),年平均氣溫為16~16.4 ℃,年平均降水量為904~1 127 mm,是典型的漁農(nóng)混作平原湖區(qū),周?chē)芗植贾咎锱c魚(yú)塘。用于本試驗(yàn)的水稻田0.667 hm2,南北方向長(zhǎng)133 m,東西方向長(zhǎng)50 m,屬于典型的水稻田。試驗(yàn)前耕層(0~20 cm)土壤化學(xué)性狀為:全氮量為1.91 g/kg、堿解氮量為69.4 mg/kg、全磷量為0.34 g/kg、速效磷量為24.1 mg/kg、全鉀量為3.34 g/kg、速效鉀量為118.7 mg/kg。養(yǎng)殖池塘位于水稻田的西北方向,池塘總面積2 000 m2,養(yǎng)殖魚(yú)類(lèi)以四大家魚(yú)為主,養(yǎng)魚(yú)總量為1 000 kg左右,每天早上8點(diǎn)左右進(jìn)行投餌喂食,餌料以雞糞、豬糞、發(fā)酵后的小麥、雜草和農(nóng)作物秸稈為主,每次投喂量約20 kg。通過(guò)輸水管和抽水泵將魚(yú)塘和稻田聯(lián)系在一起,灌溉時(shí)將魚(yú)塘水通過(guò)水泵抽到稻田各個(gè)小區(qū)內(nèi)。供試水稻品種為荊楚優(yōu)148(晚稻),6月1日播種,6月30日人工移栽,種植的株行距為25 cm×30 cm。每蔸1株,移栽前稻田內(nèi)撒施基肥,7月15日追施分蘗肥,8月24日齊穗,10月22日收割。
在0.667 hm2的水稻田內(nèi)沿東西向劃分為12個(gè)320 m2(40 m×8 m)的小區(qū),設(shè)4個(gè)處理:常規(guī)施肥量和清水灌溉(對(duì)照,CK)、常規(guī)施肥量+池塘養(yǎng)殖廢水灌溉(T1)、80%常規(guī)施肥+池塘養(yǎng)殖廢水灌溉(T2)及60%常規(guī)施肥+池塘養(yǎng)殖廢水灌溉(T3),每個(gè)處理3個(gè)重復(fù)。各小區(qū)間用高50 cm PVC隔水板隔開(kāi),防止試驗(yàn)過(guò)程中小區(qū)間的串水。稻田灌溉與池塘換水時(shí)期相結(jié)合,每次換水量約為500 m3,廢水和清水的灌水定額均為950 m3/hm2,利用水表計(jì)量,灌溉誤差在5%以?xún)?nèi),本試驗(yàn)過(guò)程中共有6次灌水,分別在分蘗期和灌漿期各2次、拔節(jié)孕穗期和抽穗揚(yáng)花期各1次,整個(gè)生育期灌溉量約為5 700 m3/hm2。各處理具體的養(yǎng)分投入量如表1所示。
試驗(yàn)開(kāi)始前采集基礎(chǔ)土樣。于水稻分蘗期、拔節(jié)孕穗期、抽穗揚(yáng)花期、灌漿期和成熟期利用土鉆采集土壤耕層樣品,采樣時(shí)間于灌溉后一周內(nèi),采樣深度為20 cm,收獲后采集土壤樣品深度為60 cm,每層20 cm,各小區(qū)樣品按“五點(diǎn)取樣法”采集。水稻成熟后每個(gè)小區(qū)收獲10株,分籽粒和秸稈2部分稱(chēng)其干質(zhì)量,風(fēng)干、脫粒、稱(chēng)量干質(zhì)量。籽粒和秸稈風(fēng)干后粉碎分別測(cè)定全氮和全磷量。各小區(qū)按處理測(cè)定其實(shí)際產(chǎn)量,并進(jìn)行考種和分析比較。
表1 各處理的養(yǎng)分投入量 kg/hm2
土樣和水稻植株樣測(cè)定指標(biāo)及方法為:全氮量采用濃H2SO4-H2O2消煮凱氏定氮法,全磷量采用濃H2SO4-H2O2消煮釩鉬黃比色法[12]。
應(yīng)用DPS15.10高級(jí)版進(jìn)行方差,采用LSD法進(jìn)行處理間多重比較,利用Microsoft Excel 2010作圖。
作物吸氮或磷量=(籽粒烘干質(zhì)量×籽粒含氮或磷量)+(秸稈干質(zhì)量×秸稈含氮或磷量)
收割后各處理水稻產(chǎn)量及其構(gòu)成見(jiàn)表2,由于各處理施肥水平不同,灌溉池塘養(yǎng)殖廢水后水稻產(chǎn)量也不盡相同。水稻產(chǎn)量以T1處理最高,為8 211 kg/hm2,相對(duì)CK處理可增產(chǎn)2.3%;以T3處理最低,為7 261 kg/hm2,相對(duì)CK處理減產(chǎn)9.6%,而T2處理其產(chǎn)量可獲得CK處理的99.4%。對(duì)各處理水稻產(chǎn)量進(jìn)行差異性比較,結(jié)果表明,各個(gè)處理間以及各處理與對(duì)照處理之間,有效穗和千粒質(zhì)量均為差異不顯著;每穗粒數(shù)、結(jié)實(shí)率和實(shí)際產(chǎn)量,T3與CK處理之間差異顯著,T2與CK處理之間相比實(shí)際產(chǎn)量有所下降,但差異不顯著。說(shuō)明池塘養(yǎng)殖污水灌溉只能在一定程度上可替代部分化學(xué)肥料,起到一定的增肥作用。
表2 池塘養(yǎng)殖廢水灌溉對(duì)水稻產(chǎn)量及構(gòu)成因素的影響
對(duì)收獲后水稻吸收氮磷量進(jìn)行顯著性分析,結(jié)果見(jiàn)圖1,氮
磷的吸收量以T1處理最高,相對(duì)于CK處理分別增加了4.4和2.0 kg/hm2,以T3處理最低,相對(duì)于CK處理降低了15.8和4.0 kg/hm2,T2處理相對(duì)于CK雖略有降低,但差異性并不顯著。表明在減量施肥的條件下,池塘養(yǎng)殖廢水中的營(yíng)養(yǎng)物質(zhì)可以補(bǔ)充水稻生長(zhǎng)所需的養(yǎng)分,80%常規(guī)施肥量條件下,用池塘養(yǎng)殖廢水進(jìn)行養(yǎng)分投入,對(duì)水稻吸收氮磷量的影響不顯著。
圖1 不同處理水稻的氮磷吸收量注:圖中不同字母表示在5%的水平上顯著差異(P<0.05),下同。
土壤全氮全磷是衡量土壤氮磷素供應(yīng)狀況的重要指標(biāo)。在水稻生長(zhǎng)的幾個(gè)關(guān)鍵生育期稻田耕層土壤氮磷養(yǎng)分狀況如圖3所示,可以看出,各處理耕層土壤氮磷含量隨生育期變化的趨勢(shì)基本一致,在分蘗期由于施肥和灌溉因素,其氮磷養(yǎng)分含量都有所提升,可見(jiàn)在水稻的營(yíng)養(yǎng)生長(zhǎng)期,其稻田養(yǎng)分供應(yīng)是相對(duì)充足的。在拔節(jié)孕穗期后,T2和T3處理全氮含量降低趨勢(shì)明顯,全磷含量除T3處理外其各處理降低趨勢(shì)相對(duì)較一致??傮w上看,各生育期稻田耕層土壤氮磷量大小表現(xiàn)為T(mén)1>CK>T2>T3,說(shuō)明灌溉池塘養(yǎng)殖廢水能夠提高稻田土壤氮磷量,但是在灌溉池塘養(yǎng)殖廢水帶入養(yǎng)分量相同的情況下,施肥量還是影響稻田土壤氮磷量的主要因素。從分蘗期至灌漿期總共計(jì)灌溉6次,在這個(gè)灌溉周期中,與CK相比,T1、T2和T3處理全氮含量分別平均增加3.8%、-2.3%和-10.1%,全磷含量分別平均增加3.3%、-4.2%和-10.7%,可見(jiàn)池塘養(yǎng)殖廢水灌溉條件下,水稻生長(zhǎng)期內(nèi)減量施肥20%稻田土壤全氮全磷含量沒(méi)有減量施肥40%降低趨勢(shì)大。
圖2 水稻不同生育期各處理耕層土壤氮磷量
水稻收獲后,不同土層全氮、堿解氮、全磷和速效磷含量如圖3所示??梢钥闯?,土壤全氮、堿解氮、全磷和速效磷含量在0~60 cm土層深度隨著土壤深度的增加而呈現(xiàn)降低的趨勢(shì),而在同一土層深度下,T1處理都高于CK處理,但會(huì)隨著施肥量的減少而減少,表明池塘養(yǎng)殖廢水灌溉配合常規(guī)施肥處理會(huì)增加土體中氮磷殘留量,說(shuō)明氮磷素供應(yīng)量超過(guò)作物需求后,將會(huì)導(dǎo)致更多的氮磷素殘留在土壤中,從而增加土壤中氮磷的累積量。各處理間差異性分析表明,土壤全氮(圖a)和堿解氮(圖b)含量分別在0~20 cm和40~60 cm土層各處理間無(wú)顯著差異,分別在20~60 cm和0~40 cm土層,T2和T3處理都顯著低于CK;土壤全磷(圖c)和速效磷(圖d)含量分別在0~20 cm和20~60 cm土層,T2和T3處理都顯著低于CK。同時(shí),由表3可知,水稻收獲后的實(shí)際產(chǎn)量、土壤全氮、堿解氮、全磷和速效磷都與氮磷的總輸入量呈顯著正相關(guān)。速效養(yǎng)分是能直接指示土壤供給植物吸收利用的,而全量養(yǎng)分可以用來(lái)衡量作物收獲后土壤中的殘留量。上述結(jié)果表明,減量施肥20%~40%條件下,池塘養(yǎng)殖廢水灌溉不會(huì)影響水稻對(duì)養(yǎng)分的吸收,同時(shí)也不會(huì)導(dǎo)致土壤剖面中氮磷的累積,從而可以降低氮磷對(duì)淺層地下水的環(huán)境風(fēng)險(xiǎn)。
圖3 不同處理土壤氮磷量的垂直分布
實(shí)際產(chǎn)量土壤全氮土壤全磷土壤速效磷土壤堿解氮總氮輸入量土壤全氮0.93*土壤全磷0.750.88*土壤速效磷0.870.810.89*土壤堿解氮0.89*0.98**0.96**0.88*總氮輸入量0.96*0.94*0.89*0.96**0.96*總磷輸入量0.95*0.99**0.90*0.89*0.99**0.97**
注:*p<0.05 ,**p<0.01。
氮磷是作物增產(chǎn)的主要限制因子,適量施入可顯著提高作物產(chǎn)量,過(guò)量的養(yǎng)分投入,降低養(yǎng)分利用率同時(shí)增加環(huán)境風(fēng)險(xiǎn)[8,13]。張興良等[14]研究表明豬場(chǎng)污水還田與化肥配施水稻、小麥產(chǎn)量與常規(guī)對(duì)照無(wú)顯著差異,杜會(huì)英等[8]研究表明麥季肥水灌溉能顯著提高冬小麥和夏玉米產(chǎn)量,但郭海剛等[15]研究表明灌溉牛場(chǎng)廢水次數(shù)過(guò)多會(huì)對(duì)冬小麥造成一定的負(fù)面影響,導(dǎo)致產(chǎn)量略有降低,本研究表明水稻產(chǎn)量以池塘養(yǎng)殖廢水灌溉常規(guī)施肥量下處理最高可達(dá)8 211 kg/hm2,相對(duì)常規(guī)處理可增產(chǎn)2.3%,當(dāng)80%常規(guī)施肥時(shí),其產(chǎn)量沒(méi)有降低,這與周元等[16]研究結(jié)果相一致。本研究中各生育期稻田耕層土壤氮磷量大小表現(xiàn)為T(mén)1>CK>T2>T3,說(shuō)明灌溉池塘養(yǎng)殖廢水能夠提高稻田土壤氮磷量,而抽穗期至灌漿期的水分養(yǎng)分管理對(duì)水稻產(chǎn)量形成的影響較大,池塘養(yǎng)殖廢水灌溉下減量施肥20%,抽穗揚(yáng)花期至灌漿期、灌漿期至成熟期稻田土壤養(yǎng)分沒(méi)有顯著下降,可見(jiàn)池塘養(yǎng)殖廢水中的營(yíng)養(yǎng)物質(zhì)能滿(mǎn)足后期水稻生長(zhǎng)的養(yǎng)分需求。同時(shí)有研究[17]表明,農(nóng)田土壤中養(yǎng)分的積累與不合理施肥有關(guān),其中因長(zhǎng)期高量施用化肥或畜禽糞便引起土壤磷、氮積累已成為許多農(nóng)區(qū)地表徑流中養(yǎng)分流失的主要原因。本研究表明,在水稻生長(zhǎng)季節(jié),池塘養(yǎng)殖廢水灌溉條件下施肥水平在常規(guī)施肥量的80%左右,能夠有效減少稻田土壤氮磷的積累,降低對(duì)農(nóng)田周?chē)w污染的風(fēng)險(xiǎn)。
稻田土壤氮磷流失主要通過(guò)地表徑流、侵蝕、淋溶(滲漏或亞表層徑流)和稻田排水進(jìn)入地表和地下水[18],這都與土壤中養(yǎng)分積累程度有關(guān)。相關(guān)研究指出[19,20],我國(guó)糧食作物的氮肥利用率只有20%~40%,當(dāng)季施用的磷至少有75% 以上以不同形態(tài)殘留在土壤中,因此,現(xiàn)代農(nóng)業(yè)對(duì)土壤肥力管理已不同于傳統(tǒng)農(nóng)業(yè),希望達(dá)到的土壤肥力既要能滿(mǎn)足農(nóng)作物正常生長(zhǎng)的需要,又不能使養(yǎng)分過(guò)度積累影響環(huán)境。有研究表明長(zhǎng)期[11]或過(guò)多[21]畜禽養(yǎng)殖污水直接灌溉可顯著地提高稻田表層土壤中氮和磷的積累,積累量隨灌溉年限的增長(zhǎng)而增加,會(huì)大大增加了對(duì)地下水污染的風(fēng)險(xiǎn)。本研究表明,土壤全氮和土壤全磷分別在20~60 cm和0~20 cm土層,T2和T3處理都顯著低于CK,說(shuō)明減量施肥20%~40%條件下,池塘養(yǎng)殖廢水灌溉不會(huì)導(dǎo)致土壤剖面中氮磷的累積,從而降低了氮磷的淋溶作用,這可能與作物品種、土壤類(lèi)型、土壤微生物、水體中營(yíng)養(yǎng)物質(zhì)的負(fù)荷和環(huán)境因素(如光照和溫度等)等因素有關(guān),且該研究主要側(cè)重于水稻生長(zhǎng)和土壤質(zhì)量的監(jiān)測(cè)評(píng)估,監(jiān)測(cè)時(shí)間僅為1 a。因此,研究長(zhǎng)期池塘養(yǎng)殖廢水還田的經(jīng)濟(jì)環(huán)境效益,還需要長(zhǎng)期試驗(yàn)進(jìn)一步深入研究確認(rèn)。
常規(guī)施肥量下池塘養(yǎng)殖廢水灌溉水稻增產(chǎn)潛能低,減量施肥20%后,在養(yǎng)分供應(yīng)減弱的情況下,保證水稻整個(gè)生長(zhǎng)期池塘養(yǎng)殖廢水灌溉,稻田土壤氮磷養(yǎng)分含量沒(méi)有顯著降低,且產(chǎn)量可達(dá)常規(guī)水平的99.4%。同時(shí),收獲后0~60 cm土層深度中氮磷量也是隨著施肥量的降低而降低。綜合考慮作物的產(chǎn)量效應(yīng)和土壤氮磷素累積的環(huán)境風(fēng)險(xiǎn),水稻全生育期池塘養(yǎng)殖廢水灌溉條件下,可節(jié)約化肥20%和清水灌溉量5 700 m3/hm2。
[1] 葉玉適, 梁新強(qiáng), 周柯錦, 等. 節(jié)水灌溉與控釋肥施用對(duì)太湖地區(qū)稻田土壤氮素滲漏流失的影響[J]. 環(huán)境科學(xué)學(xué)報(bào), 2015,35(1):270-279.
[2] Chirinda N, Cater M S, Albertb K R, et al. Emissions of nitrous oxide from arable organic and conventional cropping systems on two soil types[J]. Agric Ecosyst Environ, 2010,136:199-208.
[3] Li H, Liang X Q, Chen Y X, et al. Ammonia volatilization from urea in rice fields with zero-drainage water management[J]. Agric Water Manage, 2008,95:887-894.
[4] 劉長(zhǎng)發(fā), 綦志仁, 何 潔, 等. 環(huán)境友好的水產(chǎn)養(yǎng)殖業(yè)----零污水排放循環(huán)水產(chǎn)養(yǎng)殖系統(tǒng)[J]. 大連水產(chǎn)學(xué)院學(xué)報(bào), 2002,17(3):220-226.
[5] 陳柏湘. 稻田對(duì)池塘養(yǎng)殖廢水的凈化研究[D]. 武漢: 華中農(nóng)業(yè)大學(xué), 2009.
[6] 吳啟俠, 朱建強(qiáng), 李 谷, 等. 不同水管理模式下稻田對(duì)養(yǎng)殖肥水的凈化效果研究[J]. 灌溉排水學(xué)報(bào), 2011,30(4):5-9.
[7] 杜會(huì)英, 馮 潔, 張克強(qiáng), 等. 牛場(chǎng)肥水灌溉對(duì)冬小麥產(chǎn)量與氮利用效率及土壤硝態(tài)氮的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2016,22(2):536-541.
[8] 杜會(huì)英, 馮 潔, 郭海剛, 等. 麥季牛場(chǎng)肥水灌溉對(duì)冬小麥-夏玉米產(chǎn)量與磷吸收利用及土壤剖面分布的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2015,26(8):2 379-2 386.
[9] 杜會(huì)英, 馮 潔, 郭海剛, 等. 麥季牛場(chǎng)肥水灌溉對(duì)冬小麥-夏玉米輪作土壤氮素平衡的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào), 2015,31(3):159-165.
[10] 黃紅英, 曹金留, 常志州, 等. 豬糞沼液施用對(duì)稻、麥產(chǎn)量和氮磷吸收的影響[J]. 土壤, 2013,45(3):412-418.
[11] 章明奎, AHMED E, 鮑陳燕. 養(yǎng)殖污水直灌對(duì)稻田土壤氮、磷積累與垂直遷移的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2014,25(12):3 600-3 608.
[12] 鮑士旦. 土壤與農(nóng)業(yè)化學(xué)分析[M]. 3版. 北京: 中國(guó)農(nóng)業(yè)出版社, 2000.
[13] Ma W Q, Ma L, Li J H, et al. Phosphorus flows and use efficiencies in production and consumption of wheat, rice and maize in China[J]. Chemosphere, 2011,84:814-821.
[14] 張心良. 豬場(chǎng)污水還田與化肥配施對(duì)農(nóng)田水土環(huán)境和作物產(chǎn)量的影響[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2016,32(4):645-650.
[15] 郭海剛, 杜會(huì)英, 王 風(fēng), 等. 規(guī)?;?chǎng)廢水灌溉對(duì)土壤水分和冬小麥產(chǎn)量品質(zhì)的影響[J]. 生態(tài)環(huán)境學(xué)報(bào), 2012,21(8):1 498-1 502.
[16] 周 元, 朱建強(qiáng), 李 谷, 等. 稻田對(duì)池塘養(yǎng)殖肥水的吸收利用效果研究[J]. 灌溉排水學(xué)報(bào), 2011,30(5):78-81.
[17] Shafqat M N, Pierzynski G M. Long-term effects of tillage and manure applications on soil phosphorus fractions[J]. Communications in Soil Science and Plant Analysis, 2010,41:1 084-1 097.
[18] 呂家瓏, Fortune S, Brookes P C. 土壤磷淋溶狀況及其Olsen磷“突變點(diǎn)”研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2003,22(2):142-146.
[19] 王江立庚, 曹衛(wèi)星. 水稻高效利用氮素的生理機(jī)制及有效途徑[J]. 中國(guó)水稻科學(xué), 2002,16(3):261-264.
[20] Shafqat MN, Pierzynski G M. The effect of various sources and dose of phosphorus on residual soil test phosphorus in different soils[J]. Catena, 2013,105:21-28.
[21] 敏 鋒, 嚴(yán)正娟, 陳 碩, 等. 施用糞肥和沼液對(duì)設(shè)施菜田土壤磷素累積與遷移的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2016,35(7):1 351-1 359.