氣候變化和人類活動對水庫面源污染負荷影響

楊金尉

(遼寧省鐵嶺水文局，遼寧 鐵嶺 112000)

當前，氣候變化和人類活動對水循環的影響逐年增強。極端天氣頻發、多發，水體作為面源污染運移的主體，水循環的變化勢必然造成水體污染指標的變化?？茖W評估氣候變化和人類活動共同影響下的水庫面源污染對于水庫水環境保護十分重要。當前，單一研究氣候變化或者人類活動對水庫面源污染影響已取得一定的研究成果，但是結合氣候變化和人類活動共同影響下的水庫面源污染影響的研究還較少，而一般來說，氣候變化和人類活動是伴隨式發生的，因此單一考慮氣候變化或者人類活動，很難全面評估變化環境下的水庫面源污染負荷影響，因此需要結合人類活動和氣候變化兩種影響要素一起，分析其變化對水庫面源污染的影響。為此，本文以當前分析氣候變化和人類活動較為成熟的模擬平臺SWAT模型為依托，通過設定人類活動和氣候變化情景，分析其共同影響下的水庫面源污染響應。

1 研究方法和數據

1.1 研究方法

本文在SWAT模型適用性分析基礎上，采用SWAT模型分析氣候變化和人類活動對水庫庫區面源污染負荷的影響。結合人類活動及氣候變化貢獻率定量分離方法分析其對水庫面源污染負荷影響的具體貢獻率。

1.2 采用數據

本文主要以遼寧北部某水庫為研究實例。模型所需要輸入的為庫區所在流域的下墊面數據，包括數字高程、土地利用以及土壤類型數據，如圖1所示。從圖1中可看出區域主要利用方式為耕地和林地。此外，由于需要分析水庫面源污染負荷，對庫區主要土壤類型進行化學試驗分析，分析結果見表1。

2 模型適用性分析

首選需要分析SWAT模型在區域的適用性?；趨^域2000—2010年實測水質數據，對SWAT模型非點源污染模擬進行參數率定和驗證，參數率定結果見表2，模型驗證結果見表3—4及圖2—3。

從表3—4中可看出，SWAT模型在區域面源污染負荷方面具有較好的適用性，總氮和總磷預測誤差均較好。模型在污染負荷量上具有較好的精度外，對其污染負荷過程也具有較好的吻合度，從確定性系數指標中可看出，面源污染負荷模擬的確定性系數值在0.5以上，達到乙級面源污染負荷預測的要求。從選取的典型年面源污染負荷模擬結果也可看出，在各月份面源污染負荷中，模型模擬的面源污染負荷年內分配過程和水質監測數據也具有較好的吻合度。水庫面源污染負荷的高峰期主要出現在7—9月份，主要因為這三個月份屬于豐水期，豐水期水量的增多，使得進入庫區的面源污染負荷也增多，而SWAT模型對面源污染負模擬主要依托于水量為載體，因此在豐水期模型的吻合度明顯高于枯水期的吻合度。

圖1 模型輸入下墊面數據

單位：mg/kg

表2 模型參數率定結果

表3 SWAT模型總氮驗證結果

表4 SWAT模型總磷驗證結果

圖2 代表年份總氮模擬成果圖

圖3 代表年份總磷模擬成果圖

3 空間反演結果

結合SWAT模型模擬結果，對庫區面源污染負荷進行空間分析，空間分析成果如圖4所示。

圖4 面源污染負荷空間模擬成果

從圖4中可看出，面源污染負荷大的區域主要集中在區域北部和南部，北部區域面源污染負荷大的原因在于區域主要為農業耕地，耕地措施下農業化肥所產生面源污染負荷在降水較多的月份會進入庫區，影響庫區的面源污染。而南部面源污染主要增多的原因在于其主要為居民用地，人類活動下的點源排放量是這一區域面源污染負荷較為集中的原因。此外，區域面源污染負荷也主要受到區域降水量分布的影響，降水量較多的區域，其面源污染負荷一般較多，而降水量較少的區域，其面源污染負荷一般較少。

4 氣候變化及人類活動共同影響下的面源污染負荷分析

4.1 情景分析

情景模式設定見表5。結合該模型模擬不同變化情景模式下的面源污染負荷，并對比變化前的面源污染負荷，對比結果見表6—7。

本文將氣候變化和人類活動影響統稱為變化情景模式，在設定變化情景模式時，主要依據全球氣候變化模式以及區域土地利用遠景規劃數據。變化情景模式I、II，氣溫升高、土地利用變化主要為耕地和林地之間的轉變，從兩個情景模式分析可以看出，情景I氣溫升高且區域退耕還林，庫區總體面源污染負荷下降趨勢明顯，削減率5.2%～9.1%；情景II氣溫降低但退耕還林措施力度加大，庫區總體面源污染進一步削減，削減率明顯高于情景I，削減率為18.9%～27.6%之間。而在情景III下，降雨增加，勢必造成入庫面源污染負荷的增加，而城鎮用地增加，使得點源污染排放量增加，進一步加劇水庫面源污染的負荷。在情景IV下，降雨量的減少對水庫面源污染負荷影響超過城鎮用地面積的變化。

表5 變化情景模式設定

注：氣候模式設定依據為全球氣候變化模式，土地利用設定依據為區域土地利用遠景規劃數據。

表6 變化情景方案下的總氮響應 單位：104t

表7 氣溫變化情景方案下的總磷響應單位：104t

變化情景模式變化前模擬總磷變化后模擬總磷總磷失變化總磷變化率/%I55.149.9-5.2-8.5II55.139.9-15.2-27.6III55.164.49.315.1IV55.147.8-7.3-11.9

4.2 貢獻率分析

結合貢獻率分析方法對氣候變化及人類活動對水庫面源污染進行定量分離，結果見表8—9。

從貢獻率分離結果可看出，氣候變化和人類活動對總氮和總磷影響比例較為相似。從各個年份分離結果可看出，氣候變化對總氮影響的貢獻率均值為36%，這主要是因為氣候變化主要體現在對水循環的影響，通過影響水量的變化來影響進入庫區的面源污染負荷，而人類活動除表現為對水量的影響外，還影響產生面源污染負荷的途徑，退耕還林措施后，區域農業面源污染負荷得到控制，因此進入庫區的面源污染負荷得到減少。此外，城鎮用地面積增加，勢必造成進入庫區的點源排放量增加，也勢必增加其污染負荷，因此人類活動對面源污染影響的途徑較多，使得其對面源污染負荷影響的貢獻率均值達到64%。

表8 水庫庫區總氮影響因素分解 單位：104t

表9 水庫庫區總磷影響因素分解 單位：104t

5 結論

(1)降水及城鎮用地變化是水庫庫區面源污染負荷增多的主因，氣溫變化對面源污染負荷影響程度偏低。

(2)退耕還林措施是有效降低庫區面源污染負荷的一種有效措施，但是也會減少進入庫區的水量，水庫面源污染負荷控制還應重點注重點源排放量。

(3)本文沒有分析工業點源排放量變化對庫區面源污染負荷影響，在以后的研究中還應重點加入工業點源排放數據進行分析。