贛江下游不同樣地降雨徑流氮磷污染負(fù)荷估算

向速林, 陶術(shù)平, 王逢武

1. 華東交通大學(xué)環(huán)境工程系, 南昌 3300132. 南昌大學(xué)鄱陽(yáng)湖環(huán)境與資源利用教育部重點(diǎn)實(shí)驗(yàn)室, 南昌 330047

贛江下游不同樣地降雨徑流氮磷污染負(fù)荷估算

向速林1,2,*, 陶術(shù)平1, 王逢武1

1. 華東交通大學(xué)環(huán)境工程系, 南昌 330013
2. 南昌大學(xué)鄱陽(yáng)湖環(huán)境與資源利用教育部重點(diǎn)實(shí)驗(yàn)室, 南昌 330047

向速林, 陶術(shù)平, 王逢武. 贛江下游不同樣地降雨徑流氮磷污染負(fù)荷估算[J]. 生態(tài)科學(xué), 2016, 35(2): 62-65.

XIANG Sulin, TAO Shuping, WANG Fengwu. Nitrogen and phosphorus pollution loads from the different land use types in the Ganjiang River downstream[J]. Ecological Science, 2016, 35(2): 62-65.

選擇贛江下游地區(qū)不同樣地類型的徑流小區(qū), 遇到大降雨時(shí)監(jiān)測(cè)徑流水樣中的氮、磷含量, 并利用SCS模型計(jì)算降雨徑流量, 從而估算不同樣地類型降雨徑流中氮、磷的污染負(fù)荷。結(jié)果表明, 不同土地利用類型中氮、磷的輸出量差異顯著, 其氮、磷的污染負(fù)荷均以水田最高, 分別達(dá)到 10732.91 t·a–1與 1071.65 t·a–1, 而以林地最小, 分別為657.53 t·a–1與166.70 t·a–1, 反映了降雨徑流下水田的溶解態(tài)氮、磷流失量最大。因此, 水田的氮、磷流失是贛江下游需要重點(diǎn)考慮的因素。

贛江下游; 降雨徑流; 污染負(fù)荷

1 前言

降雨徑流對(duì)地表的擊濺以及沖刷產(chǎn)生的水力侵蝕作用, 造成了結(jié)合在土壤顆粒和土壤團(tuán)聚體中的氮、磷等營(yíng)養(yǎng)鹽的流失并進(jìn)入地表水體, 從而引起地表水體的富營(yíng)養(yǎng)化[1], 是一種重要的非點(diǎn)源污染類型[2]。流域土地利用方式及其結(jié)構(gòu)變化對(duì)徑流中營(yíng)養(yǎng)鹽含量以及輸出過程均具有極其重要的影響[3–4],不同土地利用類型決定了土壤的植被覆蓋度、物理結(jié)構(gòu)性質(zhì)以及氮、磷等營(yíng)養(yǎng)鹽含量[5]。不同來(lái)源方式的降雨徑流攜帶的污染物質(zhì)對(duì)地表水體氮、磷污染負(fù)荷具有重要的貢獻(xiàn), 可能成為地表水體富營(yíng)養(yǎng)化主要來(lái)源, 但不同土地利用類型氮、磷隨降雨徑流的輸出具有較大差異[6–7]。因此, 探討不同樣地類型下降雨徑流中氮、磷的污染負(fù)荷顯得尤為重要,對(duì)防治流域地表水體富營(yíng)養(yǎng)化具有重要意義。贛江是長(zhǎng)江八大支流之一, 起源于江西石城縣, 贛江主支于永修吳城鎮(zhèn)匯入鄱陽(yáng)湖[8]。贛江下游屬于中亞熱帶濕潤(rùn)季風(fēng)氣候, 氣候溫和濕潤(rùn), 降雨量充沛[9],根據(jù)該流域的地理?xiàng)l件與土地利用狀況, 選擇贛江下游水田、旱地、林地及其它地等四種類型的徑流小區(qū), 利用 SCS模型估算降雨徑流量, 結(jié)合降雨徑流中氮、磷含量估算其污染負(fù)荷, 為控制降雨徑流中氮、磷流失提供科學(xué)依據(jù)。

2 材料與方法

2.1 樣品采集

結(jié)合贛江下游不同土地利用方式及其特點(diǎn), 選擇四個(gè)典型徑流小區(qū)進(jìn)行取樣觀測(cè), 其降雨徑流水樣的采集點(diǎn)分別布置在贛江下游的蔣巷鎮(zhèn)(樣地類型為水田)、揚(yáng)子洲鎮(zhèn)(樣地類型為旱地)、南新鄉(xiāng)(樣地類型為林地)與南磯鄉(xiāng)(樣地類型為其它地), 于2010年5月至8月在不同季節(jié)遇較大降雨時(shí)在上述采樣點(diǎn)采集降雨徑流水樣四次(采樣時(shí)間依次為4月12日、5月8日、5月14日與6月20日), 獲取的降雨徑流水樣在經(jīng)過0.45 μm的尼龍膜過濾后裝入聚乙烯塑料瓶, 加入飽和HgC12固定, 在4 ℃下保存?zhèn)錅y(cè), 用來(lái)測(cè)量降雨徑流水體中的氮、磷濃度。

2.2 實(shí)驗(yàn)方法

降雨徑流中氮、磷的分析方法按照《湖泊富營(yíng)養(yǎng)化調(diào)查規(guī)范》(第二版)[10]中的方法進(jìn)行測(cè)定, 磷酸鹽采用鉬銻抗比色法, 氨氮采用納氏試劑比色法,硝態(tài)氮采用酚二磺酸比色法。其中, 溶解態(tài)無(wú)機(jī)氮(DIN)為硝態(tài)氮與氨氮之和, 溶解態(tài)無(wú)機(jī)磷(DIP)為溶解態(tài)的磷酸鹽。

贛江下游不同樣地類型在降雨徑流時(shí)攜帶的氮、磷負(fù)荷量可以按照降雨徑流法計(jì)算得到, 其計(jì)算公式如下[11-12]:

其中,W為農(nóng)業(yè)面源污染隨降雨徑流時(shí)攜帶的氮、磷負(fù)荷量(t·a–1);C為農(nóng)業(yè)面源污染產(chǎn)區(qū)不同土地利用類型降雨徑流中氮、磷濃度(mg·L–1)(一般采用多次實(shí)測(cè)平均值);Q為研究區(qū)間降雨徑流時(shí)產(chǎn)生的地表徑流量(mm)。

3 結(jié)果與討論

3.1 不同樣地的降雨徑流量

3.1.1 徑流量的估算過程

徑流曲線模型(Soil Conservation Service, 簡(jiǎn)稱SCS模型)是由美國(guó)農(nóng)業(yè)部土壤保持局發(fā)展起來(lái)的用于估算農(nóng)業(yè)流域降雨徑流量的流域水文模型, 具有所需資料簡(jiǎn)單易取等特點(diǎn), 能夠較為客觀地描述不同土地利用方式、土壤類型、前期的土壤含水量以及流域水文、氣象資料條件下的地表徑流形成過程, 對(duì)于小流域徑流預(yù)報(bào)具有較強(qiáng)的實(shí)用性和精確性。本研究利用SCS模型來(lái)估算贛江下游的降雨徑流量, SCS模型的降水徑流關(guān)系為[13]:

其中,Q為區(qū)域的降雨徑流量(mm);P為區(qū)域的降雨量(mm);S為降雨徑流開始后潛在的最大滯留量或稱為儲(chǔ)留系數(shù)(mm), 與流域的土壤類型、植被等因素有關(guān), 可通過徑流曲線數(shù)(CN)計(jì)算確定;Ia為初始截留量(mm), 指徑流開始前全部降雨損失, 包括表層低洼儲(chǔ)水、植被截獲、蒸發(fā)、入滲等, 其與S的關(guān)系為Ia=0.2S。

徑流曲線數(shù)(Runoff Curve Number, 簡(jiǎn)寫為CN)是反映降雨前流域特征綜合參數(shù), 它與流域前期的土壤濕潤(rùn)狀況(AMC)、植被、坡度、農(nóng)業(yè)耕作方式、水文條件、土地利用方式和土壤類型等因素有關(guān),取值范圍為 0—100[14], CN值越大, 其滲透量越小,產(chǎn)流能力越強(qiáng)[15]。為了反映土壤水分條件對(duì)CN值的影響, SCS模型根據(jù)前期降雨量大小將前期土壤的濕潤(rùn)狀況AMC(Antecedent Moisture Condition)分干旱、正常和濕潤(rùn)三級(jí), 即AMC I為干旱情況(產(chǎn)流條件最低), AMC II 為正常情況, AMC III為濕潤(rùn)情況(產(chǎn)流條件最高), 不同的前期土壤濕潤(rùn)狀況, 其CN值不同; 此外, 土壤類型按土壤水分滲透率可分為四類, 即為 A型土壤(透水性好, 產(chǎn)流能力最小,這類土壤主要是由深厚的、易大量排水的礫石組成)、B型土壤(較為透水, 產(chǎn)流能力為中等, 土壤質(zhì)地由中等細(xì)到中等粗, 土壤礫粒含量少)、C型土壤(較不透水, 其產(chǎn)流能力較好, 土層較淺, 土壤質(zhì)地常常有中到細(xì), 存在有阻礙水向下運(yùn)動(dòng)的土層)和D型土壤(接近于不透水, 其產(chǎn)流能力為最強(qiáng), 主要為含有高膨脹性的粘性土壤及幾乎不透水物質(zhì)的淺層土壤)四類。

根據(jù)贛江下游的土壤性質(zhì)(紅壤, 屬于粘壤土,即粘質(zhì)含量較高的土壤)和土地利用類型(本研究將贛江下游的土地利用類型分為水田、旱地、林地及其它地等四類), 確定該區(qū)域的土壤類型為C類, 即較不透水土壤, 故計(jì)算CN值時(shí)選擇C類土壤。參照SCS模型中提供的CN取值方法, 并結(jié)合其他學(xué)者的已有研究成果[16–19], 以及贛江下游地區(qū)的降雨徑流資料, 確定水文條件為正常含水量時(shí), 不同土地利用方式CN值, 此即CN2值, 當(dāng)土壤處于干旱或濕潤(rùn)狀況時(shí), 則需要進(jìn)行相應(yīng)地修正, 即利用正常狀態(tài)的CN2值分別計(jì)算干旱和濕潤(rùn)狀況下CN1值和CN3值。

3.1.2 徑流量的估算結(jié)果

降雨是農(nóng)業(yè)面源污染發(fā)生的源動(dòng)力, 對(duì)土壤中氮、磷等營(yíng)養(yǎng)鹽的流失具有較大的影響[20]。不同土地利用類型在降雨時(shí)均有產(chǎn)流過程, 但其特征因植被覆蓋和土壤類型等因素不同而異[20]。資料顯示,贛江下游多年平均年降水量約為1570 mm[18], 結(jié)合前述降雨徑流量的估算公式, 估算得到了贛江下游地區(qū)不同土地利用類型的年降雨徑流量, 其中水田的年降雨徑流量為962.42 mm, 旱地為552.16 mm,林地為428.75 mm, 其它地為682.45 mm, 可見贛江下游地區(qū)不同土地利用類型區(qū)域的年降雨徑流量差異相對(duì)比較明顯, 可能與不同土地利用類型所采用的農(nóng)業(yè)生產(chǎn)方式以及地表的覆蓋狀況等不同有關(guān),不同的農(nóng)業(yè)生產(chǎn)方式改變了地表的覆蓋情況, 而不同的地表覆蓋情況則影響著降雨的地表水文過程。此外, 水田、旱地、林地和其它地的降雨徑流系數(shù)分別為0.61、0.35、0.27和0.43, 反映出不同土地利用類型的徑流系數(shù)差異較大, 并以水田的徑流系數(shù)最大。

3.2 不同樣地的氮磷污染負(fù)荷

3.2.1 降雨徑流水體的氮磷含量

由于降雨徑流水樣采樣點(diǎn)的設(shè)置數(shù)量與覆蓋面、以及采樣頻次等因素的限制, 從而導(dǎo)致了本次降雨徑流水樣中氮、磷監(jiān)測(cè)結(jié)果不能完全代表整個(gè)贛江下游地區(qū)的年平均降雨徑流中的氮、磷濃度。因此, 本研究將氮、磷濃度的監(jiān)測(cè)結(jié)果近似認(rèn)為是贛江下游地區(qū)不同土地利用類型中降雨徑流中 DIP 與DIN濃度的年平均估算值。

監(jiān)測(cè)結(jié)果表明, 贛江下游地區(qū)不同樣地類型中的水田、旱地、林地及其它地降雨徑流水體中DIP平均濃度分別為0.131 mg·L–1、0.164 mg·L–1、0.108 mg·L–1、0.192 mg·L–1, DIN的平均濃度分別為1.312 mg·L–1、1.861 mg·L–1、0.426 mg·L–1、0.675 mg·L–1?？傮w上, 以水田和旱地降雨徑流水體中DIN與DIP的濃度較高, 可能與這些區(qū)域的農(nóng)業(yè)生產(chǎn)過程中施用了過多的肥料以及不同的地表覆蓋狀況有關(guān), 相比于林地和其它地, 水田和旱地在耕作時(shí)需要投入更多的肥料, 土壤中氮、磷的積累量大, 在降雨徑流時(shí)其輸出量也更大。此外, 水田和旱地頻繁的耕作造成地表植被覆蓋破壞, 從而造成植被對(duì)徑流水體中的氮、磷吸收下降。

3.2.2 降雨徑流氮磷的污染負(fù)荷

降雨徑流攜帶的氮、磷輸出與遷移方式是非常復(fù)雜的, 主要包括沉淀、沖起、吸附和解吸、轉(zhuǎn)化、衰減和生物吸收等過程。贛江下游地區(qū)降雨徑流攜帶的氮、磷等污染負(fù)荷主要與區(qū)域的土壤類型、土地利用類型等因素有關(guān), 而土地利用類型又是影響降雨徑流攜帶的氮、磷輸出負(fù)荷的重要因素。結(jié)合區(qū)域面源污染產(chǎn)區(qū)降雨徑流攜帶的氮、磷輸出負(fù)荷量的估算公式, 估算得到2010年贛江下游地區(qū)不同土地利用類型降雨徑流攜帶的DIP與DIN年負(fù)荷量,見圖1與圖2。

圖1 降雨徑流攜帶的DIP年負(fù)荷量Fig. 1 DIP loading of rainfall runoff in Ganjiang River downstream

從圖1與圖2可知, (1)通過估算區(qū)域降雨徑流量,并結(jié)合降雨徑流法測(cè)定的 2010年贛江下游地區(qū)降雨徑流攜帶的DIP年負(fù)荷總量為2389.63 t·a–1, 其中,不同土地利用類型的DIP年負(fù)荷量范圍為166.70—1071.65 t·a–1, 并以水田的DIP輸出負(fù)荷量最高, 達(dá)到1071.65 t·a–1, 旱地次之, 而以林地的DIP輸出負(fù)荷量最小, 只有166.70 t·a–1。(2)贛江下游地區(qū)降雨徑流攜帶的DIN年負(fù)荷總量為18913.06 t·a–1, 其中,不同土地利用類型在降雨徑流時(shí)攜帶的DIN輸出負(fù)荷量變化范圍為 657.53—10732.91 t·a–1, 同樣以水田的 DIN年輸出負(fù)荷量最高, 達(dá)到 10732.91 t·a–1,而以林地的DIN年負(fù)荷量最小。綜上可知, 不同土地利用類型中, 水田在降雨徑流時(shí)攜帶的氮、磷負(fù)荷量為最高, 其對(duì)贛江下游地區(qū)地表水體中氮、磷含量的貢獻(xiàn)率最大, 故而需重點(diǎn)關(guān)注水田對(duì)贛江下游地區(qū)地表水體富營(yíng)養(yǎng)化的貢獻(xiàn)。

圖2 降雨徑流攜帶的DIN年負(fù)荷量Fig. 2 DIN loading of rainfall runoff in Ganjiang River downstream

4 結(jié)論

土地利用類型是影響降雨徑流量的重要因素,贛江下游不同土地利用類型降雨徑流量差異較大,其中水田的年降雨徑流量為 962.42 mm, 旱地為552.16 mm, 林地為428.75 mm, 其它地為682.45 mm,可能與不同土地利用類型的所采用的農(nóng)業(yè)生產(chǎn)方式以及地表覆蓋狀況等不同有關(guān)。

降雨徑流攜帶的氮、磷輸出負(fù)荷同樣受到不同土地利用類型的影響, 贛江下游不同樣地類型 DIP與DIN年負(fù)荷量變化范圍分別為166.70—1071.65 t·a–1與 657.53—10732.91 t·a–1, 均以水田的輸出負(fù)荷量最高, 而以林地的輸出負(fù)荷最小, 這主要與不同土地利用類型農(nóng)業(yè)生產(chǎn)方式以及施肥量與頻次等因素有關(guān)。

[1] 宋澤芬, 王克勤, 孫孝龍, 等. 澄江尖山河小流域不同土地利用類型地表徑流氮、磷的流失特征[J]. 環(huán)境科學(xué)研究, 2008, 21(4): 109–114.

[2] 朱松, 方沛南, 藍(lán)雪春. 降雨徑流污染研究綜述[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2009, 25(12): 240–245.

[3] 張健, 陳鳳, 濮勵(lì)杰, 等. 經(jīng)濟(jì)快速增長(zhǎng)區(qū)土地利用變化對(duì)土壤質(zhì)量的影響研究[J]. 環(huán)境科學(xué)研究, 2007, 20(5): 99–104.

[4] 傅伯杰, 馬克明, 周華峰, 等. 黃土丘陵區(qū)土地利用結(jié)構(gòu)對(duì)土壤養(yǎng)分分布的影響[J]. 科學(xué)通報(bào), 1998, 43(22): 2444–2447.

[5] ARHEIMER B, LIDN R. Nitrogen and phosphorus concentrations from agricultural catchments-influence of spatial and temporal variables [J]. Journal of Hydrology, 2000, 227(1): 140–159.

[6] LIANG T, WANG H, KUNG H T, et al. Agriculture land-use effects on nutrient losses in west Tiaoxi watershed, China [J]. Agriculture Journals, 2004, 40(6): 1499–1510.

[7] 張紅愛, 張煥朝, 鐘萍. 太湖地區(qū)典型水稻土稻—麥輪作地表徑流中磷的變動(dòng)規(guī)律[J]. 生態(tài)科學(xué), 2008, 27(1): 17–23.

[8] 金顯文, 楊威, 丁建華, 等. 贛江下游及其支流輪蟲群落結(jié)構(gòu)及水質(zhì)生物學(xué)評(píng)價(jià)[J]. 生態(tài)科學(xué), 2014, 33(1): 99–105.

[9] 向速林, 王全金, 徐劉凱, 等. 贛江尾閭區(qū)農(nóng)田溝渠徑流中氮磷負(fù)荷與遷移特征[J]. 河南農(nóng)業(yè)科學(xué), 2012, 41(3): 72–74.

[10] 金相燦, 屠清瑛. 湖泊富營(yíng)養(yǎng)化調(diào)查規(guī)范[M]. 北京: 中國(guó)環(huán)境科學(xué)出版社, 1990.

[11] 金相燦. 中國(guó)湖泊環(huán)境(第二冊(cè))[M]. 北京: 海洋出版社, 1995.

[12] 屠清瑛. 巢湖富營(yíng)養(yǎng)化研究[M]. 合肥: 中國(guó)科學(xué)技術(shù)大學(xué)出版社, 1990.

[13] 史志華, 蔡崇法, 丁樹文. 基于 GIS的漢江中下游農(nóng)業(yè)面源氮磷負(fù)荷研究[J]. 環(huán)境科學(xué)學(xué)報(bào), 2002, 22(4): 473–477.

[14] 黃金良, 洪華生, 張珞平, 等. 基于GIS的九龍江流域農(nóng)業(yè)非點(diǎn)源點(diǎn)磷負(fù)荷估算研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2004, 23(5): 866–871.

[15] 羅鵬, 宋星原. 基于格柵的分布式SCS產(chǎn)流模型研究[J].水土保持通報(bào), 2010, 30(4): 138–142.

[16] 王曉輝. 巢湖流域非點(diǎn)源N、P污染排放負(fù)荷估算[D]. 合肥: 合肥工業(yè)大學(xué), 2006.

[17] 劉哲. 福州第二水源山仔水庫(kù)的農(nóng)業(yè)非點(diǎn)源污染負(fù)荷研究[D]. 福州: 福建師范大學(xué), 2006.

[18] 王曉鴻, 鄢幫有, 吳國(guó)琛. 山江湖工程[M]. 北京: 科學(xué)出版社, 2006.

[19] POLYAKOV V, FARES A, KUBO D, et al. Evaluation of a non-point source pollution model, AnnAGNPS, in a tropical watershed [J]. Environmental Modelling & Software, 2007, 22: 1617–1627.

[20] 溫海廣, 周勁風(fēng), 李明, 等. 流溪河水庫(kù)流域非點(diǎn)源溶解態(tài)氮磷污染負(fù)荷估算[J]. 環(huán)境科學(xué)研究, 2011, 24(4): 387–394.

Nitrogen and phosphorus pollution loads from the different land use types in the Ganjiang River downstream

XIANG Sulin1,2,*, TAO Shuping1, WANG Fengwu1
1.Department of Environment Engineering,East China JiaoTong University,Nanchang330013,China2.Key Laboratory of Poyang Lake Environment and Resource Utilization of MOE,Nanchang University,Nanchang330047,China

Several runoff units of various land use areas are selected to monitor nitrogen and phosphorus contents in Ganjiang River downstream. Using SCS runoff curve, runoff of each plot is calculated. Considering the actual concentration of nutrient, loading of nutrient of surface runoff is estimated. The results show that phosphorus and nitrogen outputs in runoff were significant in different land use patterns, and high nitrogen and phosphorus loads are associated with paddy land use, which are 10732.91 t·a–1and 1071.65 t·a–1, respectively. The amount of nutrient loss from woodland is the lowest, which are 657.53 t·a–1and 166.70 t·a–1. Being carrier of loss nutrients, the paddy land use is thought to be the controlling focus in Ganjiang River downstream.

Ganjiang River downstream; Rainfall runoff; Pollution load

10.14108/j.cnki.1008-8873.2016.02.010

X144

A

1008-8873(2016)02-062-04

2014-10-15;

2015-01-12

江西省教育廳科研項(xiàng)目(GJJ150540); 江西省自然科學(xué)基金項(xiàng)目(20114BAB213020)

向速林(1978—), 男, 江西人, 博士, 副教授, 研究方向?yàn)樗Y源與環(huán)境, E-mail: slxiang2001@163.com

*通信作者:向速林