東江流域非點(diǎn)源污染模擬及時(shí)空分布研究

劉 潔，陳曉宏，許振成，虢清偉，吳根義3，，王兆禮

(1.環(huán)境保護(hù)部 華南環(huán)境科學(xué)研究所，廣州510655;2.中山大學(xué) 水資源與環(huán)境研究中心，廣州510275;3.湖南農(nóng)業(yè)大學(xué)環(huán)境科學(xué)與工程系，長(zhǎng)沙410128;4.華南理工大學(xué) 土木與交通學(xué)院，廣州510641)

0 引言

隨著點(diǎn)源治理技術(shù)和管理手段的發(fā)展，非點(diǎn)源污染已成為重要的水環(huán)境污染類型之一。中國(guó)的非點(diǎn)源污染問(wèn)題也日益凸現(xiàn)，在三峽庫(kù)區(qū)以及太湖和滇池等重要湖泊，非點(diǎn)源污染已成為水質(zhì)惡化的主要原因之一［1-2］。2000年以來(lái)國(guó)內(nèi)對(duì)非點(diǎn)源污染的研究越來(lái)越受到重視，研究成果用于特定區(qū)域非點(diǎn)源污染的管理和控制取得了一定的效果。對(duì)于小尺度流域非點(diǎn)源污染的遷移轉(zhuǎn)化機(jī)理以及負(fù)荷計(jì)算等方面的研究成果頗豐［3-7］，但是在大中尺度流域非點(diǎn)源負(fù)荷模擬計(jì)算方法等方面的研究還需要進(jìn)一步的開(kāi)展，在研究氣候變化與土地利用格局對(duì)于大中尺度流域水環(huán)境的影響等方面還需要進(jìn)一步的深入。通過(guò)流域水文模型對(duì)整個(gè)流域系統(tǒng)內(nèi)部發(fā)生的復(fù)雜污染過(guò)程進(jìn)行定量描述，可幫助分析非點(diǎn)源污染產(chǎn)生的時(shí)間和空間特征，并進(jìn)一步評(píng)估氣候變化、土地利用的變化、不同的管理與技術(shù)措施等對(duì)非點(diǎn)源污染負(fù)荷和水質(zhì)的影響，為流域管理和規(guī)劃提供決策依據(jù)。

非點(diǎn)源污染的主要研究方法可以歸納為基于野外采樣布點(diǎn)的典型地塊研究［3-5］、數(shù)理統(tǒng)計(jì)分析［6-9］、水文模型模擬［10-13］。2010年后應(yīng)用國(guó)外成熟的分布式水文模型研究不同土地利用類型結(jié)構(gòu)、土地利用變化對(duì)水文水質(zhì)的影響成為研究熱點(diǎn)。從研究尺度來(lái)看，國(guó)內(nèi)的研究主要是典型小流域的水文和產(chǎn)污過(guò)程研究、不同小流域的土地利用類型對(duì)產(chǎn)流產(chǎn)污過(guò)程影響的對(duì)比研究，研究尺度較小。2010年后針對(duì)大中型流域的相關(guān)研究成果開(kāi)始更多地見(jiàn)諸于報(bào)道［12-13］。SWAT和L-THIA等分布式水文模型被應(yīng)用在我國(guó)三峽庫(kù)區(qū)流域［12］、四湖流域［13］等流域的研究中，分布式水文模型可模擬研究氣候和土地利用的長(zhǎng)期變化對(duì)流域水文水質(zhì)時(shí)空變異的影響，適用于空間差異較大的流域的研究。

總體上看，國(guó)內(nèi)對(duì)于大中尺度流域的非點(diǎn)源污染模擬與研究還較少，本研究應(yīng)用流域分布式水文模型HSPF(hydrological simulation program-fortran)，以位于中國(guó)南方濕潤(rùn)地區(qū)的東江流域?yàn)槔M(jìn)行中尺度流域的非點(diǎn)源污染的模擬研究，建立中尺度流域——東江流域的非點(diǎn)源污染模型，在此基礎(chǔ)上分析該流域非點(diǎn)源污染的時(shí)空分布，為流域非點(diǎn)源污染管理措施的制定與實(shí)施提供技術(shù)支持。研究過(guò)程和方法可為適用于中尺度流域研究的非點(diǎn)源污染模擬模型的建立以及揭示自然條件、人類活動(dòng)變化與非點(diǎn)源污染的統(tǒng)計(jì)關(guān)系提供參考。

1 研究區(qū)域概況

東江是珠江流域的主要河流之一，是廣州、惠州、東莞、深圳等地市的重要供水水源地，同時(shí)還擔(dān)負(fù)著對(duì)香港特別行政區(qū)的供水任務(wù)。東江發(fā)源于江西省尋烏縣椏髻缽，上游稱尋烏水，南流入廣東境內(nèi)，至龍川合河壩安遠(yuǎn)水后稱東江，干流由東北向西南流，河道長(zhǎng)度至石龍為520 km，至獅子洋為562 km，石龍以上流域總面積27 040 km2。

2 HSPF模型簡(jiǎn)介與基礎(chǔ)數(shù)據(jù)預(yù)處理

2.1 HSPF 模型簡(jiǎn)介

HSPF模型是由美國(guó)環(huán)境保護(hù)署開(kāi)發(fā)的水文模型，用于較大流域范圍內(nèi)自然和人工條件下的水系中水文水質(zhì)過(guò)程的連續(xù)模擬［14］。HSPF屬于半分布式水文模型，使用時(shí)先要依據(jù)河流水庫(kù)所控制的自然流域范圍將流域劃分為不同的子流域。HSPF水文響應(yīng)單元(HRU)的建立是基于將流域劃分為均質(zhì)的水文和水質(zhì)響應(yīng)單元，土壤層在每個(gè)HRU上被垂直劃分為3層(上儲(chǔ)水層、下儲(chǔ)水層和地下水儲(chǔ)水層)，每層徑流量由每個(gè)區(qū)域蓄水量的非線性函數(shù)求得。HSPF模型是較成熟的水文模型，在國(guó)外有廣泛的應(yīng)用，但在國(guó)內(nèi)應(yīng)用還較少。至今已發(fā)展到HSPF12.0，并內(nèi)嵌于美國(guó)環(huán)保署開(kāi)發(fā)的基于GIS技術(shù)的整合式平臺(tái)BASINS系統(tǒng)中。

2.2 數(shù)據(jù)預(yù)處理

2.2.1 氣象數(shù)據(jù)的預(yù)處理。降雨的空間分辨率的變化對(duì)徑流模擬結(jié)果有重要的影響［14-17］。為了降低氣象數(shù)據(jù)的空間分辨率對(duì)東江流域水文模型建模的影響，研究中采用了中國(guó)氣象數(shù)據(jù)共享服務(wù)網(wǎng)提供的中國(guó)地面降水日值0.5°×0.5°格點(diǎn)數(shù)據(jù)，作為氣象觀測(cè)站點(diǎn)數(shù)據(jù)的補(bǔ)充，數(shù)據(jù)經(jīng)交叉驗(yàn)證、誤差分析，質(zhì)量狀況良好。利用BASIN4.0系統(tǒng)中WDMUtil工具編輯生成模型運(yùn)行必需的東江流域氣象數(shù)據(jù)的時(shí)間序列文件(WDM文件)。

2.2.2 下墊面數(shù)據(jù)的預(yù)處理。對(duì)土地利用/覆被數(shù)據(jù)進(jìn)行加工處理，將東江流域土地利用類型統(tǒng)一為7種類型:耕地、園地、林地、草地、建設(shè)用地、水域和未利用地。為了充分考慮流域的空間異質(zhì)性，除考慮土地利用的差異性還需考慮土壤類型的差異性和坡度等下墊面因素。因此，將1980，2000年?yáng)|江流域的土地利用數(shù)據(jù)與土壤分布、坡度數(shù)據(jù)分別進(jìn)行疊置得到新的下墊面類型數(shù)據(jù)。每種不同的下墊面類型對(duì)應(yīng)不同的相關(guān)水文參數(shù)。

2.2.3 子流域和水文響應(yīng)單元的劃分。HSPF的空間異質(zhì)性主要考慮每個(gè)離散單元的降雨、土壤類型、土地利用類型等。土壤類型、土地利用類型、坡度的空間異質(zhì)性通過(guò)疊加下墊面類型數(shù)據(jù)體現(xiàn)，氣象數(shù)據(jù)的空間異質(zhì)性則通過(guò)流域再分割來(lái)實(shí)現(xiàn)。在子流域劃分的基礎(chǔ)上，根據(jù)氣象數(shù)據(jù)的空間分布，進(jìn)行水文響應(yīng)單元的分割處理。東江流域水文響應(yīng)單元的氣象站點(diǎn)分配根據(jù)泰森多邊形方法進(jìn)行分割。

3 非點(diǎn)源污染模擬的HSPF模型建立

3.1 模型模擬效果評(píng)估的基本方法

采用標(biāo)準(zhǔn)偏差法(Dv)［18］與NASH-Sutcliff相關(guān)系數(shù)法(ENS)［19］評(píng)估模型的模擬效果。理想情況下，模型中的ENS接近于1，Dv接近于0。根據(jù)已有的標(biāo)準(zhǔn)［20-21］，認(rèn)為Dv在觀測(cè)值的10%范圍內(nèi)的模擬已經(jīng)達(dá)到了“非常好”的標(biāo)準(zhǔn)。NASH系數(shù)取值在0～1之間，并認(rèn)為0.9以上為甲等，0.7 ～ 0.9 為乙等，0.5 ～ 0.7 為丙等［22］。

3.2 模型率定的過(guò)程與結(jié)果

3.2.1 徑流過(guò)程的模擬。使用1980—2005年?yáng)|江流域博羅水文站月平均天然徑流量資料進(jìn)行徑流的率定與校驗(yàn)，1980—1983年數(shù)據(jù)用于模型預(yù)熱。1980—1992年的徑流模擬過(guò)程采用同期的氣象數(shù)據(jù)、1980年代下墊面數(shù)據(jù)作為模型的輸入數(shù)據(jù);1993—2005年的徑流模擬過(guò)程采用同期的氣象數(shù)據(jù)、2000年代下墊面數(shù)據(jù)作為模型的輸入數(shù)據(jù)。月徑流過(guò)程在模型校正期和驗(yàn)證期都取得了較好的模擬效果(圖1、圖2)。在徑流校正和驗(yàn)證中年均徑流Dv分別為0.53%和5.00%，Dv小于10%;ENS 分別為0.926 和0.894，大于 0.7。

圖1 1984—1992年徑流模擬校準(zhǔn)結(jié)果Fig.1 Calibration results of flow simulation during 1984—1992

圖2 1993—2005年徑流模擬驗(yàn)證結(jié)果Fig.2 Verification results of flow simulation during 1993—2005

3.2.2 泥沙過(guò)程的模擬。泥沙過(guò)程的模擬是進(jìn)行非點(diǎn)源模擬的一個(gè)重要環(huán)節(jié)，泥沙模擬在校準(zhǔn)好后的水量模擬的基礎(chǔ)上進(jìn)行。模擬結(jié)果見(jiàn)圖3、圖4。泥沙模擬結(jié)果誤差在可以接受的范圍內(nèi)，率定期和驗(yàn)證期的Dv分別為12.87%和5.95%，ENS 分別為 0.616 和 0.664。泥沙模擬結(jié)果比水文過(guò)程模擬結(jié)果誤差大主要由2個(gè)原因造成:①用于水文過(guò)程率定的徑流量是對(duì)實(shí)測(cè)徑流進(jìn)行還原后的徑流量，但實(shí)測(cè)輸沙量是在實(shí)測(cè)徑流情況下的輸沙量，東江的徑流與河道泥沙輸沙過(guò)程受到水利工程的影響，這是誤差較大的原因之一;② 精確的泥沙輸移校準(zhǔn)需要懸移質(zhì)泥沙成分及顆粒大小等比較詳細(xì)的泥沙觀測(cè)斷面的資料，本研究是對(duì)東江流域長(zhǎng)期的水文和產(chǎn)匯污過(guò)程進(jìn)行模擬，用于模型建立和校準(zhǔn)的實(shí)測(cè)值為月均觀測(cè)值，所采用的參數(shù)多為經(jīng)驗(yàn)參數(shù)，因此，未能對(duì)模型進(jìn)行更加精準(zhǔn)和更小時(shí)間尺度的校準(zhǔn)，但模擬結(jié)果仍然可以代表泥沙的月變化過(guò)程，可以用于對(duì)于東江流域水文和泥沙輸送的長(zhǎng)期變化過(guò)程的研究。

圖3 1984—1992年泥沙模擬校準(zhǔn)結(jié)果Fig.3 Calibration results of sediment simulation during 1984—1992

3.2.3 污染物的模擬。河流水環(huán)境中的污染物來(lái)源包括點(diǎn)源與非點(diǎn)源，東江流域面積較大，尤其是中下游區(qū)域工業(yè)與城鎮(zhèn)生活污水集中排放等點(diǎn)源污染分布較密集，獲取與整理統(tǒng)計(jì)各子流域多年點(diǎn)源污染的排放量與入河量的數(shù)據(jù)不僅工作量很大而且較準(zhǔn)確的數(shù)據(jù)的獲取也非常困難，如果獲得的數(shù)據(jù)不準(zhǔn)確反而會(huì)對(duì)模型的建立過(guò)程產(chǎn)生負(fù)面的影響。因此，本研究在參考其他相關(guān)研究［23-26］的基礎(chǔ)上，采用了數(shù)字濾波法對(duì)東江流域的實(shí)測(cè)污染負(fù)荷系列數(shù)據(jù)進(jìn)行分割，得到非點(diǎn)源負(fù)荷數(shù)據(jù)，作為HSPF模型模擬結(jié)果的比對(duì)，該方法與實(shí)測(cè)數(shù)據(jù)的非點(diǎn)源污染負(fù)荷分割結(jié)果已經(jīng)應(yīng)用在東江流域的徑流與非點(diǎn)源污染研究［25-26］中，采用該方法能夠減小中尺度流域水文模型建模中由點(diǎn)源污染物排放量及入河量誤差引起的模擬誤差。

模擬結(jié)果(圖5、圖6)顯示，年均氨氮、總磷負(fù)荷模擬結(jié)果的 Dv分別為 9.43%和13.76%，ENS分別為0.653和0.783。對(duì)于東江流域這樣一個(gè)中尺度的流域，利用構(gòu)建的HSPF模型模擬得到的氨氮與總磷模擬結(jié)果誤差在可以接受的范圍內(nèi)，模擬誤差主要來(lái)自徑流、泥沙過(guò)程模擬的誤差累積以及實(shí)測(cè)數(shù)據(jù)本身的誤差。

圖5 2001—2005年?yáng)|江流域氨氮非點(diǎn)源污染負(fù)荷模擬結(jié)果Fig.5 Results of NH3-N load simulation in the Dongjiang River during 2001—2005

圖6 2001—2005年?yáng)|江流域總磷非點(diǎn)源污染負(fù)荷模擬結(jié)果Fig.6 Results of TP load simulation in the Dongjiang River during 2001—2005

4 非點(diǎn)源污染的時(shí)空分布特征

4.1 時(shí)間分布

1999—2005年，氨氮、總磷非點(diǎn)源負(fù)荷差異較明顯，呈現(xiàn)先升后降再升的變化趨勢(shì)。2000與2005年(豐水年)氨氮、總磷非點(diǎn)源負(fù)荷是2004年(枯水年)負(fù)荷的3倍以上，豐水年是東江流域非點(diǎn)源污染物流失的關(guān)鍵年份。

非點(diǎn)源污染的發(fā)生和輸送與年內(nèi)降雨的分布及徑流量有直接的關(guān)系。降雨量增加，徑流產(chǎn)生量多，則相應(yīng)地水土流失程度加大，氮、磷流失量增加。圖7、圖8分別表示1999—2005年平均月降雨量、氨氮和總磷非點(diǎn)源負(fù)荷的年內(nèi)分布。東江流域徑流年內(nèi)變化與降雨基本一致，1—3月變化平緩，4月大幅增大，峰值出現(xiàn)在6月份，4—9月徑流合計(jì)152億 m3，占年總徑流量的75.06%。氨氮和總磷的峰值變化明顯，并與降雨年內(nèi)分布基本一致，說(shuō)明東江流域非點(diǎn)源污染的形成分布受降雨尤其是暴雨的影響較大;氨氮非點(diǎn)源負(fù)荷峰值出現(xiàn)在4—6月，4—6月的氨氮非點(diǎn)源負(fù)荷總和占年總負(fù)荷量的57.8%;總磷非點(diǎn)源負(fù)荷峰值出現(xiàn)在4—9月，4—9月的氨氮非點(diǎn)源負(fù)荷總和占年總負(fù)荷量的91.63%。

圖7 1999—2005年氨氮月負(fù)荷模擬結(jié)果Fig.7 Average monthly NH3-N load of years(1999—2005)simulated by HSPF

圖8 1999—2005年總磷月負(fù)荷模擬結(jié)果Fig.8 Average monthly TP load of years(1999—2005)simulated by HSPF

東江流域非點(diǎn)源污染主要發(fā)生在降雨徑流期間，汛期的污染負(fù)荷占有較大比例，因此，對(duì)5—8月的徑流量、氨氮、總磷非點(diǎn)源污染負(fù)荷再進(jìn)行統(tǒng)計(jì)［27］，1999—2005年氨氮、總磷5—8月的非點(diǎn)源負(fù)荷占全年總量比例大部分均在50%以上;豐水年5—8月汛期非點(diǎn)源負(fù)荷比例高于枯水年;降雨尤其是暴雨對(duì)非點(diǎn)源污染的產(chǎn)污匯污影響很大，汛期是東江流域土壤侵蝕和氮、磷非點(diǎn)源污染防治的關(guān)鍵時(shí)期(圖9)。

4.2 空間分布

非點(diǎn)源污染負(fù)荷的產(chǎn)生輸出受流域內(nèi)降雨分布、土地利用方式等因素的綜合作用。以2001年(豐水年)為例，東江流域降雨量的空間分布總體上隨地形的變化由北向南遞增，在新豐江水庫(kù)區(qū)域以及東江下游博羅附近年降雨量最大，2001年氨氮、總磷的負(fù)荷分布總體上也呈由北向南遞增的態(tài)勢(shì)(圖10)，同時(shí)，非點(diǎn)源污染的產(chǎn)生與分布不完全取決于降雨分布，還與土地利用、植被覆蓋、地形分布等因素有關(guān)，非點(diǎn)源單位面積負(fù)荷值在東江流域西北部較小，在東江中下游分布較大，對(duì)比土地利用類型圖，在東江中下游源城區(qū)、博羅縣、惠城區(qū)、惠陽(yáng)區(qū)、龍崗區(qū)中農(nóng)業(yè)用地與建設(shè)用地集中的地塊氨氮、總磷負(fù)荷密度最大。東江流域非點(diǎn)源污染的關(guān)鍵區(qū)域主要位于流域中下游源城區(qū)、博羅縣、惠城區(qū)耕地與建設(shè)用地集中的區(qū)域。

圖10 非點(diǎn)源污染負(fù)荷的空間分布Fig.10 Spatial distribution of non-pointed NH3-N and TP load

5 降雨量與非點(diǎn)源污染的統(tǒng)計(jì)分析

東江流域非點(diǎn)源污染的時(shí)間分布受流域內(nèi)降雨的影響較大，為進(jìn)一步分析東江流域降雨與非點(diǎn)源污染的時(shí)間分布特征的關(guān)系，對(duì)東江流域2001—2010年氨氮、總磷非點(diǎn)源污染物月負(fù)荷模擬值(采用2000年下墊面數(shù)據(jù)、2001—2010年日降雨數(shù)據(jù)進(jìn)行模擬)與月降雨量進(jìn)行回歸分析。

運(yùn)用SPSS 16.0統(tǒng)計(jì)軟件進(jìn)行回歸分析。通過(guò)對(duì)回歸分析曲線的綜合比選，復(fù)合曲線可以較好地?cái)M合2001—2010年氨氮與總磷月負(fù)荷數(shù)據(jù)序列與月降雨量數(shù)據(jù)序列，擬合的復(fù)合曲線回歸方程為。

式中:Y表示氨氮月負(fù)荷(t);Z表示總磷月負(fù)荷(t);X表示月降雨量(mm)。

氨氮負(fù)荷與降雨量擬合優(yōu)度R2等于0.831，回歸方程和回歸系數(shù)通過(guò)顯著性檢驗(yàn)，回歸方程顯著性檢驗(yàn)F值等于582.085，概率P值小于0.01;總磷負(fù)荷與降雨量擬合優(yōu)度R2等于0.719，回歸方程和回歸系數(shù)都通過(guò)顯著性檢驗(yàn)，回歸方程顯著性檢驗(yàn)F值等于302.552，概率P值小于 0.01。

6 土地利用與非點(diǎn)源污染統(tǒng)計(jì)分析

將東江流域土地利用類型劃分為林地、耕地、園地、建設(shè)用地、草地、水域與未利用地7種類型，以子流域?yàn)檠芯繂卧?，統(tǒng)計(jì)各單元中各土地利用類型的面積比例與各單元非點(diǎn)源污染濃度的相關(guān)關(guān)系。結(jié)果表明各子流域單元氨氮非點(diǎn)源污染濃度與耕地、園地、林地、建設(shè)用地面積比例存在顯著的線性相關(guān)關(guān)系(概率P值小于0.01);其中，耕地、園地、建設(shè)用地與非點(diǎn)源污染氨氮濃度呈正相關(guān)(相關(guān)系數(shù)分別是 0.608，0.710，0.718)，林地面積與氨氮濃度呈負(fù)相關(guān)(相關(guān)系數(shù)為-0.719)。

各子流域總磷非點(diǎn)源污染濃度也與耕地、園地、林地、建設(shè)用地面積比例存在顯著的線性相關(guān)關(guān)系(概率P值小于0.01);其中，耕地、園地、建設(shè)用地與非點(diǎn)源污染總磷濃度呈正相關(guān)(相關(guān)系數(shù)分別是0.481，0.659，0.714)，林地面積與總磷濃度呈負(fù)相關(guān)(相關(guān)系數(shù)為-0.649)。

非點(diǎn)源污染受到各類土地利用類型所占的面積比例的影響，土地利用結(jié)構(gòu)組成和空間分布特征不同，土地利用與水質(zhì)的關(guān)系亦不同［28］。進(jìn)一步以匯水分區(qū)為研究單元，采用相關(guān)分析、逐步回歸與聚類分析等方法，分析東江流域土地利用空間格局與非點(diǎn)源污染的關(guān)系。結(jié)果表明:(1)污染物的空間分布與用地類型關(guān)系復(fù)雜，污染物濃度的分布不僅與各類土地利用類型面積比例有關(guān)，還與各類用地類型在空間上的分布有關(guān)。(2)在將各匯水分區(qū)根據(jù)土地利用空間格局進(jìn)行聚類分析后，再采用相關(guān)分析、逐步回歸方法對(duì)土地利用類型分布與非點(diǎn)源污染的關(guān)系進(jìn)行統(tǒng)計(jì)分析，結(jié)果會(huì)受到樣本數(shù)量的較大影響。(3)總體上，在上中游林地面積較大的匯水單元，林地與污染物濃度的相關(guān)系數(shù)最大，林地對(duì)于緩解由其他用地類型產(chǎn)生的污染起了積極的作用;建設(shè)用地對(duì)于氨氮、總磷水質(zhì)濃度的貢獻(xiàn)隨用地面積比例的增大有逐漸增長(zhǎng)的趨勢(shì);耕地面積比例與污染物濃度的相關(guān)關(guān)系在不同類型的匯水單元中表現(xiàn)出有較大差異，這與其他相關(guān)研究成果的統(tǒng)計(jì)結(jié)果一致［28-29］。

7 結(jié)論

(1)應(yīng)用HSPF水文模型建立了東江流域徑流與非點(diǎn)源污染模擬模型。多年徑流過(guò)程模擬的模型率定期與驗(yàn)證期Dv分別為0.53%和5.00%，ENS分別為0.926和0.894，泥沙過(guò)程模擬的率定期和驗(yàn)證期Dv分別為12.87%和 5.95%，ENS 分別為 0.616 和0.664;氨氮、總磷負(fù)荷模擬的Dv分別為9.43%和13.76%，ENS分別為0.653和0.783。徑流、泥沙與污染物模擬過(guò)程與實(shí)測(cè)過(guò)程的總體趨勢(shì)基本一致，所建立的基于HSPF的模型可以滿足東江流域徑流與非點(diǎn)源污染長(zhǎng)期變化過(guò)程與時(shí)空分布特征研究的需要。

(2)模擬結(jié)果表明東江流域不同降雨水平年的氨氮、總磷非點(diǎn)源污染負(fù)荷差異明顯，豐水年是東江流域非點(diǎn)源污染物流失的關(guān)鍵年份。東江流域徑流年內(nèi)變化與降雨基本一致;非點(diǎn)源污染的發(fā)生與輸送和年內(nèi)降雨的分布及徑流量有直接的關(guān)系，東江流域非點(diǎn)源污染負(fù)荷的年內(nèi)分布特征受流域內(nèi)降雨的年內(nèi)分布影響顯著，汛期氨氮、總磷的非點(diǎn)源污染負(fù)荷占全年負(fù)荷總量比例基本在50%以上;豐水年汛期非點(diǎn)源污染負(fù)荷比重高于枯水年。降雨尤其是暴雨對(duì)東江流域非點(diǎn)源污染的產(chǎn)污匯污影響很大，汛期是東江流域土壤侵蝕和氮、磷非點(diǎn)源污染防治的關(guān)鍵時(shí)期。

(3)非點(diǎn)源單位面積負(fù)荷值在東江流域西北部較小，在東江中下游分布較大，對(duì)比土地利用類型圖，在東江中下游源城區(qū)、博羅縣、惠城區(qū)、惠陽(yáng)區(qū)、龍崗區(qū)中農(nóng)業(yè)用地與建設(shè)用地集中的地塊氨氮、總磷負(fù)荷密度最大。東江流域非點(diǎn)源污染的關(guān)鍵區(qū)域主要位于流域中下游源城區(qū)、博羅縣、惠城區(qū)耕地與建設(shè)用地集中的區(qū)域。

(4)東江流域非點(diǎn)源污染的時(shí)間分布受流域內(nèi)降雨的影響較大。為進(jìn)一步分析東江流域降雨與非點(diǎn)源污染的時(shí)間分布特征的關(guān)系，對(duì)東江流域2001—2010年氨氮、總磷污染物月負(fù)荷模擬值與月降雨量進(jìn)行了回歸分析。經(jīng)過(guò)綜合比選，結(jié)果表明:復(fù)合回歸曲線可較好地?cái)M合污染物月負(fù)荷的變化，氨氮、總磷擬合優(yōu)度R2分別為0.831和0.719，回歸方程和回歸系數(shù)顯著性檢驗(yàn)均通過(guò)檢驗(yàn)，概率P值小于0.01。

(5)以子流域?yàn)檠芯繂卧?，?duì)各單元中各類型土地利用的面積比例與非點(diǎn)源污染濃度的相關(guān)關(guān)系進(jìn)行了統(tǒng)計(jì)分析。結(jié)果表明各子流域氨氮、總磷濃度與耕地、園地、林地、建設(shè)用地面積比例存在顯著的線性相關(guān)關(guān)系，其中，耕地、園地、建設(shè)用地與非點(diǎn)源污染濃度呈正相關(guān)，林地面積與非點(diǎn)源污染濃度呈負(fù)相關(guān)。

［1］ 付永鋒，陳文輝，趙基花.非點(diǎn)源污染的研究進(jìn)展與前景展望［J］.山西水利科技，2003(3):32-35.

［2］ 張殿發(fā)，王世杰，李瑞玲.土地利用/土地覆被變化對(duì)長(zhǎng)江流域水環(huán)境的影響研究［J］.地域研究與開(kāi)發(fā)，2003，22(1):69-72.

［3］ 李俊然，陳利頂.土地利用結(jié)構(gòu)對(duì)非點(diǎn)源污染的影響［J］.中國(guó)環(huán)境科學(xué)，2000，20(6):506-510.

［4］ 于興修，楊桂山，梁濤.西苕溪流域土地利用對(duì)氮素徑流流失過(guò)程的影響［J］.農(nóng)業(yè)環(huán)境保護(hù)，2002，21(5):424-427.

［5］ 黃云鳳，張珞平，洪華生.不同土地利用對(duì)流域土壤侵蝕和氮、磷流失的影響［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2004，23(4):735-739.

［6］ 張曉明，曹文洪，余新曉，等.黃土丘陵溝壑區(qū)典型流域徑流輸沙對(duì)土地利用/覆被變化的響應(yīng)［J］.應(yīng)用生態(tài)學(xué)報(bào)，2009，20(1):121-127.

［7］ 唐麗霞，張志強(qiáng)，王新杰，等.黃土高原清水河流域土地利用/覆蓋和降雨變化對(duì)侵蝕產(chǎn)沙的影響［J］.自然資源學(xué)報(bào)，2010，25(8):1340-1349.

［8］ 章程，袁道先.典型巖溶地下河流域水質(zhì)變化與土地利用的關(guān)系［J］.水土保持學(xué)報(bào)，2004，18(5):134-137.

［9］ 黃金良，李青生，洪華生，等.九龍江流域土地利用/景觀格局-水質(zhì)的初步關(guān)聯(lián)分析［J］.環(huán)境科學(xué)，2011，32(1):64-72.

［10］ 陳瑩，許有鵬，尹義星.土地利用/覆被變化下的暴雨徑流過(guò)程模擬分析——以太湖上游西苕溪流域?yàn)槔跩］.地理科學(xué)，2009，29(1):117-123.

［11］ 王秀娟，劉瑞民，宮永偉，等.香溪河流域土地利用格局演變對(duì)非點(diǎn)源污染的影響研究［J］.環(huán)境工程學(xué)報(bào)，2011，5(5):1194-1200.

［12］ 陳媛，郭秀銳，程水源，等.基于SWAT模型的三峽庫(kù)區(qū)大流域不同土地利用情景對(duì)非點(diǎn)源污染的影響研究［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2012，31(4):798-806.

［13］ You Y Y，Jin W B，Xiong Q X，et al.Simulation and Validation of Non-point Source Nitrogen and Phosphorus Loads under Different Land Uses in Sihu Basin，Hubei Province，China［J］.Procedia Environmental Sciences，2012，13:1781-1797.

［14］ 徐宗學(xué).水文模型［M］.北京:科學(xué)出版社，2009.

［15］ 郝振純，李麗，王加虎，等.分布式水文模型理論與方法［M］.北京:科學(xué)出版社，2010.

［16］ Arnaud P，Bouvier C，Cisneros L，et al.Influence of Rainfall Spatial Variability on Flood Prediction［J］.Journal of Hydrology，2002，260(1-4):216-230.

［17］ Lee S，Ni-mesister W，Toll D.Assessing the Hydrologic Performance of the EPA’s Nonpoint Source Water Quality Assessment Decision Support Tool Using North American Land Data Assimilation System(NLDAS)Products［J］.Journal of Hydrology，2010，387(3):212-220.

［18］ Nash J E，Sutcliffe J V.River Flow Forecasting through Conceptual Models—A Discussion of Principles［J］.Journal of Hydrology，1970，10(3):282-290.

［19］ ASCE Task Committee.The ASCE Task Committee on Definition of Criteria for Evaluation of Watershed Models of the Watershed Management Committee Lrrigation and Drainage Division.Criteria for Evaluation of Watershed Models［J］.Journal of Lrrigation and Drainage Engineering，1993，119(3):429-442.

［20］ Donigian A S，Imhoff J C，Bicknell B R，et al.EPA-600/3-84-065 Application Guide for Hydrological Simulation Program:FORTRAN(HSPF)［S］.Athens GA:U.S.Environmental Protection Agency，1984.

［21］ 林誠(chéng)二，村上正吾，渡邊正孝，等.基于全球降水?dāng)?shù)據(jù)估計(jì)值的地表徑流模擬——以長(zhǎng)江上游地區(qū)為例［J］.地理學(xué)報(bào)，2004，59(1):125-135.

［22］ 董延軍，李杰，鄭江麗，等.流域水文水質(zhì)模擬軟件(HSPF)應(yīng)用指南［M］.鄭州:黃河水利出版社，2009.

［23］ 林凱榮，陳曉宏，江濤，等.數(shù)字濾波進(jìn)行基流分割的應(yīng)用研究［J］.水力發(fā)電，2008(6):28-30.

［24］ 李懷恩.估算非點(diǎn)源污染負(fù)荷的平均濃度法及其應(yīng)用［J］.環(huán)境科學(xué)學(xué)報(bào)，2000，20(4):397-400.

［25］ 黎坤.華南地區(qū)坡面溶解氮流失機(jī)理及流域非點(diǎn)源負(fù)荷解析方法研究［D］.廣州:中山大學(xué)，2009.

［26］ 黎坤，林凱榮，江濤，等.數(shù)字濾波法在點(diǎn)源和非點(diǎn)源污染負(fù)荷分割中的應(yīng)用［J］.環(huán)境科學(xué)研究，2010，23(3):298-303.

［27］ 李明濤，王曉燕，段淑懷，等.潮河流域非點(diǎn)源污染負(fù)荷的時(shí)空分布模擬研究［C］//武漢大學(xué).Proceedings of Conference on Environmental Pollution and Public Health.武漢:Scientific Research Publishing，2012.

［28］ 郭青海，馬克明，張易.城市土地利用異質(zhì)性對(duì)湖泊水質(zhì)的影響［J］.生態(tài)學(xué)報(bào)，2009，29(2):776-787.

［29］ 劉松波.滇池流域土地利用變化與入湖河流水質(zhì)關(guān)系研究［J］.環(huán)境科學(xué)導(dǎo)刊，2013，32(2):42-44.