稻麥兩熟農(nóng)田稻季養(yǎng)分徑流流失特征

郭智，肖敏，陳留根，鄭建初

江蘇省農(nóng)業(yè)科學院農(nóng)業(yè)資源與環(huán)境研究所，江蘇 南京 210014

氮、磷含量過高是導致水體富營養(yǎng)化的根本原因。據(jù)報道，農(nóng)田徑流中氮、磷流失是水網(wǎng)富營養(yǎng)化的重要面源污染之一[1-3]。同時，稻田是太湖地區(qū)河網(wǎng)的主要農(nóng)田，減少稻田徑流中的氮、磷排放，對減緩水網(wǎng)的富營養(yǎng)化、改善水質(zhì)具有重要意義[4-5]。然而，目前尚無稻田生產(chǎn)過程中氮、磷等養(yǎng)分流失的確切數(shù)據(jù)以及其對周邊水體富營養(yǎng)化貢獻率的相關(guān)報道。因此，農(nóng)田排放氮、磷的規(guī)律及其對水體富營養(yǎng)化貢獻是亟待解決的關(guān)鍵問題。利用人工控制的徑流小區(qū)試驗是定量研究農(nóng)田養(yǎng)分流失特征的常用方法[6-8]。大量研究通過人工模擬降雨法探討了農(nóng)田氮、磷的輸出特征及機理[9-12]。然而，在自然降雨條件下通過定位試驗長期觀測對比研究秸稈還田、肥料運籌等不同處理方式下稻麥兩熟制農(nóng)田的養(yǎng)分輸出特征及其季節(jié)性變化的報道極少?；谌缟戏治?，本研究通過定性定量跟蹤分析稻麥兩熟制農(nóng)田養(yǎng)分在不同處理條件下地表徑流的流失規(guī)律，探討稻麥輪作模式下農(nóng)田養(yǎng)分排放對水體富營養(yǎng)化的貢獻，為水體富營養(yǎng)化人工治理提供一定理論依據(jù)。同時，為進一步提出有效控制養(yǎng)分流失策略提供技術(shù)支撐。

1 材料與方法

1.1 試驗材料及試驗點概況

試驗點設(shè)在江蘇省南京市溧水縣白馬鎮(zhèn)江蘇省農(nóng)業(yè)科學院植物科學基地，位于南京市溧水縣東南部，屬北亞熱帶向中亞熱帶的過渡區(qū)，年平均氣溫15.5 ℃，年平均日照2145.8 h，年平均降雨量1036.9 mm，年無霜期237 d。試驗時間為2009年05月至10月水稻生長季節(jié)。試驗田土壤屬黃棕壤，其基本理化性狀為：pH (6.21±0.53) (H2O, 1:5)，有機質(zhì)(16.62±3.15) g·kg-1，全 N(0.87±0.01) g·kg-1，全P(0.24±0.01) g·kg-1，速效氮(35.16±1.58) mg·kg-1，速效磷(11.84±2.23) mg·kg-1，速效鉀(89.23±3.84)mg·kg-1。水稻品種為南粳44號，5月15日播種，6月15日移栽，栽插規(guī)格為30 cm×13 cm，每穴2～3株。

1.2 試驗設(shè)施

選擇一地勢平坦田塊劃分為試驗小區(qū)，小區(qū)規(guī)格為15 m×4 m，為了防止小區(qū)間發(fā)生水分和養(yǎng)分的交換，每個試驗小區(qū)周圍垂直埋入70 cm深的雙層塑料膜將各個小區(qū)隔開。14個小區(qū)并排排列（兩側(cè)兩個小區(qū)為保護小區(qū)），一側(cè)設(shè)置灌水溝，每個小區(qū)靠近灌水溝一側(cè)有進水口，其高度保持一致，以保證各小區(qū)的灌水量保持一致。小區(qū)另一側(cè)為徑流收集池，其規(guī)格為4 m×1.5 m×1.5 m（長×寬×高），在靠近試驗小區(qū)一側(cè)距徑流池底1.0 m高度處開有Φ=15 cm的徑流收集管口，用PVC管與試驗小區(qū)相連，PVC管保持水平，以利于徑流水的收集，且每個徑流池上蓋有蓋板，防止雨水進入徑流池。

1.3 試驗處理

本試驗按照秸稈還田、肥料施用量及運籌三個因素設(shè)置 4個處理，分別為：常規(guī)施肥（T0），麥秸不還田，常規(guī)施肥量，基蘗肥︰穗肥=6︰4；秸稈還田（T1），秸稈全量還田，常規(guī)施肥量，基蘗肥︰穗肥=6︰4；還田減肥（T2），秸稈全量還田，減少 20%施肥量，基蘗肥︰穗肥=6︰4；肥料運籌（T3），秸稈不還田，常規(guī)施肥量，基蘗肥︰穗肥=4︰6。常規(guī)施肥量分別為 N 300 kg·hm-2，P 90 kg·hm-2和K 135 kg·hm-2，P肥全作基肥，K肥基肥施用50%，促花肥施50%；在試驗開始之前，各小區(qū)施用有機肥 3750 kg·hm-2。麥秸還田量為 7500 kg·hm-2。每個處理重復3次，共12個試驗小區(qū)，隨機排列。

1.4 徑流樣品的采集與分析

分別于降雨和烤田期間利用水泥徑流池收集徑流量，記錄稻田的排水量并采樣。每次采集徑流水樣時，先測量徑流水深，用以計算徑流水量。然后將徑流池內(nèi)的水攪拌均勻，采集徑流水樣500 ml，過濾分離土壤顆粒。

供試土壤的基本理化性質(zhì)：按照V(去離子水)︰m(土)=5︰1的比例浸提，用原位pH計（型號）測定pH[13]；全N、全P、速效K和速效P參照鮑士旦主編的土壤農(nóng)化分析中的方法測定[14]；有機質(zhì)采用重鉻酸鉀容量法-稀釋熱法[14]。分別采用紫外分光光度法、紫外消解鉬藍比色法和火焰光度法測定徑流水中全N、全P、速效K的質(zhì)量濃度[13]。

1.5 數(shù)據(jù)處理

用Microsoft Office Excel (2003)和SPSS for windows (11.5)軟件進行數(shù)據(jù)處理，文中所列數(shù)據(jù)均為3次重復平均值。

2 結(jié)果與分析

2.1 稻季徑流水量變化規(guī)律

本試驗過程中，共產(chǎn)生7次徑流，分別為6月22日、6月30日、7月7日、7月22日、7月28日、8月3日和8月11日，對應(yīng)7次降雨過程，其累計降雨量分別為17.8、117.9、106.8、96.1、71.1、86.2和75.3 mm。由圖1可以看出，稻季徑流水量主要集中在6月30日和7月7日兩次強降雨過程，分別達到8.92和6.23 m3(以小區(qū)計，60 m2)，即徑流量達1487.55和1039.35 m3·hm-2，兩者之和占稻季總徑流量（5705.55 m3·hm-2）的44.3%。同時，通過回歸相關(guān)分析可知，降雨量與徑流量呈顯著的線性相關(guān)關(guān)系，其回歸方程為y=0.8782x-17.327(R2=0.9415)。

圖1 徑流水量變化規(guī)律Fig.1 The quantities of runoff water

2.2 徑流水養(yǎng)分質(zhì)量濃度變化規(guī)律

由圖2（A）可以看出，稻季徑流水全氮（TN）質(zhì)量濃度隨著水稻生育期而呈逐漸下降的趨勢，6月22日雖然徑流量很小，但是徑流水中TN質(zhì)量濃度極高，極顯著高于其它幾次徑流過程（P＜0.01），各處理條件下，TN質(zhì)量濃度由高到低依次為常規(guī)施肥(8.42 mg·L-1)、肥料運籌(6.25 mg·L-1)、秸稈還田(6.19 mg·L-1)和還田減肥(5.04 mg·L-1)。6月30日和7月7日兩次強降雨過程產(chǎn)生徑流水中TN質(zhì)量濃度次之，8月3日產(chǎn)生徑流TN質(zhì)量濃度最低。同時，常規(guī)施肥條件下，幾次徑流過程平均TN濃度較秸稈還田、肥料運籌和還田減肥處理高19.87%、26.40%和28.04%(P＜0.01)。

由圖2（B）可以看出，稻季徑流水全磷（TP）質(zhì)量濃度隨著水稻生育期而呈現(xiàn)先升后降的趨勢，7月7日強降雨過程產(chǎn)生徑流中TP質(zhì)量濃度最高，顯著高于其他幾次徑流(P＜0.05)。常規(guī)施肥條件下，TP質(zhì)量濃度高達0.054 mg·L-1。同時，常規(guī)施肥條件下，幾次徑流過程平均TP質(zhì)量濃度較肥料運籌、秸稈還田和還田減肥處理高13.37%、28.12%和34.84%(P＜0.01)。

由圖 2(C)可以看出，稻季徑流水速效鉀(AK)質(zhì)量濃度在6月22日這次降雨過程中最高，且極顯著高于其他幾次徑流AK質(zhì)量濃度(P＜0.01)。各處理條件下，AK質(zhì)量濃度由高到低依次為常規(guī)施肥(56.92 mg·L-1)、秸稈還田(37.36 mg·L-1)、肥料運籌(35.58 mg·L-1)和還田減肥(31.54 mg·L-1)。同時，常規(guī)施肥條件下，幾次徑流過程平均K質(zhì)量濃度較秸稈還田和還田減肥處理高 32.03%和 45.71%(P＜0.01)。

2.3 養(yǎng)分徑流量變化規(guī)律

由圖3（A）可以看出，稻季全氮(TN)徑流主要集中在生育前期，6月30日TN徑流流失量極顯著高于其它幾次徑流過程(P＜0.01)。同時，常規(guī)施肥條件下，6月30日和7月7日兩次強降雨過程產(chǎn)生徑流量占總徑流流失量的60.8%；由圖3(B)可以看出，稻季全磷(TP)徑流主要集中在生育前期，常規(guī)施肥條件下，6月30日和7月7日兩次強降雨過程產(chǎn)生徑流流失量占總流失量的59.6%。同時，由圖3(C)可以看出，稻季速效鉀(AK)徑流也集中在生育前期，6月30日AK徑流流失量極顯著高于其它幾次徑流過程(P＜0.01)。常規(guī)施肥條件下，6月30日和7月7日兩次強降雨過程產(chǎn)生徑流量之和占總流失量的69.1%。

2.4 不同處理條件下養(yǎng)分徑流總量及流失率

圖2 徑流水養(yǎng)分質(zhì)量濃度變化規(guī)律Fig.2 The concentration of total N (A), total P (B)and K (C) in runoff water

從稻季TN總流失量角度考量(圖4A)，各處理由高到低依次為常規(guī)施肥、秸稈還田、還田減肥和肥料運籌，流失量分別達到 11.29、9.77、9.31和7.92 kg·hm-2，較常規(guī)施肥處理而言，秸稈還田、還田減肥和肥料運籌處理減少 TN流失量分別達13.48%、17.55%和 29.91%(P＜0.01)；各處理條件下，稻季TP總流失量由高到低依次為常規(guī)施肥、肥料運籌、秸稈還田和還田減肥，流失量分別達到0.19、0.17、0.14 和 0.13 kg·hm-2，較常規(guī)施肥處理而言，秸稈還田、還田減肥和肥料運籌處理減少TP流失量分別達25.00%、31.25%和12.50%(P＜0.01)；AK總流失量為常規(guī)施肥處理最高，達 13.22 kg·hm-2，秸稈還田次之，流失量分別達 10.22 kg·hm-2；還田減肥最低，為 6.15 kg·hm-2，較常規(guī)施肥處理而言，秸稈還田和還田減肥處理減少AK流失量分別達22.69%和53.48%(P＜0.01)。

圖3 養(yǎng)分徑流量變化規(guī)律Fig.3 The quantities of nutrient loss

圖4 不同處理條件下養(yǎng)分徑流總量（A）及流失率（B）Fig.4 The quantities (A) and rate (B) of nutrient loss during the whole paddy season

由圖4（B）可以看出，稻麥兩熟制農(nóng)田稻季養(yǎng)分徑流AK流失率最高，TN次之，TP最低。常規(guī)施肥條件下，AK流失率最高，達9.8%，秸稈還田和還田減肥處理次之，分別達到7.6%和5.7%。較常規(guī)施肥處理而言，秸稈還田和還田減肥處理降低AK流失率分別達22.45%和41.84%(P＜0.01)；TN流失率由高到低依次為還田減肥、常規(guī)施肥、秸稈還田和肥料運籌，分別達到 3.9%、3.8%、3.3%和2.6%，較常規(guī)施肥處理而言，秸稈還田、還田減肥和肥料運籌處理降低 TN流失率分別達 13.16%、-2.63%、31.58%；TP流失率由高到低依次為常規(guī)施肥、肥料運籌、還田減肥和秸稈還田，分別達到0.21%、0.19%、0.18%和 0.16%，較常規(guī)施肥處理而言，秸稈還田、還田減肥和肥料運籌處理降低TP流失率分別達23.81%、14.29%、9.52%(P＜0.01)。

3 討論

稻麥兩熟是江蘇省面積最大的典型種植模式，雖然作物年產(chǎn)量明顯增加，但由于長期連續(xù)種植，作物消耗的土壤養(yǎng)分也日趨增多，尤其是近幾年在增施氮素化肥的同時大幅度減施有機肥，更加劇了土壤養(yǎng)分的失調(diào)，導致土壤肥力下降。如何充分利用兩熟地區(qū)豐富的稻麥秸稈資源，以克服土壤肥力下降與作物穩(wěn)產(chǎn)高產(chǎn)的矛盾，并減輕環(huán)境污染，已成為政府和廣大農(nóng)業(yè)科技工作者極為關(guān)注的課題。從本試驗研究結(jié)果可以看出，秸稈還田和還田減肥處理較常規(guī)施肥處理均能夠有效減少稻季養(yǎng)分徑流流失總量和流失率（圖4）。究其原因，可能與還田后麥秸吸附氮、磷及麥秸腐解時微生物吸收環(huán)境中的氮、磷，從而顯著降低施肥期間稻田水氮、磷濃度有關(guān)[15]，然而，麥秸還田降低稻田水氮、磷濃度的具體機制有待于進一步研究?？梢?，秸稈還田作為一項重要的農(nóng)業(yè)管理措施值得大力推廣，從而能夠減少稻田對水網(wǎng)的潛在面源污染。

兩熟制農(nóng)田的基本特點是高投入高產(chǎn)出，由此帶來的農(nóng)田土壤養(yǎng)分流失的加劇不僅對區(qū)域土壤質(zhì)量產(chǎn)生巨大的威脅，成為進一步發(fā)展江蘇農(nóng)業(yè)生產(chǎn)的主要障礙，而且對流域水體富營養(yǎng)化的氮磷貢獻率較大[16-17]。有關(guān)農(nóng)田養(yǎng)分流失的研究報道差異很大，稻季氮素流失量相差可達百倍(0.5～54.3 kg·hm-2)，磷素相差達 5 倍(0.43～2.39 kg·hm-2)（以P2O5計）[16-19]。實際上，農(nóng)田養(yǎng)分流失特征和流失量受降雨過程、種植制度、農(nóng)作物類型、植被覆蓋度、起壟方式、地下水位、地形條件、排灌水方式、土壤類型、土壤氮磷有效養(yǎng)分含量、施肥方式及肥料運籌方式等因素的影響[20-21]。同時，從本試驗結(jié)果可以看出，肥料施用時間會直接影響徑流水中的養(yǎng)分濃度變化。距離施肥的時間越長，徑流水中氮素和速效鉀的質(zhì)量濃度就越低。如果在施肥后的短時期內(nèi)發(fā)生降水事件時，氮素的徑流損失就比較大，即氮素的徑流損失和施肥時間與降水時間及強度直接相關(guān)。若徑流發(fā)生在剛剛施肥后，徑流水中大量的氮素對附近水域的水體將構(gòu)成嚴重的威脅，而且造成了氮肥的大量損失。因此應(yīng)大力加強氣象預測和施肥管理，盡量避免在大的降水發(fā)生前夕施肥，以此來減少稻田氮素徑流損失，提高氮肥的利用率，減輕由此引起的環(huán)境負擔。

徑流水量和徑流水中養(yǎng)分質(zhì)量濃度是養(yǎng)分徑流流失的兩個重要構(gòu)成因素。本試驗中，麥秸還田、還田減肥和肥料運籌處理條件下，其徑流水量與常規(guī)施肥處理保持一致。通過回歸相關(guān)分析，徑流水量與徑流養(yǎng)分質(zhì)量濃度之間呈現(xiàn)顯著的非線性相關(guān)性。常規(guī)施肥、秸稈還田、還田減肥和肥料運籌條件下，徑流水量和徑流水中TN質(zhì)量濃度間的回歸方程分別為y=0.0018x2-0.2252x+8.4377(R2=0.9915)、y=0.0014x2-0.1653x+6.1681(R2=0.987)、y=0.0011x2-0.133x+5.03(R2=0.9481)和y=0.0011x2-0.1558x+6.2446 (R2=0.9573)；徑流水量和徑流水中TP質(zhì)量濃度間的回歸方程分別為y=-8E-07x3+0.0001x2-0.0048x+0.0413(R2=0.8233)、y= -5E-07x3+ 9E-05x2- 0.0034x+0.029(R2= 0.8274)、y=-6E-07x3+9E-05x2-0.0032x+0.0208(R2=0.6738)和y=-7E-07x3+0.0001x2-0.0044x+0.0362(R2=0.7682)；徑流水量和徑流水中AK質(zhì)量濃度間的回歸方程分別為y=0.0116x2-1.6509x+56.585(R2=0.9826)、y=0.0076x2-1.0772x+37.144(R2=0.9813)、y= 0.0062x2-0.9055x+31.367(R2=0.9828)和y=0.007x2-1.0165x+35.378(R2=0.9823)。同時，各處理條件下，徑流水中養(yǎng)分質(zhì)量濃度和養(yǎng)分流失量的變化規(guī)律不太一致（圖2，圖3），可能是徑流水量和徑流養(yǎng)分質(zhì)量濃度二者相互作用的結(jié)果。當徑流水中的養(yǎng)分質(zhì)量濃度較低而徑流水量較大時，養(yǎng)分徑流流失量也較大。通過回歸相關(guān)分析，徑流水量和徑流養(yǎng)分質(zhì)量濃度間呈現(xiàn)顯著的非線性相關(guān)。常規(guī)施肥、秸稈還田、還田減肥和肥料運籌條件下，徑流水量和徑流TN流失量間的回歸方程分別為y=0.0439x2-1.1412x+6.8577(R2=0.971)、y=0.0402x2-1.105x+4.1934(R2= 0.9892)、y=0.038x2-1.0205x+3.546 (R2=0.9891)和y=0.0178x2+0.1181x+3.1954 (R2=0.8942)；徑流水量與徑流TP流失量間的回歸方程分別為y=-4E-05x3+0.0071x2-0.2287x+0.0751 (R2=0.9381)、y=-3E-05x3+0.0047x2-0.1562x+0.0541 (R2=0.9441)、y=-3E-05x3+0.0051x2-0.1637x+0.0468 (R2=0.887)和y=-4E-05x3+0.0064x2-0.2071x+0.0664 (R2=0.9212)；徑流水量和徑流AK流失量間的回歸方程分別為y=0.0957x2-4.558x+19.43(R2=0.9519)、y=0.0646x2-2.8575x+13.349(R2=0.9001)、y=0.037x2-1.6442x+10.895(R2=0.8277)和y=0.0387x2-1.5573x+11.634(R2=0.8421)。

同時，通過回歸相關(guān)分析發(fā)現(xiàn)，徑流水中養(yǎng)分質(zhì)量濃度與徑流養(yǎng)分流失量間也存在顯著的非線性相關(guān)性。常規(guī)施肥、秸稈還田、還田減肥和肥料運籌條件下，徑流水中TN質(zhì)量濃度與徑流TN流失量間的回歸方程分別為y=-25.646x2+241.65x-211.88(R2=0.9323)、y=-41.177x2+290.98x-222.61(R2=0.9594)、y=-57.803x2+341.81x-249.82(R2=0.9287)和y=-22.58x2+152.41x-68.474(R2=0.7854)；徑流水中TP質(zhì)量濃度與徑流TP流失量間的回歸方程分別為y=-3E+11x5+4E+10x4-3E+09x3+ 9E+07x2-1E+06x+7797(R2=0.9979)，y=-7E+09x5+8E +08x4-3E+07x3+606221x2-5002.9x+15.075(R2=0.874)、y=-1E+10x5+1E+09x4-3E+07x3+485032x2-2997.8x+6.7602 (R2=0.9667)和y=-4E+09x5+5E+08x4-2E+07x3+496177x2-4359.2x+14.077(R2=0.6152)。同時，養(yǎng)分徑流流失總量包括徑流水中所含養(yǎng)分和流失泥沙中所富集的養(yǎng)分。本研究結(jié)果僅包含徑流水養(yǎng)分流失部分，流失泥沙對養(yǎng)分富集特征及各處理對徑流泥沙量的影響另文報道。

4 結(jié)論

(1) 稻麥兩熟農(nóng)田稻季徑流水量達 5705.55 m3·hm-2；常規(guī)施肥條件下，稻田徑流流失養(yǎng)分總量分別為全N 11.29 kg hm-2、全P 0.19 kg hm-2、速效K 13.22 kg hm-2；養(yǎng)分流失率分別為3.8%、0.21%和9.8%。

(2) 稻麥兩熟農(nóng)田稻季徑流水中全N和速效K濃度隨距施肥時間的延長而呈現(xiàn)逐漸下降趨勢，全P濃度先升后降。

(3) 較常規(guī)施肥處理而言，秸稈還田和還田減肥處理不僅能夠有效降低徑流水中全N、全P和速效K質(zhì)量濃度，而且能夠減少稻季養(yǎng)分徑流流失總量，分別減少 13.48%、17.55%，25.00%、31.25%和 22.69%、53.48%，并降低養(yǎng)分流失率，分別達13.16%、-2.63%，23.81%、14.29%和 22.45%、41.84%。

[1]SHARPLEY A N, DANIEL T, SIMS T, et al.Agricultural phosphorus and eutrophication[M].Second edition.US Department of Agriculture,Agricultural Research Service, ARS.2003: 149

[2]BOESCH D F, BRINSFIELD R B, MAGNIEN R E.Chesapeake bay eutrophication: scientific understanding, ecosystem restoration, and challenges for agriculture[J].Journal of Environmental Quality, 2000,30(2), 303-320.

[3]TABUCHI T.Consideration of water quality conservation of lakes[R].Gihodoshuppan Co.Ltd., Tokyo, Japan, 2005.

[4]張志劍, 王珂, 朱蔭湄, 等.水稻田表水磷素的動態(tài)特征及其潛在環(huán)境效應(yīng)的研究[J].中國水稻科學, 2000, 14(1): 55-57.ZHANG Zhijian, WANG Ke, ZHU Yinmei, et al.Dynamic characteristics of phosphorus in surface water of paddy field and its potential environmental impact[J].Chinese Rice Science, 2000, 14(1): 55-57.

[5]張志劍, 董亮, 朱蔭湄.水稻田面水氮素的動態(tài)特征、模式表征及排水流失研究[J].環(huán)境科學學報, 2001, 21(4): 475-480.ZHANG Zhijian, DONG Liang, ZHU Yinmei.The dynamic characteristics and modeling of nitrogen in paddy field surface water and nitrogen loss from field drainage[J].Acta Scientiae Circumstantiae,2001, 21(4): 475-480.

[6]KUMAR R, ARNBASHT R S, SRIVASTAVA A, et al.Reduction of nitrogen losses through erosion byLeonotis nepetaefoliaandSida acutain simulated rain intensities[J].Ecological Engineering, 1997, 8:233-239.

[7]ZOBISCH M A, RICHTER C, HEILIGTAG B, et al.Nutrient losses from cropland in the Central Highlands of Kenya due to surface runoff and soil erosion[J].Soil and Tillage Research, 1999, 33: 109-116.

[8]李恒鵬, 劉曉玫, 黃文鈺.太湖流域浙西區(qū)不同土地類型的面源污染產(chǎn)出[J].地理學報, 2004, 59(3): 403-408.LI Hengpeng, LIU Xiaomei, HUANG Wenyu.The non-point output of different landuse types in Zhexi hydraulic region of Taihu basin[J].Acta Geographica Sinica, 2004, 59(3): 403-408.

[9]單保慶, 尹澄清, 白穎, 等.小流域磷污染物非點源輸出的人工降雨模擬研究[J].環(huán)境科學學報, 2000, 20(1): 33-37.SHAN Baoqing, YIN Chengqing, BAI Ying, et al.Study on phosphorus load from a watershed with rainfall simulation method[J].Acta Scientiae Circumstantiae, 2000, 20(1): 33-37.

[10]黃滿湘, 章申, 張國梁, 等.北京地區(qū)農(nóng)田氮素養(yǎng)分隨地表徑流流失機理[J].地理學報, 2003, 58(1): 147-154.HUANG Manxiang, ZHANG Shen, ZHANG Guoliang, et al.Losses of nitrogen nutrient in overland flow from farmland in Beijing under simulated rainfall conditions[J].Acta Geographica Sinica, 2003, 58(1):147-154.

[11]康玲玲, 朱小勇, 王云璋, 等.不同雨強條件下黃土性土壤養(yǎng)分流失規(guī)律研究[J].土壤學報, 1999, 36(4): 536-543.KANG Lingling, ZHU Xiaoyong, WANG Yunzhang, et al.Research on nutrient loss from a loessial soil under different rainfall intensities[J].Acta Pedologica Sinica, 1999, 36(4): 536-543.

[12]陳皓, 章申.黃土地區(qū)氮磷流失的模擬研究[J].地理科學, 1991,11(2): 142-148.CHEN Hao, ZHANG Shen.Nitrogen and phosphorus movement from loess tillage under modeling the rainfall-runoff erosion processes[J].Scientia Geographica Sinica, 1991, 11(2): 142-148.

[13]魯如坤.土壤農(nóng)業(yè)化學分析方法[M].北京: 中國農(nóng)業(yè)科技出版社,2000.LU Rukun.Analytical method of soil agrochemistry[M].Beijing:China Agricultural Science and Technology Press, 2000.

[14]鮑士旦.土壤農(nóng)化分析(第三版)[M].北京: 中國農(nóng)業(yè)出版社, 2000.BAO Shidan.Soil and agricultural chemistry analysis[M].3Th edition.Beijing: China Agricultural Press, 2000.

[15]朱普平, 常志州, 孫麗, 等.麥草還田及耕作方式對稻田水氮、磷濃度和水稻產(chǎn)量的影響[J].江蘇農(nóng)業(yè)科學, 2004, 6: 151-153.ZHU Puping, CHANG Zhizhou, SUN Li, et al.Effects of wheat straw returning and soil tillage on nitrogen, phosphorus concentrations in the surface water of paddy field and rice yield[J].Jiangsu Agricultural Sciences, 2004, 6: 151-153.

[16]秦祖平, 徐琪.水循環(huán)對太湖地區(qū)稻田生態(tài)系統(tǒng)中養(yǎng)分平衡的影響[J].農(nóng)村生態(tài)環(huán)境, 1989, (2): 22-25.QIN Zuping, XU Qi.Effects of hydrological cycle on the balance of nutrients in rice ecosystems in the Taihu lake region[J].Rural Eco-Environment, 1989, (2): 22-25.

[17]馬立珊.蘇南太湖水系農(nóng)業(yè)非點源氮污染及其控制對策研究[J].應(yīng)用生態(tài)學報, 1992, 3(4): 346-354.MA Lishan.Nitrogen pollution from agricultural non-point sources and its control in water system of Taihu lake[J].Chinese Journal of Applied Ecology, 1992, 3(4): 346-354.

[18]張維理, 徐愛國, 冀宏杰, 等.中國農(nóng)業(yè)面源污染形勢估計及控制對策 III.中國農(nóng)業(yè)面源污染控制中存在問題分析[J].中國農(nóng)業(yè)科學, 2004, 37(7): 1026-1033.ZHANG Weili, XU Aiguo, JI Hongjie, et al.Estimation of agricultural non-point source pollution in China and the alleviating strategies.A Ⅲreview of policies and practices for agricultural non-point source pollution control in China[J].Scientia Agricultura Sinica, 2004, 37(7):1026-1033.

[19]張煥朝, 張紅愛, 曹志洪.太湖地區(qū)水稻土磷素徑流流失及其Olsen磷的“突變點”[J].南京林業(yè)大學學報(自然科學版), 2004,28(5): 6-10.ZHANG Huanchao, ZHANG Hongai, CAO Zhihong.Research on phosphorus runoff losses from paddy soils in the Taihu lake region and its Olsen-P “Change-point”[J].Journal of Nanjing Forestry University:Natural Sciences Edition, 2004, 28(5): 6-10.

[20]VAN DER MOLEN D T, BREEUWSMA A, BOERS P C M.Agricultural nutrient losses to surface water in the Netherlands: impact, strategies, and perspectives[J].Journal of Environmental Quality, 1998, 27:4-11.

[21]RITTER W F, SHIRMOHAMMADI A.Agricultural non-point source pollution: Watershed management and hydrology[M].CRC Press,Printed in the United states of America, 2001.