田皓予, 佟 玲, 余國安, Arika Bridhikitti

(1.中國農(nóng)業(yè)大學(xué) 水利與土木工程學(xué)院, 北京 100083; 2.中國科學(xué)院 地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點實驗室, 北京100101; 3.瑪希隆大學(xué) 北碧府校區(qū)跨學(xué)科研究學(xué)院 環(huán)境工程與災(zāi)害管理系, 泰國 71150)

1 研究背景

隨著經(jīng)濟社會快速發(fā)展，地表水質(zhì)退化作為全球環(huán)境問題之一，日益受到高度關(guān)注。土地利用是影響流域水質(zhì)變化的重要人為驅(qū)動力。土地利用變化不僅影響流域水文過程[1]，還將引起水體營養(yǎng)鹽、細菌、有毒有害物含量、溶解氧變化等化學(xué)響應(yīng)[2]。楊琴等[3]通過Spearman相關(guān)分析發(fā)現(xiàn)淮河流域的水質(zhì)情況與土地利用類型百分比有關(guān)。農(nóng)業(yè)和城市用地的增加被認為是全球淡水生態(tài)系統(tǒng)營養(yǎng)和沉積物增加的最大貢獻者之一[4-5]。養(yǎng)分(P, N)的過度使用會影響河流水質(zhì)。Turner等[6]研究表明，密集而廣泛的耕作，特別是化肥的使用，對營養(yǎng)負荷有顯著的影響。農(nóng)業(yè)非點源污染導(dǎo)致氮磷營養(yǎng)物、農(nóng)藥等流入河流，引起水體富營養(yǎng)化，毒害物質(zhì)含量增加，破壞水體環(huán)境，進一步影響人類生產(chǎn)生活甚至威脅人體健康[7]。

農(nóng)業(yè)種植、城市化進程和人口迅速增長是影響泰國地表水質(zhì)的主要原因。泰國污染管制署(Pollution Control Department, Thailand,簡稱PCD)自1994年起監(jiān)測全國48條主要河流的水質(zhì)，并將觀測指標(biāo)作為管理河流水環(huán)境的基線[8-11]。

水質(zhì)評價有很長的歷史，早期的評價主要是定性描述水的感官特性，如顏色、氣味和渾濁度。隨著科學(xué)技術(shù)的不斷發(fā)展以及研究的深入，對水體各項理化性質(zhì)的認識更為深刻，進而出現(xiàn)了不同的水質(zhì)評價方法。如單因子評價法[12]、水污染指數(shù)法[13-14]、水質(zhì)標(biāo)識指數(shù)法[15]、模糊評價法[16-17]、主成分分析[18-19]和人工神經(jīng)網(wǎng)絡(luò)法[20-22]等。

水質(zhì)指數(shù)(Water quality index，簡稱WQI)最初由Horton于1965年提出，1972年Brown小組開發(fā)了一種類似于Horton指數(shù)的新的WQI[23]。后來許多科學(xué)家和專家對WQI進行了修改(如美國國家衛(wèi)生基金會水質(zhì)指數(shù)(NSFWQI)[24]、加拿大環(huán)境質(zhì)量部長理事會水質(zhì)指數(shù)(CCMEWQI)[25]、俄勒岡水質(zhì)指數(shù)(OWQI)[26])。WQI的確定有一個標(biāo)準的步驟，首先由行業(yè)專家、相關(guān)機構(gòu)或政府部門從對水質(zhì)有較大影響的因素(如氧含量、富營養(yǎng)化、健康方面、物理特性和溶解物質(zhì))中選取代表性指標(biāo)；然后確定各指標(biāo)的適合度曲線；最后根據(jù)各指標(biāo)對當(dāng)?shù)厮|(zhì)影響程度進行加權(quán)聚合得到最終的WQI值。WQI法將大量的數(shù)據(jù)轉(zhuǎn)化為一個有效反映整體水質(zhì)狀況的數(shù)值，為管理者和決策者提供了一種簡單高效的了解水質(zhì)狀況的方法[27-28]。

在“一帶一路”大背景下，大湄公河次區(qū)域(GMS)水環(huán)境管理對話合作正成為“瀾湄合作”的重要內(nèi)容，其中，區(qū)域水環(huán)境分析評估作為水環(huán)境管理的基礎(chǔ)日益受到關(guān)注。本研究結(jié)合蒙河干流1997-2017年水環(huán)境指標(biāo)持續(xù)監(jiān)測結(jié)果，分析蒙河流域代表性水質(zhì)指標(biāo)和整體水環(huán)境的時空分異特征，了解在人類活動影響下，水環(huán)境變化是否呈現(xiàn)趨勢性特征，并嘗試分析影響水環(huán)境時空變化的潛在原因，為改善蒙河流域水環(huán)境提供科學(xué)依據(jù)，并希冀對自然條件相近且以農(nóng)業(yè)種植為主的大湄公河次區(qū)域的其他中小流域水環(huán)境管理有所借鑒。

2 研究區(qū)域概況

蒙河是湄公河最大支流，位于泰國東北部，大致位于北緯14°～16°，東經(jīng)101° 30′～105° 30′之間，干流河長約800 km，流域面積約8.2×104km2，蒙河流域及其水質(zhì)、水文站點分布見圖1。

圖1 蒙河流域及其水質(zhì)監(jiān)測和水文站點分布圖

蒙河流域?qū)賮啛釒駶櫄夂?，受亞熱帶季風(fēng)影響，流域干濕季氣候差異分明，其河流水文過程(汛期和非汛期)的季節(jié)差異顯著。流域內(nèi)年平均溫度不低于18℃，年降雨量為1 300～1 500 mm，降雨一般在5月中旬-10月中旬，主要降雨過程在8-9月，10月下旬至5月中旬降雨稀少，且2月中旬至5月中旬較為炎熱[29-30]。流域主要土壤類型為沙質(zhì)壤土、壤質(zhì)砂土等，結(jié)構(gòu)較為松散，抗蝕力較低。這里是泰國重要的農(nóng)業(yè)區(qū)，主要種植作物包括水稻、木薯、甘蔗等。

3 數(shù)據(jù)與方法

3.1 數(shù)據(jù)來源

泰國污染管制署(PCD)自20世紀90年代中期開始在蒙河流域布置水質(zhì)站點，開展水質(zhì)監(jiān)測。目前，蒙河干流沿河分布21個水質(zhì)監(jiān)測站點(圖1)。正常情況下，PCD每年在各監(jiān)測站點進行4次水質(zhì)采樣及測試，監(jiān)測時間一般分別為2月(有時在3月)、5月、8月和11月下旬(有時在12月份或1月上旬)，以反映和評估旱季(非汛期)和雨季(汛期)的河流水質(zhì)狀況。其中，8月份的水質(zhì)監(jiān)測反映汛期水質(zhì)狀況，其他3個時間段的監(jiān)測反映非汛期水質(zhì)狀況。本文基于PCD在1997-2017年間對蒙河水環(huán)境的持續(xù)監(jiān)測數(shù)據(jù)，選擇8個典型水質(zhì)指標(biāo)開展分析，包括溶解氧(DO)、生化需氧量(BOD)、總大腸菌群(TCB)、糞大腸菌群(FCB)、總磷(TP)、硝態(tài)氮(NO3—N)、氨氮(NH3—N)和懸浮物(SS)。

3.2 分析方法

季節(jié)性Kendall檢驗是基于Mann-Kendall檢驗并在此基礎(chǔ)上推廣的一種單調(diào)趨勢分析方法。該方法可以降低數(shù)據(jù)未檢出、漏測以及奇異值對水質(zhì)趨勢的影響，且不存在豐水期與非豐水期水質(zhì)資料不可比所帶來的影響，所以在水質(zhì)趨勢分析中具有廣泛的應(yīng)用[31]。在進行水質(zhì)趨勢分析時，數(shù)據(jù)序列的長短對檢驗結(jié)果有很大影響。過短的水質(zhì)序列在趨勢分析中的準確性得不到保證，而水質(zhì)序列過長可能會導(dǎo)致趨勢被掩蓋。針對各研究區(qū)及監(jiān)測數(shù)據(jù)的特點，不同學(xué)者選擇不同的時間序列長度。其中Chang[32]使用了8 a的水質(zhì)數(shù)據(jù)，Zhai等[33]使用了12 a的水質(zhì)監(jiān)測數(shù)據(jù)。本研究將1997-2017共21年的水質(zhì)監(jiān)測序列劃分為兩個時段(即1997-2007和 2008-2017)進行分析。該方法的原理是對于n年p月的水質(zhì)監(jiān)測資料，將每年相同月份的數(shù)據(jù)(x11，…，xnp)進行比較，計算比較結(jié)果之和S:

(1)

p月內(nèi)可作比較的差值數(shù)據(jù)組個數(shù)m為：

(2)

在零假設(shè)下，p月S的均值為：E(S)=0

當(dāng)ni個水質(zhì)序列非漏測值中有t個數(shù)相同時，方差為：

(3)

當(dāng)n>10時，S同樣服從正態(tài)分布，其標(biāo)準方差Z為：

(4)

Kendall檢驗統(tǒng)計量τ為：

(5)

若τ為正說明水質(zhì)序列具有上升趨勢，τ為負說明具有下降趨勢，τ為零說明無變化趨勢[31]。一般情況下，取顯著性水平P為0.05和0.01，相應(yīng)的標(biāo)準方差Z值為1.96和2.58。即|Z|大于2.58說明具有高度顯著變化趨勢，介于1.96和2.58之間說明具有顯著變化趨勢，小于1. 96無顯著性。

利用PCD于2006年對代表性水質(zhì)指標(biāo)開發(fā)的適合度曲線，將各水質(zhì)指標(biāo)的觀測平均值轉(zhuǎn)換為子指標(biāo)得分。為了得到適合度曲線，PCD向泰國約100名水質(zhì)專家發(fā)放了問卷，要求他們獨立選擇重要水質(zhì)指標(biāo)，并對所選指標(biāo)的取值范圍進行評分(從最差的0分到最好的100分)。PCD利用這些建議為每個指標(biāo)繪制適合度曲線。為更好地評估泰國河流的水質(zhì)狀況，使用泰國45條河流的746個采樣站點的相關(guān)水質(zhì)指標(biāo)觀測值進行為期一年的測試，然后做出相關(guān)調(diào)整，得到了最終的適合度曲線的評分公式[34-36]。

本文采用非加權(quán)調(diào)和平方平均WQI的方法評估流域整體水環(huán)境的優(yōu)劣。該方法被認為是對算術(shù)加權(quán)、幾何加權(quán)以及非加權(quán)平均方法的改進[37]，其特點是側(cè)重考慮超標(biāo)嚴重的水質(zhì)指標(biāo)對整體水質(zhì)評價結(jié)果的影響，并且認為不同的水質(zhì)變量在不同的時間和地點對整體水質(zhì)的影響是不同的[26]，其表達式如下：

(6)

式中：SIi為各水質(zhì)指標(biāo)對應(yīng)的適合度曲線轉(zhuǎn)換后的水質(zhì)指標(biāo)i的得分;i為研究所選的水質(zhì)指標(biāo)(本文為DO、BOD、TCB、FCB、TP、NO3—N、NH3—N、SS);m為所選水質(zhì)指標(biāo)的總數(shù)(本研究中m=8)?；诟鱾€站點的WQI計算值，參考泰國國家環(huán)境局(National Environment Board, Thailand)[38]劃分的5類地表水特征及相應(yīng)的WQI值范圍(表1)可對該站點的地表水水質(zhì)進行分類。

考慮到水質(zhì)監(jiān)測站點MU08與MU09之間有大型支流齊河(Chi River)匯入，本文將MU09以上河段(MU09-MU20)作為一個整體，稱為上游，MU08以下(MU01-MU08)為一個整體，稱為下游。

4 結(jié)果與分析

4.1 主要水環(huán)境指標(biāo)時空分異

對1997-2017年蒙河流域的21個水質(zhì)監(jiān)測站點的數(shù)據(jù)分析表明，主要水質(zhì)指標(biāo)含量具有顯著的時空分異特征。選取4個典型監(jiān)測站點(MU17、MU12、MU06、MU01)的代表性水質(zhì)指標(biāo)(DO、BOD、TP、NO3—N)進行分析，結(jié)果見圖2。由圖2可看出，各指標(biāo)含量在1997-2017年雖無明顯增加或減小趨勢，但隨時間變化波動劇烈，且空間差異明顯?；旧希嫌蜠O含量較下游低，而BOD、TP含量較下游高。按照泰國地表水質(zhì)標(biāo)準[38](泰國地表水質(zhì)標(biāo)準中無TP分類標(biāo)準，參考中國《地表水環(huán)境質(zhì)量標(biāo)準(GB3838-2002)》中TP分類標(biāo)準)，上游DO含量主要為3～4級，下游為2～3級；上游BOD含量主要為2～4級，下游為2～3級；上游TP含量主要為3～5級，下游為1～3級；整個河段NO3—N含量主要為3級及優(yōu)于3級。

表1 泰國5類地表水特征及相應(yīng)WQI值范圍

圖2 1997-2017年蒙河4個監(jiān)測站點代表性水質(zhì)指標(biāo)隨時間變化

蒙河流域?qū)儆趤啛釒駶櫄夂?，該地區(qū)汛期和非汛期降雨和徑流差異十分顯著，而汛期和非汛期河流水文過程(流量、流速)差異對河流水質(zhì)有很大影響。利用季節(jié)性Kendall檢驗得到1997-2007、2008-2017年兩個時段18個水質(zhì)站點代表性水質(zhì)指標(biāo)的時間變化趨勢，計算結(jié)果見表2、3。

表2、3中季節(jié)性Kendall檢驗結(jié)果顯示，兩個分析時段BOD、TCB、FCB、SS含量變化趨勢不同。1997-2007年TCB、FCB含量呈顯著上升趨勢，而近10年(2008-2017)呈顯著下降趨勢。21世紀初水體的主要污染來自糞大腸菌群的污染，高生化需氧量導(dǎo)致低溶解氧和高氨水平，近年來相關(guān)環(huán)境保護部門大力控制水體中細菌類污染物含量，直到2015年蒙河水體TCB達標(biāo)率為96%，F(xiàn)CB達標(biāo)率為99%。而SS含量由1997-2007年的幾乎無顯著變化到2008-2017年大部分水質(zhì)站點呈下降趨勢(特別是上游)反映了土地利用變化對于SS含量的影響。

蒙河上游與下游DO、BOD含量的趨勢性變化具有明顯的差異，對于大者為優(yōu)的收益性指標(biāo)DO，上游河段呈上升趨勢，且2008-2017時段較1997-2007時段的變化更為明顯。

1997-2007時段下游無顯著變化趨勢，但2008-2017時段大部分站點呈下降趨勢。對于小者為優(yōu)的成本性指標(biāo)BOD，上游除個別站點(MU16，MU17)在2008-2017時段有下降趨勢外，大部分站點的BOD含量在兩個分析時段(1997-2007；2008-2017)未監(jiān)測到顯著趨勢性變化，但下游多個站點監(jiān)測結(jié)果表明BOD含量呈上升趨勢。說明齊河(Chi River)的匯入對于蒙河水體中的溶解氧含量和生化需氧量有很大的影響。

對于氮類營養(yǎng)物含量(NO3—N、NH3—N)流域內(nèi)大部分站點呈現(xiàn)顯著上升趨勢，特別是2008-2017年流域下游所有站點均呈顯著上升趨勢，可能是由于近些年來農(nóng)藥化肥的施用量增多導(dǎo)致蒙河下游水體中氮類營養(yǎng)物質(zhì)含量呈現(xiàn)上升趨勢，且由于人口增多，排入河流中的生活污水增加也會導(dǎo)致營養(yǎng)物質(zhì)含量的升高。比如下游MU06站點附近是蒙河下游最大的城市烏汶(Ubon Ratchathani)，近些年該城市規(guī)模擴大，人口增加，排入蒙河的生活污水量上升，而生活污水中通常含有大量的氮，該結(jié)果與流域?qū)嶋H情況吻合。

表2 1997-2007年18個站點代表性水質(zhì)指標(biāo)的季節(jié)性Kendall趨勢檢驗

注：↓和↑分別表示減小和增加趨勢，*和**分別表示通過0.05和0.01顯著性水平檢驗。

表3 2008-2017年18個站點代表性水質(zhì)指標(biāo)的季節(jié)性Kendall趨勢檢驗

注：↓和↑分別表示減小和增加趨勢，*和**分別表示通過0.05和0.01顯著性水平檢驗。

對1997-2017年4個代表性水質(zhì)指標(biāo)按汛期和非汛期、上游和下游進行分析對比，結(jié)果見圖3。由圖3可以看出，各代表性水質(zhì)指標(biāo)具有明顯的時空分異特征。蒙河汛期的DO濃度低于非汛期，而BOD、TP和NO3—N濃度高于非汛期，這與其他河流汛期水質(zhì)一般好于非汛期水質(zhì)的情況有所不同，可能與蒙河流域主要為農(nóng)業(yè)種植區(qū)且降雨集中有關(guān)，蒙河的沙質(zhì)壤土和壤質(zhì)砂土在汛期強降雨時易受侵蝕，土壤養(yǎng)分淋溶導(dǎo)致汛期營養(yǎng)物質(zhì)進入水體。在空間分布上，蒙河上游DO濃度低于下游，而BOD、TP濃度高于下游。

為進一步探究代表性水質(zhì)指標(biāo)的空間變化，對1997-2017年蒙河干流河段21個水質(zhì)監(jiān)測站點監(jiān)測數(shù)據(jù)進行分析，結(jié)果見圖4。由圖4上看，流域水體的DO濃度呈下游向上游遞減趨勢，而BOD、TP濃度則呈下游向上游遞增趨勢，且TP濃度在MU08以下和MU09以上差異明顯，可能是由于齊河(Chi River)匯入的影響。NO3—N的空間變化則具有季節(jié)性差異，在汛期一般上游河段含量較高，而在非汛期通常下游河段含量較高?？赡茉蚴巧嫌沃饕獮檗r(nóng)業(yè)種植區(qū)，汛期正值莊稼生長季，化肥施用成為NO3—N濃度高的主要因素，而下游城市鄉(xiāng)鎮(zhèn)較多，非汛期生活污水的排入成為河流中NO3—N濃度高的主要原因。

4.2 整體水環(huán)境時空變化

4個典型站點(MU01、MU06；MU12、MU17)WQI值計算結(jié)果顯示，1997-2017年蒙河WQI值波動劇烈。不過，僅MU01站點的WQI值年際變化在0.01的顯著性水平上有減小趨勢，其他3個站點沒有顯著的趨勢性變化。4個站點的WQI值在年內(nèi)均表現(xiàn)出顯著的季節(jié)性差異，與非汛期相比，汛期的WQI值通常較低(圖5)，說明蒙河流域汛期的水質(zhì)較非汛期差。通過五類地表水對應(yīng)的WQI值范圍(表1)可以看出非汛期蒙河流域水質(zhì)基本處于Ⅲ類或優(yōu)于Ⅲ類(特別是下游)。

圖3 1997-2017年蒙河代表性水質(zhì)指標(biāo)濃度上下游、汛期非汛期對比

注：圖中橫坐標(biāo)1～20分別代表水質(zhì)站點MU01～MU20圖4 1997-2017年蒙河干流各水質(zhì)站點代表性水質(zhì)指標(biāo)汛期和非汛期平均濃度

圖5 1997-2017年蒙河干流4個典型站點汛期非汛期WQI值變化對比

與單項水質(zhì)指標(biāo)相似，1997-2017年蒙河不同站點多年平均WQI值也具有顯著的空間差異(圖6)。蒙河下游WQI值通常較高，說明下游水質(zhì)好于上游。

圖6 1997-2017年蒙河多年平均WQI值空間變化

4.3 蒙河水環(huán)境時空分異特征形成原因初步分析

1997-2017年蒙河下游NO3—N含量明顯增加，這可能與蒙河的水文氣象條件、土壤物理屬性以及下游城鎮(zhèn)規(guī)模擴大有關(guān)，隨著人口不斷增多，生活污水成為蒙河NO3—N含量增加的主要原因。

蒙河流域年降水量變化較大，最大年降水量在1549.9～2194.0 mm之間，最小年降水量在604.9～999.4 mm之間，汛期降水量約占全年降水總量的92%～98%。且降水分布不均勻，年降水量空間分布表現(xiàn)為由上游向下游梯次增加。圖7展示了蒙河上游春蓬(Chumphon)和下游烏汶(Ubon)兩個水文站1980-2014年汛期(6-10月)和非汛期(11-次年5月)平均流量的變化。由圖7可知，蒙河流域徑流過程波動劇烈，但并沒有顯著的趨勢性變化。徑流量變化過程存在明顯的年內(nèi)季節(jié)性變化和空間差異。春蓬站多年汛期平均流量為99.62 m3/s，是非汛期42.17 m3/s的2.4倍，烏汶站多年汛期平均流量為1 140.80 m3/s，是非汛期282.40 m3/s的4倍。流域上下游徑流量差異顯著，同一年份同一時期烏汶站與春蓬站平均流量之比最高為33∶1，這主要是由于春蓬站位于上游海拔較高，烏汶站位于下游海拔較低且匯集了更多的支流，加之降雨的影響使得徑流量呈現(xiàn)該種時空分異特點。

圖7 1980-2014年春蓬、烏汶水文站汛期、非汛期月平均徑流量變化

一般情形下，在點源、面源入河排污量相對一定時，河流水質(zhì)的優(yōu)劣與徑流量相關(guān)性較強，其汛期水質(zhì)優(yōu)于非汛期，(如我國長江流域枯/平水期總磷污染較重，豐水期污染較輕[39])，因為汛期流量大，有助于稀釋和降低污染物濃度。但蒙河情況有所不同，其非汛期水質(zhì)往往優(yōu)于汛期，其原因一方面與流域所處地域的氣候環(huán)境有關(guān)，蒙河流域干濕季分明，全年降雨主要集中在雨季(8-10月)，非汛期降雨稀少且氣溫高、蒸發(fā)強烈，造成部分小支流、塘圩同干流的聯(lián)系暫時中斷，部分營養(yǎng)物(如磷、氮)蓄積在支流和塘圩中。進入雨季后，小支流和塘圩與干流重新聯(lián)通，蓄積的污染物隨支流匯入干流，導(dǎo)致干流汛期的營養(yǎng)物含量高于非汛期。雨季正值作物生長旺盛階段，農(nóng)藥化肥的施用也是導(dǎo)致水體營養(yǎng)物質(zhì)等污染物濃度含量高的原因。另一方面，蒙河流域的水質(zhì)變化也與該流域的土壤結(jié)構(gòu)有關(guān)，蒙河主要土壤類型包括沙質(zhì)壤土、砂土、壤質(zhì)砂土，通常透水能力較強而抗蝕能力較弱，汛期強降雨條件下，土體中細顆粒泥沙和營養(yǎng)物易被沖刷、淋溶進入水體，導(dǎo)致汛期水體中營養(yǎng)物和懸浮物含量高于非汛期。

5 結(jié) 論

本文以泰國東北部蒙河流域為例，基于季節(jié)性Kendall法和水質(zhì)指數(shù)(WQI)法分析了蒙河流域1997-2017年沿干流21個水質(zhì)監(jiān)測站點代表性水質(zhì)指標(biāo)的時空變化，并評估流域整體水環(huán)境時空分異特征。研究結(jié)果表明：

(1) 蒙河水質(zhì)總體尚好，大部分為Ⅲ類或優(yōu)于Ⅲ類。1997-2017年蒙河水質(zhì)在時間上波動劇烈，總體上無明顯趨勢性變化，但存在明顯的年內(nèi)季節(jié)性變化。與非汛期(2、5、11月)相比，汛期(8月)DO含量較低，BOD、TP、NO3—N含量較高，WQI值通常最低，說明汛期水質(zhì)情況較差。

(2) 蒙河代表性水質(zhì)指標(biāo)含量及整體水環(huán)境存在顯著的空間差異。近10年來上下游DO、BOD含量存在相反的變化趨勢，下游NO3—N、NH3—N含量呈現(xiàn)顯著的上升趨勢，但始終都低于地表水Ⅲ類標(biāo)準限值。與下游相比，上游DO含量較低，BOD含量較高，上游TP含量是下游的3～5倍。下游水體WQI值較高，水質(zhì)類別主要為Ⅱ類；上游WQI值較低，主要為Ⅲ類水，說明下游水質(zhì)好于上游。

(3) 蒙河流域是泰國重要的農(nóng)業(yè)區(qū)，農(nóng)業(yè)種植(尤其水稻種植)中化肥和農(nóng)藥的施用不可避免地影響流域水環(huán)境。由于蒙河流域上游有大面積的農(nóng)業(yè)種植區(qū)，其土壤侵蝕及營養(yǎng)物淋溶是造成上游水質(zhì)相對下游水質(zhì)較差的重要原因，應(yīng)控制化肥農(nóng)藥過度施用和城鎮(zhèn)生活污水直接排放。