張 燕，高 翔，張 洪

（南京大學(xué) 地理與海洋科學(xué)學(xué)院，南京210023）

湖泊水體富營養(yǎng)化常與流域農(nóng)業(yè)資源的低效利用及其排放有關(guān)，不少學(xué)者對這方面進(jìn)行了研究［1－7］，但流域農(nóng)業(yè)資源投入及利用效率與湖泊水質(zhì)之間的關(guān)聯(lián)性存在與否及其強(qiáng)弱的定量研究卻不多見。

太湖是我國第三大淡水湖。太湖流域位于東經(jīng)119°11′—121°53′，北緯30°28′—32°15′，地處長江三角洲南翼，北濱長江，南瀕錢塘江，東臨東海，西以天目山、茅山等山區(qū)為界；湖區(qū)面積約2 000km2，流域面積36 895km2；行政區(qū)劃分屬江蘇、浙江、上海和安徽［8］。2009年太湖水體總體營養(yǎng)狀態(tài)為中度富營養(yǎng)［9］。有研究認(rèn)為，非點源污染是太湖水體富營養(yǎng)化的主要原因［7，10］。

作為我國古老的農(nóng)業(yè)區(qū)之一，太湖流域農(nóng)業(yè)發(fā)達(dá)，墾殖指數(shù)高，農(nóng)業(yè)利用集約，高投入、高產(chǎn)出是當(dāng)今太湖地區(qū)種植業(yè)的主要特點。尤其無錫市人多地少，人均耕地面積僅0．02hm2／人；土地利用強(qiáng)度高，2010年復(fù)種指數(shù)達(dá)172．6%；糧食、油料、蔬菜產(chǎn)量為80．44×107，1．16×107，126．07×107kg；當(dāng)年無錫市化肥平均施用量［11］（表1）高于同期全國化肥平均用量（折純量）346．1kg／hm2［12］，更遠(yuǎn)高于國際公認(rèn)的化肥安全施用上限225kg／hm2［13］。因此，農(nóng)業(yè)活動，尤其農(nóng)業(yè)資源的低效利用很可能是水體的主要非點源污染源。

表1 2010年無錫單位播種面積化肥施用量 kg／hm2

為更有針對性地對太湖水體營養(yǎng)化進(jìn)行控制和治理，本文根據(jù)無錫市不同地區(qū)的農(nóng)業(yè)資料，估算各地區(qū)多年的化肥利用效率；用多元相關(guān)分析探討無錫市各地區(qū)農(nóng)業(yè)資源投入與太湖水體中總磷、總氮濃度（CTP、CTN）的關(guān)聯(lián)性，以便掌握不同地區(qū)農(nóng)業(yè)活動與水體富營養(yǎng)化聯(lián)系的緊密程度，并由此判別太湖水體CTP、CTN的主要農(nóng)業(yè)來源區(qū)域。

1 材料與方法

1．1 指標(biāo)、因子選擇

1．1．1 水質(zhì)指標(biāo) 本文以太湖的營養(yǎng)物質(zhì)總磷與總氮濃度CTP、CTN為水質(zhì)指標(biāo)。

1．1．2 農(nóng)業(yè)投入因子 岸基物質(zhì)影響湖泊水體水質(zhì)需同時具備兩類因子：源物質(zhì)和傳輸它們的動力。具體到農(nóng)業(yè)生產(chǎn)，導(dǎo)致農(nóng)田面源污染通常是因為投入的資源有富余，多余的資源由徑流或流失土壤攜帶匯入了水體，可見，有剩余才有可供流失的物質(zhì)，即相對于產(chǎn)出而言，未被完全利用的資源是引起水環(huán)境問題的根源?？紤]到農(nóng)業(yè)投入對水體水質(zhì)的影響主要應(yīng)來自于化肥的殘留（剩余），途徑則主要為農(nóng)田降雨徑流與灌溉排水引起的流失。其中，化肥剩余是影響水質(zhì)的物質(zhì)來源；而排水則起著輸送源物質(zhì)的作用，但因缺少農(nóng)業(yè)生產(chǎn)中實際排水量的資料，本文采用有效灌溉面積間接代替農(nóng)田排水量。因此，本文選擇的影響水質(zhì)的農(nóng)業(yè)資源因子為化肥剩余、有效灌溉面積。

化肥剩余的計算?；适寝r(nóng)作物生長所需營養(yǎng)元素的主要來源之一，它提供的營養(yǎng)元素會轉(zhuǎn)化成農(nóng)作物最終產(chǎn)物的組成部分，因此，化肥施用量的多少并不是造成環(huán)境問題的根源，只有它的利用效率才會引發(fā)環(huán)境問題：過量施用化肥可能導(dǎo)致營養(yǎng)元素未被作物吸收而在土壤表層富積，化肥中這部分未被作物吸收從而沒能轉(zhuǎn)化成農(nóng)作物最終產(chǎn)物的營養(yǎng)元素（化肥剩余）才可能成為湖泊的潛在“營養(yǎng)物質(zhì)源”。

式中：rs（t）——t年未被利用的化肥剩余量；r（t）——t年化肥投入量；η（t）——t年農(nóng)作物對化肥的利用效率；η（t）×r（t）——t年化肥投入量中已經(jīng)轉(zhuǎn)化成農(nóng)作物最終產(chǎn)物的部分。

式（1）中的η（t）很難直接測量，因而缺少相關(guān)數(shù)據(jù)，但作物對資源的利用效率與生產(chǎn)力成正比，故可用產(chǎn)出與資源投入的對比（資源的生產(chǎn)力）來間接度量作物對資源的利用效率：單位資源形成的綜合產(chǎn)出Q／r越多，則作物對資源的利用效率越高；t時期的生產(chǎn)力Q（t）／r（t）對應(yīng)此時期的效率η（t），基準(zhǔn)期t0的生產(chǎn)力Q（t0）／r（t0）則對應(yīng)當(dāng)時的效率η（t0），即有：

式中的綜合產(chǎn)出Q用下文1．3所述TOPSIS方法計算。

于是，用資源生產(chǎn)力的相對變化代替效率的相對變化，t年未被利用的剩余量rs（t）為：

有效灌溉面積的確定。要有動力，農(nóng)田投入資源中未被利用的剩余物質(zhì)才會流失并進(jìn)入湖泊，太湖地區(qū)輸送農(nóng)田剩余物質(zhì)的主要動力便是排水量。

總用水量W取決于單位面積用水量w和有效灌溉面積S，即W＝w×S。如果w不變，則W∝S。所以，農(nóng)田排水量主要取決于有效灌溉面積S。在缺少排水量直接數(shù)據(jù)的條件下，有效灌溉面積S來間接度量排水量。

1．2 研究方法

源物質(zhì)和傳輸它們的動力共同作用才會對湖泊水質(zhì)產(chǎn)生影響，因而，本文采用復(fù)相關(guān)分析方法，借助相關(guān)系數(shù)數(shù)值的比較，大致推斷各地的相關(guān)因子對湖泊水體CTP、CTN的影響及其大小。

因為留在農(nóng)田中的化肥剩余要經(jīng)過一定時間的遷移才會進(jìn)入水體，且受各種條件限制，各地化肥剩余遷移的時間有長有短，故本文用“滯后相關(guān)系數(shù)”來反映各地化肥剩余在遷移時間上的差別，即選擇不同時間序列化肥剩余的數(shù)據(jù)計算相關(guān)系數(shù)。

此外，本文中還采用“排除法”：（1）化肥剩余是水體營養(yǎng)物的來源之一，剩余越多則污染水體的風(fēng)險越大；有效灌溉面積的擴(kuò)大增加了農(nóng)田排水量，可能增加營養(yǎng)元素的流失，從而提高水體富營養(yǎng)化的風(fēng)險。因此，根據(jù)常識，各項農(nóng)業(yè)投入因子要對水體中的CTP、CTN產(chǎn)生影響，它們之間的相關(guān)性就應(yīng)表現(xiàn)為：化肥剩余、有效灌溉面積應(yīng)該與CTP、CTN呈正相關(guān)性。如果不是如此，只能認(rèn)為相應(yīng)農(nóng)業(yè)投入因子對CTP、CTN沒有影響。（2）如果化肥剩余與CTP或CTN呈負(fù)相關(guān)性，即使有效灌溉面積與CTP、CTN具有正相關(guān)性，由于沒有了“物源”，亦應(yīng)認(rèn)為該地對太湖水體CTP、CTN無影響。（3）農(nóng)田排水則是攜帶“源物質(zhì)”的“動力”，如果有效灌溉面積與CTP或CTN呈負(fù)相關(guān)性，表明當(dāng)?shù)厥チ伺潘當(dāng)y帶“源物質(zhì)”的能力，也認(rèn)為對應(yīng)地區(qū)對太湖水體中CTP、CTN的影響可以不計。

1．3 用TOPSIS法計算農(nóng)業(yè)綜合產(chǎn)出Q

用 TOPSIS法（Technique for Order Preference by Similarity to Ideal Solution）［14］計算農(nóng)業(yè)綜合產(chǎn)出Q的思路是：以各地的農(nóng)業(yè)產(chǎn)出作為樣本，先建立由m個樣本（即m個地區(qū)）和n個指標(biāo)（即n種農(nóng)業(yè)產(chǎn)出）組成的討論對象，則n個指標(biāo)構(gòu)成一個n維空間，每個樣本對應(yīng)此空間的一個特定點，此空間中共有m個對應(yīng)點；再從此空間中找兩個參照點：所有指標(biāo)的最大值和最小值分別構(gòu)成的點X＋（最大參照樣本）、X－（最小參照樣本）；然后分別計算每個樣本與參照樣本X＋、X－間的距離（即空間兩點間的距離）d＋、d－；最后計算指數(shù)Q作為農(nóng)業(yè)綜合產(chǎn)出。具體步驟是：

（1）選擇指標(biāo)，建立矩陣，并做數(shù)據(jù)標(biāo)準(zhǔn)化處理。根據(jù)研究目的和指標(biāo)的可獲取性，建立一個由m個樣本n個指標(biāo)構(gòu)成的矩陣X＝（xij）m×n，xij為指標(biāo)實際值，i為樣本序數(shù)，j為指標(biāo)序數(shù)。

為排除xij的量綱及數(shù)量級差異的影響，對矩陣進(jìn)行極差標(biāo)準(zhǔn)化處理，處理后的矩陣為X′＝（x′ij）m×n。

（2）用矩陣X′中各元素的最大值和最小值確定最大、最小參照樣本X＋、X－。對產(chǎn)出而言，最大參照樣本對應(yīng)最理想的樣本狀態(tài)，最小參照樣本對應(yīng)最不理想的樣本狀態(tài)。

（3）計算各樣本與X＋、X－間的加權(quán)歐氏距離d＋、d－。

因為不同作物的產(chǎn)出量不同，從而對綜合產(chǎn)出的貢獻(xiàn)不同，而產(chǎn)出量與播種面積成比例，于是，權(quán)系數(shù)wj就根據(jù)第j種作物播種面積Aj占總播種面積A的比例來確定，即wj＝Aj／A［15］。

一個樣本所有指標(biāo)均達(dá)最大值，它與最大參照樣本的距離便最小（d＋＝0），與最小參照樣本的距離最大（d－＝1）；而一個樣本所有指標(biāo)均為最小值時，它與最大參照樣本的距離便最大（d＋＝1），與最小參照樣本的距離最?。╠－＝0）。對產(chǎn)出而言，d＋i越小表示i樣本與期望越接近，d－i越大表示i樣本越遠(yuǎn)離最不理想狀態(tài)，則i樣本狀態(tài)越好。

（4）用指數(shù)Qi綜合d＋i與d－i。對產(chǎn)出來說，Qi越大，表示i樣本狀態(tài)越優(yōu)。若樣本各指標(biāo)均處于最優(yōu)（即d＋＝0）則Qi＝1；若樣本各指標(biāo)均處于最劣（即d－＝0）則Qi＝0。

2 結(jié)果與分析

2．1 無錫市化肥利用效率和化肥剩余量的時空演變

利用無錫各區(qū)市農(nóng)業(yè)資料［11］，由式（5）—（8）綜合糧食、油料與蔬菜得到農(nóng)業(yè)綜合產(chǎn)出Q。再根據(jù)我國化肥利用效率的已有研究［16－17］及與發(fā)達(dá)國家化肥施用量［18］對比得到的間接結(jié)果，由式（2）、（4）計算無錫各區(qū)市化肥利用效率和化肥剩余量，其中，化肥折純、氮肥與磷肥的η（t0）分別取0．35，0．40，0．15，Q（t0）／r（t0）取2000—2010年Q（t）／r（t）的平均值，結(jié)果見表2、表3。

由表2可見，2000—2010年，市轄區(qū)與江陰化肥利用效率變化的總趨勢是波動中快速下降，只是近兩年才有稍許提高；宜興化肥利用效率在2003年前呈逐年下降，之后在波動中逐步提高，近年的提高更加顯著。就各區(qū)市的對比來看，各地的化肥利用效率大致相同，宜興稍高一些。2001—2002年無錫各區(qū)市的化肥利用效率都相對較高，是因為這一期間各地的農(nóng)業(yè)綜合產(chǎn)出均較高；2003年各區(qū)市的化肥利用效率都較低（宜興尤其明顯），可能是因當(dāng)年降水量少［11］，大面積受災(zāi)造成農(nóng)業(yè)產(chǎn)出低所致；2004—2005年宜興化肥利用效率低是因為當(dāng)年化肥施用量提高而產(chǎn)出卻未同步增長；2008年市轄區(qū)與江陰的農(nóng)業(yè)綜合產(chǎn)出下降，導(dǎo)致這一期間的化肥利用效率也相對低下。農(nóng)業(yè)在農(nóng)村經(jīng)濟(jì)中占比高的地區(qū)，化肥利用效率也相對高，化肥利用效率從高至低為宜興、江陰、市轄區(qū)，這可能與農(nóng)技普及水平有關(guān)。

表2 2000－2010年無錫各區(qū)市化肥利用效率 %

表3 2000－2010年無錫各區(qū)市化肥剩余量107 kg

由表3可見，化肥剩余量存在區(qū)域差異，總體來看，2003年以前化肥剩余量逐年增加，且增速較快，剩余量高的地區(qū)是宜興與市轄區(qū)；而2003年以后化肥剩余量逐年降低，尤其是近幾年宜興的降速較快，化肥剩余量的區(qū)域差異也相對縮??；按化肥剩余量從多到少的排序為宜興、市轄區(qū)、江陰。化肥剩余量的變化趨勢也存在空間差異：市轄區(qū)為先逐年快速增加，而后趨于穩(wěn)定并逐步下降；宜興先逐年緩慢增加，而后快速下降；江陰相對穩(wěn)定，年際間變化不大，且化肥剩余量最低。

引起各地化肥剩余量空間與時間變化的主要因素是農(nóng)作物播種面積，播種面積大的地區(qū)化肥剩余總量相對高，歷年播種面積的總體降序排列為宜興、市轄區(qū)、江陰，排位靠前的宜興與市轄區(qū)的化肥剩余量就相對高些。其次是化肥施用強(qiáng)度，經(jīng)濟(jì)相對發(fā)達(dá)、且產(chǎn)業(yè)結(jié)構(gòu)中農(nóng)業(yè)所占比重高的地區(qū)對流域化肥剩余量的貢獻(xiàn)更大一些，這與各地經(jīng)濟(jì)不斷發(fā)展、化肥施用逐步增加的趨勢有關(guān)，經(jīng)濟(jì)相對發(fā)達(dá)的市轄區(qū)的單位播種面積化肥施用量為547．3kg／hm2；第一產(chǎn)業(yè)占比相對高的地區(qū)的單位播種面積化肥施用量也相對高，如宜興化肥施用量為434．3kg／hm2［11］。另外還與化肥施用效率的變化有關(guān)，如2003年各區(qū)市的化肥利用率低，當(dāng)年的化肥剩余量高，再如江陰2008年的化肥利用率相對低，則化肥剩余量相對高，2007—2010年宜興的化肥利用率大幅提高，對應(yīng)年份的化肥剩余量急劇下降。化肥剩余量還與當(dāng)年化肥施用的絕對量有關(guān)，盡管2001—2002年各區(qū)市的化肥利用效率均較高，但因這一期間各區(qū)市的化肥施用量高，化肥剩余量也相對高。

2．2 太湖水質(zhì)的農(nóng)業(yè)影響源識別

表4為2000—2009年太湖3個湖區(qū)（梅梁湖、五里湖與西部沿岸）的CTP、CTN［19］。據(jù)此計算不同湖區(qū)CTP、CTN與各地農(nóng)業(yè)投入因子的偏相關(guān)系數(shù)以及復(fù)相關(guān)系數(shù)，結(jié)果見表5。

表4 2000－2009年太湖水質(zhì)指標(biāo) mg／L

表5 太湖水體CTP、CTN與無錫農(nóng)業(yè)投入因子的偏相關(guān)系數(shù)與復(fù)相關(guān)系數(shù)

（1）在一些地區(qū)，水體CTP或CTN只與當(dāng)年磷肥剩余或氮肥剩余呈正相關(guān)，如西部沿岸CTP與宜興磷肥剩余的關(guān)系，可能是由于宜興地處西部沿岸上游，且為丘陵地區(qū)，這樣的地形有利于磷肥的流失。在另一些地區(qū)，水體CTP或CTN與當(dāng)年磷肥或氮肥剩余不存在相關(guān)關(guān)系，但與滯后1a的磷肥剩余或氮肥剩余的相關(guān)系數(shù)顯著相關(guān)，如梅梁湖、五里湖的CTP與江陰滯后1a的磷肥剩余有正的相關(guān)性，這可能與江陰地處市轄區(qū)以北，而附著于土壤顆粒上的磷的輸送需要一定時間有關(guān)；梅梁湖、五里湖的CTN與市轄區(qū)的氮肥剩余的相關(guān)性也是如此。此外，有些地區(qū)，水體CTP或CTN同時與當(dāng)年及滯后1a的磷肥剩余或氮肥剩余呈正相關(guān)性，如梅梁湖的CTP與宜興的磷肥剩余，五里湖的CTP與市轄區(qū)的磷肥剩余，及西部沿岸的CTN與宜興的氮肥剩余，表明這些地區(qū)農(nóng)業(yè)生產(chǎn)中施用磷肥或氮肥產(chǎn)生的剩余對水體CTP或CTN有1a的滯后影響期，原因可能是化肥剩余中的營養(yǎng)元素從陸地農(nóng)田遷移至水體要經(jīng)歷一定的時間，現(xiàn)在施用的化肥或者要過1a才對湖水中的CTP或CTN產(chǎn)生影響，或者這種影響能夠持續(xù)很長時間，1a以后還在繼續(xù)產(chǎn)生影響。

應(yīng)該說明的是，此處1a滯后期與采用數(shù)據(jù)的時間尺度有關(guān)，即探討農(nóng)業(yè)投入與太湖水體CTP或CTN間的關(guān)系時，受資料間隔時間限制，時間尺度為“年”，如果有更精細(xì)、時間更短的數(shù)據(jù)，化肥剩余與水體CTP或CTN間的滯后關(guān)系不一定是1a，也許只有若干月。

（2）應(yīng)用上文1．2的“排除法”，根據(jù)表5，江陰當(dāng)年及滯后1a的氮肥剩余與水體CTN均無合理的相關(guān)關(guān)系（即它們之間為負(fù)相關(guān)），原因可能是我國施用的氮肥以碳酸氫銨和尿素為主，而這種銨態(tài)氮肥施入土壤后易揮發(fā)［20］，且氮素流失還受降水量、氣溫以及農(nóng)事活動的影響［4］，加之所處地理位置距太湖較遠(yuǎn)，這些因素使得江陰對太湖CTN無影響，因此，從“物源”角度，可以排除江陰的氮肥施用對梅梁湖、五里湖及西部沿岸水體CTN的影響；與之類似，還可以排除市轄區(qū)當(dāng)年的氮肥剩余對西部沿岸水體CTN、以及市轄區(qū)和江陰滯后1a的磷肥剩余對西部沿岸水體CTP的影響。從“動力”角度，則可排除宜興對五里湖CTP、CTN及對梅梁湖CTN的影響，也可排除市轄區(qū)和江陰對西部沿岸水體CTP、CTN的影響。

（3）作為我國東部著名農(nóng)業(yè)區(qū)，太湖流域農(nóng)業(yè)增產(chǎn)的主要措施之一是高化肥施用強(qiáng)度，但利用效率低致使化肥剩余多，加上灌溉方式不合理造成用水量大及隨之產(chǎn)生的大量排水，導(dǎo)致高強(qiáng)度的N、P元素流失，太湖地區(qū)稻麥輪作系統(tǒng)中，氮、磷的徑流流失量分別占當(dāng)年施肥量的6．36%和1．93%［2］，說明灌排水也是影響太湖水體CTP、CTN的重要因素之一。由表5可知，市轄區(qū)與江陰的有效灌溉面積對梅梁湖、五里湖CTP有顯著影響，市轄區(qū)和宜興的有效灌溉面積分別對五里湖及西部沿岸水體CTN有顯著影響。因此，無錫農(nóng)業(yè)生產(chǎn)中有必要推廣先進(jìn)灌溉技術(shù)，以求節(jié)約水資源并降低太湖水體富營養(yǎng)化的風(fēng)險。

（4）對照表3與表5可以看出，不一定化肥剩余量多的地區(qū)對湖水水質(zhì)的影響就大，如化肥剩余量多的宜興對五里湖CTP、CTN的影響就可以排除，原因是五里湖位于宜興的東北，將剩余的磷肥與氮肥輸送到五里湖的動力可能不足；反而是磷肥剩余量相對少的江陰對梅梁湖與五里湖CTP有較大影響，可能是其地形利于化肥剩余向太湖輸送，這正反映了農(nóng)業(yè)投入資源對湖泊水質(zhì)的影響是綜合作用的結(jié)果?！笆Ｓ唷倍?，只表明潛在的影響大，但只要能從“動力源”上加以適當(dāng)控制，就可以使之不轉(zhuǎn)化為現(xiàn)實的風(fēng)險，即控制農(nóng)業(yè)投入資源對水環(huán)境的影響，應(yīng)從“物源”與“動力”兩方面入手。在“物源”方面，要減少化肥施用量，采用科學(xué)的施肥技術(shù)，在恰當(dāng)時間、合理位置施用數(shù)量適當(dāng)?shù)姆柿希贿^量施肥，如測土配方施肥。在“動力”方面，則要減少水土流失，主要措施不外乎：改造農(nóng)田為節(jié)水灌溉型農(nóng)田，以減少灌排水量。在“物源”難以立刻減少（因為有農(nóng)作物增產(chǎn)的要求）的現(xiàn)狀下，從“動力”方面努力，不失為一種能較快見效的辦法。以達(dá)到既保護(hù)水環(huán)境，又節(jié)約農(nóng)業(yè)活動中的資源投入量的目的。

（5）復(fù)相關(guān)系數(shù)的含義是：太湖CTP、CTN受多種農(nóng)業(yè)投入因素影響，是多因素共同作用的結(jié)果。由表5可見，各農(nóng)業(yè)投入因子的偏相關(guān)系數(shù)并不大，但復(fù)相關(guān)系數(shù)卻較大。根據(jù)復(fù)相關(guān)系數(shù)（并結(jié)合從“物源”和“動力”角度排除其影響的情況），各地農(nóng)業(yè)活動對湖區(qū)水質(zhì)影響較大的是：市轄區(qū)對五里湖與梅梁湖CTP和對五里湖CTN，江陰對五里湖與梅梁湖CTP，宜興對梅梁湖CTP和對西部沿岸CTN?？紤]到相對高的化肥施用強(qiáng)度與化肥剩余量，故應(yīng)特別重視市轄區(qū)對五里湖與梅梁湖CTP、CTN和宜興對西部沿岸CTP、CTN的影響。再結(jié)合表5中偏相關(guān)系數(shù)達(dá)到顯著性水平的農(nóng)業(yè)投入因子的數(shù)目判斷，市轄區(qū)的農(nóng)業(yè)活動對五里湖水質(zhì)的影響較大，宜興的農(nóng)業(yè)活動對西部沿岸水質(zhì)的影響較大。

總體來看，無錫各區(qū)市農(nóng)業(yè)活動對太湖不同湖區(qū)的水質(zhì)有不同影響，江陰主要影響五里湖與梅梁湖的CTP，而宜興主要影響西部沿岸與梅梁湖的CTP，特別是市轄區(qū)的農(nóng)業(yè)活動對五里湖與梅梁湖的CTP、CTN都有影響，因此，這些地區(qū)的農(nóng)業(yè)活動應(yīng)是太湖水環(huán)境管理重點關(guān)注的對象。

3 結(jié) 論

（1）農(nóng)業(yè)面污染源的物質(zhì)常源自化肥利用效率低下產(chǎn)生的化肥剩余，本文提出了一種有一定物理含義的計算化肥剩余量的簡明方法。由于化肥施用量多或／和化肥利用效率相對較低，無錫的宜興、市轄區(qū)的農(nóng)業(yè)活動中產(chǎn)生的化肥剩余量較大。

（2）無錫各地農(nóng)業(yè)活動對太湖水質(zhì)的影響是農(nóng)業(yè)投入因子協(xié)同作用的結(jié)果。在“物源”與“動力”的共同作用下，市轄區(qū)的農(nóng)業(yè)活動對五里湖與梅梁湖的CTP、CTN有較大影響，宜興的農(nóng)業(yè)活動對西部沿岸與梅梁湖的CTP有較大影響，江陰的磷肥施用對五里湖與梅梁湖的CTP影響較大。而因“動力”不足，宜興的農(nóng)業(yè)活動雖然產(chǎn)生的化肥剩余量最多，但對五里湖CTP、CTN的影響并不大，同樣，市轄區(qū)的農(nóng)業(yè)活動雖然產(chǎn)生較多的化肥剩余，但對西部沿岸水體CTP、CTN的影響也不大；因“物源”不足，江陰的氮肥剩余對五里湖與梅梁湖的CTN的影響不大。

（3）為降低農(nóng)業(yè)活動對水環(huán)境的影響，控制并改善太湖水體污染程度，同時滿足資源節(jié)約型發(fā)展的要求，無錫農(nóng)業(yè)生產(chǎn)中，有必要減少化肥施用量，并要特別注意提高化肥利用效率，以減少化肥剩余；還要改進(jìn)灌溉方式，建設(shè)節(jié)水型農(nóng)田，包括增加節(jié)水灌溉面積、減少灌溉水量。這樣從“源頭”從手，才可能實現(xiàn)提高經(jīng)濟(jì)效益、節(jié)約資源、控制并改善水環(huán)境多重功效的管理目標(biāo)。

［1］ 黃玉潔，張銀龍，李海東，等．太湖人工恢復(fù)濕地區(qū)植物群落建植對沉積物中氮、磷空間分布的影響［J］．水土保持研究，2011，18（5）：161－165．

［2］ 楊林章，王德建，夏立忠．太湖地區(qū)農(nóng)業(yè)非點源污染特征及控制途徑［J］．中國水利，2004，20（1）：29－30．

［3］ 唐佐芯，王克勤，李秋芳，等．等高反坡階對坡耕地產(chǎn)流產(chǎn)沙和氮磷遷移的作用研究［J］．水土保持研究，2013，20（1）：1－8．

［4］ 王鵬，高超，姚琪，等．太湖丘陵地區(qū)農(nóng)田氮素遷移的時空分布特征［J］．環(huán)境科學(xué)，2006，27（8）：1671－1675．

［5］ Bahar M M，Ohmori H，Yamamuro M．Relationship between river water quality and land use in a small river basin running through the urbanizing area of Central Japan［J］．Limnology，2008，9（1）：19－26．

［6］ Lee S W，Hwang S J，Lee S B，et al．Landscape ecological approach to the relationships of land use patterns in watersheds to water quality characteristics［J］．Landscape and Urban Planning，2009，92（2）：80－89．

［7］ 孫偉，陳雯，陳誠．水環(huán)境協(xié)同約束分區(qū)與產(chǎn)業(yè)布局引導(dǎo)研究：以江蘇省為例［J］．地理學(xué)報，2010，65（7）：819－827．

［8］ 王蘇民，竇鴻身．中國湖泊志［M］．北京：科學(xué)出版社，1998．3－21．

［9］ 國家統(tǒng)計局與環(huán)境保護(hù)部．中國環(huán)境統(tǒng)計年鑒2011［M］．北京：中國統(tǒng)計出版社，2011．

［10］ 張維理，武淑霞，冀宏杰，等．中國農(nóng)業(yè)面源污染形勢估計及控制對策I．21世紀(jì)初期中國農(nóng)業(yè)面源污染的形勢估計［J］．中國農(nóng)業(yè)科學(xué)，2004，37（7）：1008－1017．

［11］ 無錫統(tǒng)計局．2001—2011無錫統(tǒng)計年鑒［M］．北京：中國統(tǒng)計出版社，2001－2011．

［12］ 中華人民共和國國家統(tǒng)計局．中國統(tǒng)計年鑒2011［M］．北京：中國統(tǒng)計出版社，2012．

［13］ Syers J K，Johnston A E，Curtin D．Efficiency of Soil and Fertilizer Phosphorus Use［M］．Rome：Food and Agriculture Organization of the United Nations，2008．

［14］ Hwang C L，Yoon K．Multiple Attribute Decision Making：Method and Application［M］．New York：Springverlag，1981．

［15］ 張燕，高翔，張洪．農(nóng)業(yè)綜合產(chǎn)出對投入資源的敏感性研究［J］．中國農(nóng)業(yè)大學(xué)學(xué)報，2011，16（6）：169－173．

［16］ 朱昌雄，郭萍．農(nóng)業(yè)源污染與微生物修復(fù)技術(shù)［J］．現(xiàn)代化工，2008，28（10）：1－6．

［17］ 陳同斌，曾希柏，胡清秀．中國化肥利用率的區(qū)域分異［J］．地理學(xué)報，2002，57（5）：531－538．

［18］ FAO （Food and Agriculture Organization）．Fertilizer Use by Crop ［M］．Rome：Food and Agriculture Organization of the United Nations，2006．

［19］ 中國環(huán)境年鑒編輯委員會．中國環(huán)境年鑒2001—2010［M］．北京：中國環(huán)境年鑒社，2001－2010．

［20］ 張智峰，張衛(wèi)峰．我國化肥施用現(xiàn)狀及趨勢［J］．磷肥與復(fù)肥，2008，23（6）：9－12．