大伙房水庫調水對上游流域土壤侵蝕變化的響應

王之中

(撫順市水利技術推廣中心，遼寧 撫順 113008)

隨著城市的快速發展，人類活動導致土壤侵蝕加劇，我國水環境污染問題日益突出[1]。水質安全關系到居民飲用水安全，保證水源地水質安全對社會的發展具有重要意義。土地利用對水環境影響的相關研究不斷得到重視，Hwang[2]認為土壤侵蝕程度對地表水的質量起著決定性作用。巨拓[3]研究認為水庫水質污染主要來自土壤侵蝕造成的非點源污染。陶亮[4]對南水北調工程水源區進行了生態泥沙和環境泥沙及泥沙資源化研究，認為土壤侵蝕與水環境的惡化有著密切聯系。國內外學者對水質評價與影響因素之間的研究，取得了一些成果，但結合上游土地利用格局變化對調水輸入水質影響的研究少有涉獵。本文以大伙房水庫上游流域土壤侵蝕數據為基礎，將土壤侵蝕變化過程與水庫水質進行關聯分析，明確了土壤侵蝕程度與水質間關聯程度，為更有效保護水庫水質提供理論指導，有著重要的現實意義。

1 研究方法

1.1 水質數據

大伙房水庫固定兩個水質觀測斷面，包括入庫河流入水口觀測斷面(北雜木斷面)和水庫渾河出水口處觀測斷面(撫順取水口斷面)，數據觀測以自動觀測儀器與人工取樣測定分析相結合的方式。水質觀測指標包括大腸菌群(Coli-group)、生化需氧量(BOD)、高錳酸鉀指數(CODmn)、溶解氧(DO)、總磷(TP)、總氮(TN)、氨氮(NH3-N)，調水后增加了Chl-a。

1.2 土壤侵蝕與土地利用數據

根據GB/T 21010—2007《土地利用現狀分類》，對土地利用類型進行劃分，土地利用數據以水土保護部門多年統計數據和實地調查為準[5- 6]。為保證數據的可靠性，利用遙感普查的統計數據進行校核驗證[7- 8]。不同侵蝕類型占總侵蝕面積的比例如圖1所示。

圖1 不同侵蝕類型面積占總侵蝕面積的比例

調水前2008年與調水后2017年大伙房水庫上游土壤侵蝕程度變化見表1。可以看出，上游土壤侵蝕的面積有所減少，其中輕度和中度侵蝕面積減小最多，變化率分別為-25.74%和-9.37%。而強烈侵蝕、極強烈侵蝕面積有所增加，變化幅度分別為317.32%和34.32%，劇烈侵蝕面積也呈現出增加的趨勢，占現在總侵蝕比例的0.28%。

表1 大伙房水庫上游土壤侵蝕程度變化情況

2 調水前水質對上游土壤侵蝕的響應

2.1 北雜木斷面水質的響應

上游流域土壤侵蝕程度與北雜木斷面水質相關分析結果見表2。

可以看到，上游土壤侵蝕面積的百分比與水質的總氮、總磷以及大腸桿菌呈正相關關系，說明土壤的侵蝕導致土壤中的N、P等元素流失，造成水庫中水體和下游水系中含N、P污染物増加，導致大伙房水庫的水質惡化，下游水體富營養化；土壤侵蝕的面積比和BOD呈現顯著的負相關，說明土壤侵蝕面積的增加并不會對下游的BOD產生影響。上游土壤侵蝕主要是輕度侵蝕，其百分比與總氮、總磷和大腸桿菌是正相關的關系，說明下游氮和磷含量的增多主要受土壤輕度侵蝕的影響[9]。土壤輕度侵蝕面積和BOD的相關性呈顯著的負相關，說明土壤輕度侵蝕面積的增加并不會使下游水體中BOD含量發生變化。極強烈侵蝕、強烈侵蝕和中度侵蝕面積所占比例與總磷含量的相關性為負，說明這三種侵蝕類型的面積并沒有對下游水系輸送過多含氮磷的物質，對于下游水系輸送磷、氮的貢獻不大。極強烈、強烈和中度侵蝕程度還與大腸桿菌呈現一種負相關的關系，說明這三種侵蝕面積的增加并不會使下游大腸桿菌含量增加。

大伙房水庫上游劇烈侵蝕面積最小，其所占的比例與幾種指標并沒有明顯的相關性。輕度侵蝕分布區域是造成氮磷元素流失的主要部分，是下游氮磷元素的主要來源，表層土壤的流失是造成下游水質惡化的主要因素。

2.2 撫順取水口斷面水質的響應

上游流域土壤侵蝕程度與撫順取水口斷面水質相關分析結果見表3。

可以看出，大伙房水庫上游土壤侵蝕面積的百分比和水中總氮、總磷以及大腸桿菌的含量呈現正相關關系，說明土壤侵蝕導致N、P等元素流失，造成大伙房水庫及下游河水中N、P元素的增加，使得下游水體富營養化，大腸桿菌大量繁殖，水庫水質惡化。土壤侵蝕的面積比與BOD呈顯著負相關，說明土壤侵蝕面積的增加并不會使下游水體中BOD的含量發生變化。

表2 上游流域土壤侵蝕程度與水質相關分析

表3 上游流域土壤侵蝕程度與水質相關分析

上游流域土壤輕度侵蝕所占的比例與總磷以及大腸桿菌數量呈顯著正相關，說明輕度侵蝕是造成氮磷等營養元素流失到大伙房水庫中的主要因素。土壤輕度侵蝕和BOD的相關性為顯著負相關，說明土壤輕度侵蝕面積的增加并不使下游的BOD含量發生變化。

極強烈侵蝕、強烈侵蝕和中度侵蝕的面積比與總磷、大腸桿菌含量的相關性為負，說明這三種侵蝕類型對下游水系輸送磷、氮的貢獻不大。極強烈、強烈和中度侵蝕程度的面積百分比與BOD呈正相關的關系，說明這三種侵蝕面積的增加會使撫順斷面水中BOD含量上升。劇烈侵蝕與總氮、高錳酸鹽含量指數以及BOD呈正相關關系，說明劇烈侵蝕對下游水域氮元素的貢獻也較大，水體中BOD和高錳酸鹽含量增加。

3 調水后水質對上游流域土壤侵蝕的響應

3.1 北雜木斷面水質的響應

上游流域土壤侵蝕程度與北雜木斷面水質相關分析結果見表4。

可以看出，調水后上游土壤侵蝕面積比與水中的總氮、總磷以及葉綠素a含量具有顯著的正相關，說明隨著土壤侵蝕面積的增加，土壤中的氮磷元素會逐步流失，使得下游大伙房水庫中含氮磷元素超標，水體富營養化，促進了藻類植物的繁殖，同時也造成了下游水體中葉綠素a含量的增加，使大伙房水庫水質進一步下降。

輕度侵蝕所占比例與總氮、總磷及葉綠素a的相關性為正，輕度侵蝕造成大伙房水庫上游表層土壤的流失，土壤中氮磷元素被輸送到大伙房水庫下游，水域水體富營養化，促進了藻類繁殖，水體中葉綠素a的增加，使下游水體水質惡化加劇。

極強烈侵蝕、強烈侵蝕和中度侵蝕的面積比例與總磷、總氮含量的相關性為負，說明這三種侵蝕類型對下游水系輸送磷、氮的貢獻不大，一定程度上抑制了藻類的繁殖，水體中葉綠素a含量降低。也表明表層土壤的流失是造成氮磷等營養元素流失的主要原因，隨著侵蝕慢慢發展到強烈侵蝕，土壤中氮磷等營養元素流失減少。極強烈侵蝕、強烈侵蝕和中度侵蝕所占面積與大腸桿菌呈負相關關系，說明這三種侵蝕面積的增加不會造成北雜木斷面水質中大腸桿菌含量增加。

劇烈侵蝕面積與水庫水流中的氨氮含量呈正相關性，說明調水后劇烈侵蝕所占面積比的增加引起了下游水體中氨氮含量的增加，劇烈侵蝕的面積百分比也與高錳酸鉀指數呈正相關，說明調水后劇烈侵蝕所占面積比的增加引起了下游水體中高錳酸鉀指數的增加。

3.2 撫順取水口斷面水質的響應

上游流域土壤侵蝕程度與撫順取水口斷面水質相關分析結果見表5。

可以看出，調水后大伙房水庫上游土壤侵蝕面積比與水質中的總氮、總磷以及葉綠素a呈顯著的正相關，說明隨著土壤侵蝕面積的增加，土壤中的氮磷元素會逐漸流失，導致下游大伙房水庫中氮磷元素超標，水體富營養化嚴重，同時下游水體中葉綠素a含量的增加。土壤侵蝕面積比與大腸桿菌含量呈負相關，說明調水后大伙房水庫上游侵蝕面積的增加對撫順取水口斷面水中大腸桿菌的含量沒有直接影響。

輕度侵蝕所占面積的比例與總氮、總磷以及葉綠素a呈顯著的正相關，與高錳酸鹽指數以及大腸桿菌含量也是正相關，說明輕度侵蝕造成了大伙房水庫上游表層土壤的流失，使得土壤中的氮磷元素被輸送到了大伙房水庫下游，造成了下游水域水體的富營養化，促進了藻類的繁殖，進而使得水體中葉綠素a的增加，同時高錳酸鉀和大腸桿菌含量的增加進一步使得下游水體水質惡化。

表4 上游流域土壤侵蝕程度與水質相關分析

表5 上游流域土壤侵蝕程度與水質相關分析

極強烈侵蝕、強烈侵蝕和中度侵蝕的面積與古樓斷面上水質中的總磷、總氮以及葉綠素a具有顯著的負相關，與大腸桿菌也是負相關的關系，說明這三種侵蝕區域流失土壤中所含的氮、磷元素較少，使得下游流域和古樓水域中葉綠素a和大腸桿菌含量減少。中度侵蝕、強烈侵蝕和極強烈侵蝕所占的比例與高錳酸鉀指數的相關性為正，說明這三種侵蝕面積的增加，會引起下游撫順取水口斷面水中高錳酸鹽含量增加。

劇烈侵蝕與撫順取水口斷面水質指標之間并沒有顯著的相關性，因此劇烈侵蝕并未對該斷面水質有明顯的影響。

4 結論

以大伙房水庫上游流域土壤侵蝕數據為基礎，將土壤侵蝕變化過程與水庫水質進行關聯分析，明確了土壤侵蝕程度與水質間關聯程度，得出以下結論：

(1)調水前上游流域的土壤侵蝕面積與各水質指標之間具有明顯的相關性。土壤侵蝕面積百分比與氮磷含量呈顯著正相關，表明土壤侵蝕是造成下游水質惡化的重要因素。輕度侵蝕和營養素氮、磷之間具有顯著正相關，表明輕度侵蝕區域，氮磷元素的流失最多，對下游水域水體富營養化貢獻最大，而以深層土壤流失為主的其它幾種侵蝕方式對下游水體惡化影響不大。

(2)調水后上游流域的土壤侵蝕面積和各水質指標之間也有明顯的相關性。土壤侵蝕面積百分比與氮、磷和葉綠素a含量呈顯著正相關，表明土壤侵蝕是造成下游水質惡化的重要因素。輕度侵蝕和營養素氮、磷和葉綠素a之間具有顯著正相關，表明輕度侵蝕區域，氮磷元素的流失最多，對下游水域水體富營養化貢獻最大，其它幾種侵蝕方式影響不大。