李 巧,周金龍,2,賈瑞亮

(1.新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院,新疆 烏魯木齊 830052;2.中國地質(zhì)大學(xué)環(huán)境學(xué)院,湖北 武漢 430074)

在我國有 7億多人口飲用地下水,有半數(shù)以上的城市、鄉(xiāng)村的經(jīng)濟(jì)發(fā)展和居民生活用水以地下水為主要供水水源。但地下水也與其它環(huán)境資源一樣其污染問題越來越嚴(yán)重。目前全國 25%的地下水體遭到污染,平原區(qū)約有 54%的地下水不符合生活用水水質(zhì)標(biāo)準(zhǔn)。因此地下水的污染問題必須引起高度關(guān)注[1]。而造成地下水污染最普遍、也最難治理的就是農(nóng)業(yè)面源污染。它具有形成過程隨機(jī)性大、影響因子復(fù)雜;分布范圍廣、影響深遠(yuǎn);過程復(fù)雜、機(jī)理模糊;潛伏周期長、危害大等特點(diǎn),因此研究和控制面源污染難度較大。地下水農(nóng)業(yè)面源污染的研究主要用于對(duì)地下水污染的預(yù)防與控制,減少由地下水農(nóng)業(yè)面源污染對(duì)人類健康和生態(tài)和諧造成的影響。下面主要就地下水農(nóng)業(yè)面源污染近十多年的研究做一較為全面的綜述。

1 地下水農(nóng)業(yè)面源污染的危害

1.1 威脅人類健康

農(nóng)業(yè)面源污染對(duì)地下水水環(huán)境安全構(gòu)成了嚴(yán)重威脅,有毒有害物質(zhì)及一些營養(yǎng)鹽經(jīng)過淋溶作用進(jìn)入地下水體造成飲用水源的污染,尤其是一些有毒有害物質(zhì)具有致癌、致突變、致畸形的性質(zhì)對(duì)人體健康構(gòu)成嚴(yán)重威脅。國內(nèi)外大量研究已表明地下水硝酸鹽含量上升最主要的原因是農(nóng)區(qū)大量施用 N素化肥和牧區(qū)畜禽糞便的分解淋失,而飲水中的硝酸鹽會(huì)導(dǎo)致“蘭嬰”綜合癥和胃癌、結(jié)直腸癌、淋巴瘤等癌癥發(fā)病率升高[2]。

1.2 破壞農(nóng)田生態(tài)平衡

地下水農(nóng)業(yè)面源污染不但使農(nóng)產(chǎn)品品質(zhì)下降而且使害蟲抗藥性增強(qiáng),害蟲天敵減少、土壤有機(jī)質(zhì)消耗加快、團(tuán)粒結(jié)構(gòu)減少保水保肥能力下降、土壤耕作層板結(jié)耕性變劣肥力下降、再生產(chǎn)能力降低。這些都會(huì)使得農(nóng)田的生態(tài)平衡遭到破壞[3]。

1.3 污染地表水體

農(nóng)業(yè)面源污染物不僅隨地表徑流直接進(jìn)入河流、湖泊或近海污染地表水體而且滲入土壤或巖石污染地下水。而被污染的地下水最終亦流入河流、湖泊或近海污染地表水體。地下水和地表水相互作用緊密相聯(lián)。在通常情況下地下水接受降雨或融雪入滲補(bǔ)給向位于地形低處的河流、湖泊或海洋排泄。除地表徑流和直接降水外,地下水是河流與湖泊水的主要來源。特別是在枯水季和枯水年河流與湖泊的水往往全部來自地下水。當(dāng)?shù)叵滤晃廴竞蠼邮艿叵滤a(bǔ)給的河流、湖泊和近海亦會(huì)遭受同樣甚至更加嚴(yán)重的污染[4-5]。

2 地下水農(nóng)業(yè)面源污染的主要影響因素

2.1 化肥施用量大,利用率低

過去十多年中全國化肥的施用量從 1990年的 2 590萬噸增加到 2005年的 4 897萬噸,占全世界平均消費(fèi)量的 1/4,達(dá) 400 kg/hm2,遠(yuǎn)遠(yuǎn)超過國際上為防止水體污染而設(shè)置的225 kg/hm2的化肥使用安全上限。中國每年農(nóng)田養(yǎng)分被植物利用的部分很少,氮肥的利用率僅為 30%～35%,磷肥為10%～25%,中國化肥的平均施用量是發(fā)達(dá)國家化肥安全施用上限的 2倍,平均利用率僅為 40%左右[6]。剩余的養(yǎng)分中的氮、磷、鉀元素隨水流進(jìn)入溝渠,再匯集江、河、湖、水庫及近海域,使水體中的氮、磷等營養(yǎng)元素富集,導(dǎo)致水質(zhì)的惡化;進(jìn)入土壤,改變?cè)型寥赖慕Y(jié)構(gòu)和特性,造成板結(jié);還有很大一部分會(huì)通過淋溶等途徑進(jìn)入地下水環(huán)境污染地下水[7-9]。

2.2 農(nóng)藥使用不科學(xué),形成流失

中國是世界生產(chǎn)和使用農(nóng)藥的第一大國,2000年全國農(nóng)藥使用量約為 128萬 kg,平均值為 14 kg/hm2,比發(fā)達(dá)國家高出 1倍,高毒、高殘留的有機(jī)磷、有機(jī)氯農(nóng)藥仍占很大比例。全國農(nóng)藥噴施約有 60%～70%進(jìn)入環(huán)境中,僅有約 30%被農(nóng)作物吸收。據(jù)調(diào)查,若噴灑液體農(nóng)藥,僅有 20%左右附著在植物體上,1%～4%接觸到目標(biāo)害蟲,5%～30%的藥劑飄游于空中,其余 40%～60%降落到地面,這些剩余的有害物通過大氣、地下水和地表徑流遷移,污染大氣、地下水以及地表水[10]。

2.3 禽畜糞便和農(nóng)村居民生活垃圾污染

中國的畜禽養(yǎng)殖業(yè)近年來發(fā)展迅速,畜禽糞便主要污染物化學(xué)耗氧量(COD)和生物耗氧量(BOD)的流失量逐年增加,畜禽糞便如隨意堆放且不及時(shí)妥善處理,在降雨季節(jié),污染物隨地表徑流和農(nóng)田滲漏進(jìn)入地表和地下水體會(huì)形成大面積污染。農(nóng)村普遍缺乏基本排水和生活垃圾處理系統(tǒng),生活污水大多隨意潑灑,大量蔬菜、秸稈等生產(chǎn)垃圾與生活垃圾四處堆放,在雨水的沖刷下使大量的滲濾液排入地下水體,形成地下水農(nóng)業(yè)面源污染。

2.4 農(nóng)田灌溉排水和污水灌溉

長期灌溉使土壤中的鹽類淋溶和肥料的淋濾,而導(dǎo)致土壤的鹽堿化,地下水含鹽量上升[11-12]。一些地區(qū)用污水灌溉,也使大量污染物積累。雖然污水中的營養(yǎng)元素可以為農(nóng)作物吸收,但如果施用量過大或時(shí)間不恰當(dāng),許多污水未經(jīng)過農(nóng)作物和土壤的自然凈化而直接進(jìn)入水體,同樣會(huì)導(dǎo)致土壤和地表水及地下水體污染[13]。

2.5 土地利用方式和水環(huán)境惡化

國內(nèi)外許多研究表明土地利用方式和土壤耕作措施不當(dāng)會(huì)引起嚴(yán)重的水土流失并加重地下水農(nóng)業(yè)面源污染的程度。近年來我國雖然加大了生態(tài)保護(hù)力度,但由于資金有限和人們思想上還沒有引起足夠重視,水土流失治理力度不夠,山區(qū)水土流失仍較嚴(yán)重。發(fā)生水土流失后在較大范圍內(nèi)泥沙攜帶水中的化學(xué)物質(zhì)逐步匯集并通過河流等通道注入各種水體,其結(jié)果是一方面形成水體中懸浮物,增加了水的濁度構(gòu)成物理污染;另一方面來自于土壤的泥沙作為人類施用肥料(化肥和有機(jī)肥)的載體也將大量污染物帶入水體增加了地表水的富營養(yǎng)化程度及地下水的硝態(tài)氮含量構(gòu)成地表水和地下水的化學(xué)污染?？梢娪捎谕恋乩媒Y(jié)構(gòu)不合理導(dǎo)致的水土流失是引起地下水面源污染的重要方面[14]。

3 地下水農(nóng)業(yè)面源污染物的測(cè)定方法

地下水農(nóng)業(yè)面源污染的主要污染物為氮和磷。

3.1 氮的測(cè)定

可采用國標(biāo)法(GB 11894-89)來測(cè)定地下水樣中的氮含量。總氮：堿性過硫酸鉀消解——紫外分光光度法;硝態(tài)氮：紫外分光光度法;亞硝態(tài)氮：重氮偶合比色法;銨態(tài)氮：水楊酸——次氯酸鹽光度法。我國研究者對(duì)如何對(duì)地下水中硝酸鹽氮、亞硝酸鹽氮的含量進(jìn)行符合精度要求的測(cè)定也進(jìn)行了大量的研究,提出的測(cè)定地下水中硝酸鹽氮的方法主要有：離子色譜法、一階導(dǎo)數(shù)紫外分光光度法、還原 -紫外分光光度法和 2、6-二甲酚分光光度法等。提出的測(cè)定地下水中亞硝酸鹽氮的方法有：催化光度法測(cè)定痕量亞硝酸根——嗅酸鉀 -結(jié)晶紫體系、催化光度法測(cè)定痕量亞硝酸根——嗅酸鉀 -維多利亞藍(lán)體 B系、非平衡流動(dòng)注射光度法、硝酸根電極法、新紫外光度法和堿性品紅紫外光度法等[15]。

3.2 磷的測(cè)定

水樣關(guān)于磷的分析所需指標(biāo)(總磷、顆粒態(tài)磷、水溶性總磷、磷酸鹽)的測(cè)定方法均可采用國標(biāo)。地下水樣不需過濾,可直接用鉬銻抗比色法測(cè)定磷酸鹽(PO3-4)含量,總磷(TP)含量采用鉬酸銨分光光度法測(cè)定(GB 11893-89)[16]。

4 地下水農(nóng)業(yè)面源污染的數(shù)值模擬

數(shù)值模擬是研究地下水環(huán)境最有效的方法之一,現(xiàn)已開發(fā)出 Visua lMODFLOW、GMS等專用于地下水模擬的系統(tǒng),在國內(nèi)外得到廣泛認(rèn)可和應(yīng)用。但是,如何降低模型的不確定性成為數(shù)值模擬面臨的主要挑戰(zhàn)[17]。

隨著計(jì)算機(jī)技術(shù)的快速發(fā)展,利用數(shù)學(xué)方法來解決面源污染問題成為越來越值得研究的重點(diǎn)。以氮負(fù)荷為例：地表的面源氮負(fù)荷污染地下水需要經(jīng)過 2個(gè)過程——首先經(jīng)降雨淋洗進(jìn)入包氣帶,在這一過程中會(huì)發(fā)生氮的礦化、固定化、硝化、反硝化、被植物吸收等作用,引起實(shí)際進(jìn)入到地下水中的氮負(fù)荷的變化;穿過包氣帶,氮素垂向遷移進(jìn)入飽和帶,經(jīng)過對(duì)流、擴(kuò)散并發(fā)生一系列反應(yīng)后,此時(shí)的氮負(fù)荷才是污染地下水的真實(shí)值。因此,對(duì)地表氮負(fù)荷污染地下水的數(shù)值模擬可以分為兩個(gè)部分來考慮,非飽和帶中垂向遷移模型和地下水污染遷移模型。

從 20世紀(jì) 70年代人們開始對(duì)面源模型進(jìn)行研究以來,已經(jīng)開發(fā)了很多模型可以用來模擬農(nóng)業(yè)面源污染。具有代表性的模型有：由美國普度大學(xué)農(nóng)業(yè)工程系 Beasley和 Huggins于 1981年提出的基于降水事件的分布式參數(shù)模型 ANSWERS;由Hanson等于 1981年提出的一個(gè)允許用水動(dòng)力學(xué)和沉積化學(xué)共同作用來模擬路面和土壤污染物徑流過程的物理分布型綜合模型 HSPF;由美國農(nóng)業(yè)部農(nóng)業(yè)研究署(USDA-ARS)和土壤保護(hù)署于 1986年聯(lián)合開發(fā)的基于降水事件的分布式參數(shù)模型AGNPS;由丹麥水工試驗(yàn)所(Danish Hydraulic Institute)1982年開發(fā)的 MIKE SHE;由W illiams等和 Arnold等開發(fā)提供的物理分布模型 SWRRB/SWAT等等[18]。但是,這些模型大部分是模擬研究面源污染在水平方向上的橫向遷移,很少或沒有考慮農(nóng)業(yè)面源污染的垂向遷移及其對(duì)地下水的影響。在現(xiàn)有的面源模型中,涉及到污染物在包氣帶中遷移轉(zhuǎn)化的模擬建模的理論基礎(chǔ)都是由達(dá)西定律和連續(xù)方程導(dǎo)出的土壤水分子運(yùn)動(dòng)基本方程以及污染物運(yùn)動(dòng)的水動(dòng)力彌散方程[19],從而建立數(shù)學(xué)模型進(jìn)行模擬和預(yù)測(cè)。常用的模型有 GLEAMS、HYDRUS-1D等,均可以在包氣帶中得到較好的模擬效果[20-21]。飽和帶中氮的模擬是通過建立地下水動(dòng)力學(xué)方程,并在實(shí)驗(yàn)室測(cè)定模型需要的參數(shù),對(duì)地下水的水質(zhì)水量等進(jìn)行模擬。

數(shù)值模擬法以其有效性、靈活性和相對(duì)廉價(jià)性逐漸成為地下水研究領(lǐng)域的一種不可或缺的重要方法,并受到越來越多的重視和廣泛的應(yīng)用[22]。目前較常使用的是 Visua lMODFLOW、GMS、FEFLOW等模型,在一些地區(qū)已經(jīng)有了較成功的模擬范例[23-24]。其中,GMS不僅可以模擬地下水中污染物質(zhì)的遷移轉(zhuǎn)化過程,還可以將包氣帶中污染物的垂向遷移考慮進(jìn)來,共同模擬[25-26]。

5 RS和 GIS技術(shù)在地下水農(nóng)業(yè)面源污染研究中的應(yīng)用

農(nóng)業(yè)面源污染是一種間歇發(fā)生的,隨機(jī)性、不確定性很強(qiáng)的復(fù)雜過程,估算其發(fā)生負(fù)荷的難度大。RS可適時(shí)獲取對(duì)地觀測(cè)的空間信息,GIS具有強(qiáng)大的空間數(shù)據(jù)管理、分析和制圖能力。因此,它們可以廣泛應(yīng)用于具有空間特征的農(nóng)業(yè)面源污染研究。RS和 GIS在面源污染研究中的應(yīng)用本質(zhì)上是 RS及 GIS與面源模型的結(jié)合問題。而具體的應(yīng)用主要有以下幾個(gè)方面。

5.1 農(nóng)業(yè)面源污染中地下水質(zhì)的評(píng)價(jià)與預(yù)測(cè)

非點(diǎn)源污染具有廣泛性、隨機(jī)性、潛在威脅性及空間相關(guān)性的特點(diǎn),對(duì)其進(jìn)行分析和控制較點(diǎn)源污染更為復(fù)雜。將GIS與專業(yè)模型有機(jī)結(jié)合,通過空間分析和統(tǒng)計(jì)為環(huán)境模型提供基礎(chǔ)數(shù)據(jù)與參數(shù),進(jìn)行水質(zhì)評(píng)價(jià)與預(yù)測(cè),以輔助水污染控制規(guī)劃決策。楊悅鎖等[27]將 G IS與 DRASTIC方法結(jié)合開發(fā)了一個(gè)動(dòng)態(tài)風(fēng)險(xiǎn)評(píng)價(jià)方法,并將其運(yùn)用于英國北愛爾蘭Upper Bann流域中的一個(gè)小流域。此方法也有能效地幫助決策者在流域范圍內(nèi)開展農(nóng)業(yè)面源地下水污染預(yù)防措施。

5.2 農(nóng)業(yè)面源污染的時(shí)空變化的研究和預(yù)測(cè)

GIS能存貯、處理、顯示大量與非點(diǎn)源污染有關(guān)的數(shù)據(jù)資料,RS技術(shù)是一種測(cè)定溶質(zhì)傳輸模型所需要的物理、化學(xué)、生物特性及其它數(shù)據(jù)資料最經(jīng)濟(jì)有效的方法。非點(diǎn)源污染模型將這些方法有機(jī)地結(jié)合來模擬非點(diǎn)源污染物的時(shí)間與空間變化。蔣勇軍等[28]以云南省瀘西縣小江典型巖溶農(nóng)業(yè)流域?yàn)檠芯繂卧?利用 1982年和 2004年的地下水質(zhì)數(shù)據(jù)及1982年的航片和 2004年的 TM影像,在 GIS支持下,研究其20年來的地下水質(zhì)的時(shí)空變化及原因。

5.3 地下水體污染物的空間分布、含量及其變化的模擬

將地下水模型和 GIS技術(shù)結(jié)合可以有效的模擬水流和水體中面源污染物的變化。王明新等[29]針對(duì)于面源污染機(jī)理模型在實(shí)際運(yùn)用中的限制,將人工神經(jīng)網(wǎng)絡(luò)引入地下水非點(diǎn)源污染格局的模擬和預(yù)報(bào)中,建立了基于 GIS的 BP神經(jīng)網(wǎng)絡(luò)模型用以模擬分析農(nóng)區(qū)淺層地下水 NO-3-N含量及其空間分布特征。

毫無疑問,RS與 GIS正推動(dòng)溶質(zhì)資料搜集、傳輸模擬、空間數(shù)據(jù)庫三者有機(jī)結(jié)合,以評(píng)價(jià)地下水面污染問題。然而值得注意的是,利用 GIS模擬的結(jié)果也不能完全取代田間試驗(yàn)和調(diào)查。

6 地下水農(nóng)業(yè)面源污染的防治對(duì)策

6.1 加強(qiáng)管理,依法治污

加強(qiáng)對(duì)農(nóng)藥、化肥生產(chǎn)管理,切實(shí)減少對(duì)環(huán)境有長期影響的有害物質(zhì)的使用量,嚴(yán)格按照國家環(huán)境保護(hù)部制訂的有關(guān)法規(guī)控制規(guī)模化畜禽養(yǎng)殖場糞污排放并實(shí)行嚴(yán)格監(jiān)管制度。對(duì)確實(shí)超標(biāo)排放的單位可加大處罰力度。加大輿論宣傳力度,提高人們特別是廣大決策者和農(nóng)民對(duì)面源污染的認(rèn)識(shí),引導(dǎo)農(nóng)民科學(xué)施肥、合理使用農(nóng)藥等,盡量減少由于農(nóng)事活動(dòng)的不科學(xué)而造成資源浪費(fèi)和環(huán)境中殘余污染物的增加。健全的農(nóng)產(chǎn)品監(jiān)測(cè)和市場準(zhǔn)入制度及農(nóng)產(chǎn)品認(rèn)證制度有助于有效防治農(nóng)業(yè)面源污染。

6.2 合理施肥,科學(xué)用藥

我國農(nóng)田化肥、農(nóng)藥施用量普遍偏高特別是氮素化肥施用量;磷肥施用也處于過量水平而鉀肥嚴(yán)重不足。研究表明：農(nóng)田投入養(yǎng)分過大,盈余部分并未起作用而是最終進(jìn)入土壤和水環(huán)境造成地表水和地下水環(huán)境的污染。因而防治重點(diǎn)應(yīng)放在減量、提質(zhì)、增效上,將“末端治理”改為“源頭治理”,研制高效、多功能復(fù)合肥,并從技術(shù)上指導(dǎo)農(nóng)民嚴(yán)格控制氮肥的使用量,平衡氮肥、磷肥、鉀肥的施用比例,減少流失量。同時(shí)通過增加秸稈還田、增施有機(jī)肥來增加鉀肥。鑒于目前農(nóng)藥的不可替代性,提倡農(nóng)業(yè)綜合防治策略。包括利用耕作、栽培、育種等農(nóng)事措施來防治農(nóng)作物病蟲害;利用生物技術(shù)和基因技術(shù)防治農(nóng)業(yè)有害生物;應(yīng)用光、電、微波、超聲波、輻射等物理措施來控制病蟲害;減量使用;推廣廣譜、高效環(huán)保型的生物農(nóng)藥[30]。

6.3 科學(xué)布局,合理規(guī)劃,發(fā)展復(fù)合農(nóng)業(yè)

農(nóng)業(yè)面源污染主要來自于農(nóng)村的農(nóng)業(yè)生產(chǎn)活動(dòng),是由不合理的耕作制度、粗放的施肥措施、農(nóng)業(yè)模式的不合理組合、畜禽養(yǎng)殖的不規(guī)范化和農(nóng)村宣傳力度不夠等因素造成的。復(fù)合農(nóng)業(yè)生態(tài)系統(tǒng)則強(qiáng)化農(nóng)業(yè)模式的合理優(yōu)化組合,在加強(qiáng)農(nóng)民生產(chǎn)生活知識(shí)的同時(shí),大力宣傳和認(rèn)識(shí)復(fù)合農(nóng)業(yè),合理利用有限的土地資源,創(chuàng)造更大的經(jīng)濟(jì)效益;鼓勵(lì)人們發(fā)展生態(tài)農(nóng)業(yè),減少養(yǎng)分流失,保證下游地區(qū)的農(nóng)業(yè)生產(chǎn)和生活安全;從源頭減少污染物的排放量,使其達(dá)到國家排放標(biāo)準(zhǔn)。這些措施為最終治理地下水農(nóng)業(yè)面源污染打下良好的基礎(chǔ)。

7 展望

(1)我國地下水農(nóng)業(yè)面源污染的問題很突出,但目前的管理方式還集中在肥料管理上,特別是還沒有污染源等級(jí)劃分標(biāo)準(zhǔn),這將導(dǎo)致管理上的盲目性。我國幅員遼闊,產(chǎn)生面源污染的時(shí)空變化情況比較復(fù)雜,影響因素較多,因此,在等級(jí)標(biāo)準(zhǔn)的確定、數(shù)學(xué)模型的建立等方面還需進(jìn)行大量細(xì)致而深入的工作[31]。

(2)現(xiàn)今,科學(xué)研究工作者仍在努力嘗試并找出解決地下水農(nóng)業(yè)面源污染的經(jīng)濟(jì)、合理且易于推廣實(shí)施的生態(tài)防治措施,但其效果并不明顯。因此,在未來很長一段時(shí)間內(nèi),地下水農(nóng)業(yè)面源污染的治理仍將是科研工作者的重中之重。通過對(duì)復(fù)合農(nóng)業(yè)生態(tài)系統(tǒng)原理、方法及機(jī)理的進(jìn)一步研究,總結(jié)合理的農(nóng)業(yè)管理模式,加快其推廣應(yīng)用,以復(fù)合農(nóng)業(yè)生態(tài)系統(tǒng)治理地下水農(nóng)業(yè)面源污染將是主要趨勢(shì)之一[29]。

(3)利用 “3S”集成技術(shù) (GIS、GPS、RS),針對(duì)地下水農(nóng)業(yè)面源污染的空間特性,可有效地監(jiān)控和治理地下水農(nóng)業(yè)面源污染,為開展全國范圍的面源污染調(diào)查提供有效的工具,為制定控制策略提供全面可靠的信息,為各級(jí)政府在發(fā)展規(guī)劃中農(nóng)業(yè)面源污染評(píng)價(jià)指標(biāo)體系的建立提供有效的依據(jù)。3S集成技術(shù)也將更多應(yīng)用于地下水農(nóng)業(yè)面源污染防治的研究中[32]。

如果說,農(nóng)業(yè)面源污染控制的成敗直接關(guān)系到我國流域污染治理的成敗,那么,地下水污染防治的成敗則關(guān)系到我國農(nóng)業(yè)面源污染控制的成敗。與地表水相比,地下水流速慢、循環(huán)周期長,一旦被污染,則很難凈化,這大大增加了地下水污染防治的難度。故從策略上來講,更應(yīng)防患于未然,極早防治地下水農(nóng)業(yè)面源污染,否則會(huì)事倍功半。地下水農(nóng)業(yè)面源污染是危害全球生態(tài)環(huán)境的一個(gè)最重要因素。應(yīng)借鑒其他發(fā)達(dá)國家的經(jīng)驗(yàn),根據(jù)控制區(qū)域自然、經(jīng)濟(jì)、技術(shù)條件及地下水面源污染特征的不同,因地制宜地制定多種控制與管理措施,最終實(shí)現(xiàn)環(huán)境、經(jīng)濟(jì)和社會(huì)的協(xié)調(diào)發(fā)展。

[1]董少杰,董哲等.大沽河中下游地區(qū)地下水環(huán)境特征與污染機(jī)理研究[J].水資源與水工程學(xué)報(bào),2008,19(6)：111-114.

[2]劉光棟,吳文良等.華北農(nóng)業(yè)高產(chǎn)糧區(qū)地下水面源污染特征及環(huán)境影響研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2005,13(2)：125-128

[3]張維理,武淑霞等.中國農(nóng)業(yè)面源污染形勢(shì)估計(jì)及控制對(duì)策[J].中國農(nóng)業(yè)科學(xué),2004,37(7)：1008-1017

[4]薛禹群,張幼寬.地下水污染防治在我國水體污染控制與治理中的雙重意義[J].環(huán)境科學(xué)學(xué)報(bào),2009,29(3)：474-481

[5]黃東風(fēng),王果.菜地土壤氮磷面源污染現(xiàn)狀、機(jī)制及控制技術(shù)[J].應(yīng)用生態(tài)學(xué)報(bào),2009,20(4)：991-1001

[6]王建兵,程磊.農(nóng)業(yè)面源污染現(xiàn)狀分析[J].江西農(nóng)業(yè)大學(xué)學(xué)報(bào),2008,(3)：35-38.

[7]劉僑博,劉薇等.非點(diǎn)源污染影響分析及防治措施[J].環(huán)境科學(xué)與管理 2010,35(6)：106-110

[8]王道涵梁成華.農(nóng)業(yè)磷素流失途徑及控制方法研究進(jìn)展[J].土壤與環(huán)境,2002,11(2)：183-188

[9]張杰麗,林濤.農(nóng)村面源污染治理建議與對(duì)策[J].遼寧城鄉(xiāng)環(huán)境科技,2006,26(4)：11-12

[10]張青松,劉飛等.農(nóng)業(yè)化肥面源污染現(xiàn)狀及對(duì)策[J].亞熱帶水土保持,2010,22(2)：44-45,52

[11]周祖澄.固體氮肥施入旱田土壤中去向的研究[J].環(huán)境科學(xué),1985,(6)：2-7

[12]李慶逵.現(xiàn)代磷肥研究的進(jìn)展[J].土壤學(xué)進(jìn)展,1986,(2)：1-7

[13]賀纏生,傅伯杰等.非點(diǎn)源污染的管理及控制[J].環(huán)境科學(xué),1998(9)：87-98.

[14]郭淑敏,程序等.京郊農(nóng)業(yè)活動(dòng)非點(diǎn)源污染現(xiàn)狀及防治對(duì)策[J].Environmental Protection,2004,(4)：32-36

[15]陳建耀,王亞等.地下水硝酸鹽污染研究綜述[J].地理科學(xué)進(jìn)展,2006,25(1)：34-44

[16]黃紅艷,高揚(yáng)等.上海市都市農(nóng)業(yè)區(qū)域地下水磷素非點(diǎn)源污染特征研究[J].水土保持學(xué)報(bào),2010,24(1)：101-113

[17]虞子婧,王明玉等.農(nóng)業(yè)生態(tài)系統(tǒng)地下水氮面源污染數(shù)值模擬[J].中國科學(xué)院研究生院學(xué)報(bào),2007,(9)：774-780

[18]張秋玲,陳英旭等.非點(diǎn)源污染模型研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào),2007,18(8)：1886-1890.

[19]湯亞飛,王焰新等.農(nóng)藥在包氣帶中遷移過程的數(shù)值模擬[J].安全與環(huán)境工程,2005,12(2)：11-14.

[20]王吉蘋,曹文志.應(yīng)用GLEAMS模型評(píng)估我國東南地區(qū)農(nóng)業(yè)小流域硝態(tài)氮的滲漏淋失[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2007,23(1)：28-32.

[21]畢經(jīng)偉,張佳寶等.應(yīng)用 HYDRUS-1D模型模擬農(nóng)田土壤水滲漏及硝態(tài)氮淋失特征[J].農(nóng)村生態(tài)環(huán)境,2004,20(2)：28-32.

[22]張洪霞,宋文.地下水?dāng)?shù)值模擬的研究現(xiàn)狀與展望[J].水利科技與經(jīng)濟(jì),2007,13(11)：794-796.

[23]徐艷杰,常利武等.FEFLOW在地下水?dāng)?shù)值模擬中的應(yīng)用[J].華北水利水電學(xué)院學(xué)報(bào),2009,30(2)：86-88.

[24]王慶永,賈忠華等.V isual MODFLOW及其在地下水模擬中的應(yīng)用[J].水資源與水工程學(xué)報(bào),2007,18(5)：90-92.

[25]曹淵,王鐵良等.GMS在飽和—非飽和地下水流動(dòng)及核素遷移模擬中的應(yīng)用[J].遼寧工程技術(shù)大學(xué)學(xué)報(bào)(自然科學(xué)版),2009,(28)：190-192.

[26]張倩,蘇保林等.農(nóng)田面源氮素垂向運(yùn)移負(fù)荷估算方法綜述[J].南水北調(diào)與水利科技 2010,8(3)：35-38

[27]Y.S.Yang,J.L.Wang.GIS-based Dynamic Risk Assessment for Groundwater Nitrate Pollution from Agricultural Diffuse Sources[J].Journal of Jilin University(Earth Science Edition),2007,37(2)：311-318

[28]蔣勇軍,袁道先等.典型巖溶農(nóng)業(yè)區(qū)地下水質(zhì)與土地利用變化分析[J].地理學(xué)報(bào) 2006,61(5)：471-481

[29]王明新,吳文良等.基于 G IS和 BP神經(jīng)網(wǎng)絡(luò)的農(nóng)區(qū)地下水硝態(tài)氮含量分布特征分析[J].農(nóng)業(yè)工程學(xué)報(bào),2006,22(12)：39-43

[30]王道涵,梁成華.農(nóng)業(yè)磷素流失途徑及控制方法研究進(jìn)展[J].土壤與環(huán)境,2002,11(2)：183-188

[31]唐蓮,張衛(wèi)兵.寧夏水環(huán)境保護(hù)與農(nóng)業(yè)非點(diǎn)源污染[J].水土保持研究,2007,14(5)：239-241

[32]A.P.Karen,L.B.Barbara.GIS-based non–point source pollution modeling：Considerations for wetlands[J],Journal of Soil and W ater Conservation,1995,(6)：613-619