楊佳磊,張 瑞

(江蘇海洋大學 海洋技術與測繪學院，江蘇 連云港 222005)

0 引言

隨著長三角地區(qū)經濟的快速發(fā)展，太湖流域水環(huán)境污染問題日益嚴重，引起了人們的廣泛關注[1-2]。目前工業(yè)點源污染逐步得到有效控制，而非點源污染已經成為當前水環(huán)境污染的最主要來源[3]。非點源污染的影響因素較多，土地利用被認為是影響非點源污染的關鍵因素[4-6]。土地利用方式與非點源污染密切相關，其中尤以耕地和林地與其關系最為密切。研究表明，農田排水是造成水體富營養(yǎng)化的主要原因[7]，林地則對非點源污染具有顯著的截留作用[8]，土地利用格局的優(yōu)化配置對非點源污染的控制有著重要的意義[9]。

太湖流域經濟發(fā)達，人口密集，河網復雜。快速的城鎮(zhèn)化進程導致土地利用類型發(fā)生了較大轉變[10]。同時農業(yè)生產中化肥的大量使用，導致流域生態(tài)環(huán)境污染負荷日益加重，直接或間接地影響流域的水環(huán)境質量[11-13]。在流域非點源污染控制的計算方法中，輸出負荷系數算法因其對數據的要求相對較低，計算結果精度較高，在較大尺度的流域污染負荷估算中有著較好的實際應用效果，目前已在國內外得到廣泛應用[14]。荊延德等[15]、李政道等[16]、郝桂珍等[17]分別對南四湖、紅楓湖、清水河流域的非點源污染負荷進行了估算及變化規(guī)律分析；李思思等[18]引入產污因子和截留因子校正非點源氮磷流失的空間格局；龐樹江等[19]則通過回歸樹模型確定了主要污染源；Matias等[20]將輸出系數模型應用于干旱地區(qū)，探究了磷來源及從陸地到水體的輸出系數；Frink等[21]整理了美國東北部地區(qū)不同土地利用的N和P輸出系數，并驗證了模型的精度。

本文利用輸出系數法，通過對太湖流域1980—2018年間土地利用類型變化及同期非點源氮磷負荷的影響進行研究，探討土地利用變化對非點源氮磷負荷的影響機制，從而為太湖流域水環(huán)境治理提供借鑒。

1 材料與方法

1.1 研究區(qū)概況

太湖流域北濱長江，南接錢塘江，東臨東海，西以天目山、茅山等山區(qū)為界，流域總面積約3.68萬km2。地形呈周邊高、中間低的碟狀分布，流域內最大高程500 m左右，西南低山丘陵區(qū)高程在12～32 m之間，中部平原區(qū)域高程一般低于5 m，沿江及濱海平原區(qū)域高程一般在5～10 m。太湖流域內河網密布，上游水系有苕溪、南溪和洮滆水系等，經望虞河、太浦河等水閘流出，與長江、黃浦江水道相通。流域內土地類型以耕地為主，主要農作物是水稻、小麥及其他經濟作物。該研究區(qū)域地理位置如圖1所示。

圖1 太湖流域地理位置

1.2 數據來源及處理

太湖流域邊界數據及1980年、1990年、2000年、2010年、2018年土地利用數據均來自中國資源科學與數據中心(https://www.resdc.cn)。其中土地利用數據分辨率為30 m，在ArcGIS軟件中按照太湖流域邊界裁剪而來。按照全國土地分類標準，將研究區(qū)域的土地利用類型分為6類：耕地、林地、草地、水域、建設用地和未利用地。

1.3 數據分析方法

1.3.1 動態(tài)度分析方法 土地利用相對動態(tài)度不僅可以直觀地反映某一土地利用類型相對于本身的變化速度和劇烈程度，也可以反映各土地利用類型的變化在區(qū)域變化中的貢獻大小，衡量某土地利用類型的變化是否成為區(qū)域變化的主要類型。其計算公式為

(1)

式中：Bi為第i種土地利用相對動態(tài)度；Pi為研究初期第i種土地利用類型占研究區(qū)土地總面積的比率；S(i,t1)和S(i,t2)分別表示基期、現期某一土地利用類型面積。

1.3.2 輸出系數法 輸出系數模型主要是通過污染源的輸出系數來估算污染負荷，該模型對大尺度非點源污染負荷的測算有著較高的精度[22-23]。輸出系數模型表達式為

(2)

式中：Lj為污染物j在該流域的總負荷量，單位為kg/km2；i為流域內土地利用類型，包括耕地、林地、草地、水域、建設用地和未利用地共計6類；Eij為污染物j在第i種土地利用類型下的輸出系數，單位kg/(km2·a)；Ai為第i類土地利用類型的面積，單位km2。氮、磷是影響太湖流域表水質的主要致污因子，因此本研究選取氮、磷作為污染物。

同種土地利用類型的氮、磷輸出系數具有巨大的時空變異性，這不僅與降水量、徑流量和施肥量有關，而且與坡度、距河遠近、土壤類型等因素息息相關。不同土地利用類型的氮、磷輸出系數隨著降雨量或徑流量的增加而呈現出增加的趨勢[24]。因此，確定合理的輸出系數是構建模型的關鍵。國內外對輸出系數模型作了大量研究，相關的方法體系也較為成熟。本研究參照太湖流域已有的相關研究成果[25-29]，給出不同土地利用類型的污染輸出系數，具體如表1所示。

表1 太湖流域不同土地利用類型氮、磷輸出系數

2 結果與分析

2.1 土地利用變化特征

2.1.1 土地利用結構變化 從1980年、1990年、2000年、2010年、2018年5個時期土地利用變化(見圖2)來看，太湖流域各個階段變化最大的土地利用類型主要為耕地和建設用地，耕地表現為減少，建設用地表現為增加，這與改革開放以來長三角地區(qū)工業(yè)發(fā)展和城鎮(zhèn)化進程有很大關系。從土地利用結構(見圖3)來看，太湖流域主要土地利用類型為耕地，5個時期均占比45%以上，草地和未利用地面積占比最小，不足流域總面積的1%。耕地的面積占比由1980年的65.66%減少到2018年的45.15%，而建設用地的面積占比由8.51%增加到27.7%。

圖2 1980年、1990年、2000年、2010年和2018年太湖流域土地利用類型變化

圖3 1980年、1990年、2000年、2010年和2018年各土地利用類型占比情況

2.1.2 動態(tài)度變化 基于各時期不同土地利用類型的面積，運用式(1)計算得到各個時期不同土地利用類型的相對動態(tài)度，結果見表2。由表2可知，4個時期相對動態(tài)度絕對值最高的土地利用類型均是耕地，相對動態(tài)度絕對值分別為0.149 5，0.529 2，1.088 1和0.559 6。各個時期土地利用變化劇烈程度不同，1980—1990年土地利用變化較為緩和，2000—2010年土地利用變化最劇烈，6大地類相對動態(tài)度總和達到2.199 7，其中變化最明顯的土地利用類型為耕地和建設用地。2000年以來，太湖流域經歷了高速城鎮(zhèn)化過程，已經成為我國城鎮(zhèn)化水平較高的地區(qū)之一[30]。高速城鎮(zhèn)化過程伴隨著建設用地的快速擴張。2000—2010年，太湖流域各大城市建成區(qū)面積都有不同程度擴大。如上海建成區(qū)面積由2000年的857.94 km2增加到2010年的1 514.52 km2，城區(qū)面積擴大了接近1倍；蘇州建成區(qū)面積由2000年的137.81 km2增加到2010年的494.03 km2，城區(qū)面積擴大了2.58倍；2010年常州建成區(qū)面積是2000年的2.46倍[31]。而耕地變化劇烈，一方面是由于建設用地需求增加，尤其是在江蘇省，蘇南鄉(xiāng)鎮(zhèn)企業(yè)的快速發(fā)展占用了大量寶貴的耕地資源；另一方面，在國家取消農民糧食定購任務的政策出臺后，農業(yè)結構發(fā)生了很大變化。

表2 1980—2018年各土地利用類型相對動態(tài)度分析

從4個時期太湖流域土地利用轉移矩陣(表3～表6)可以看出，2000—2010年土地轉變最大，約占總量的51.07%。耕地和建設用地轉變數量最多，分別占變化總量的49.15%和44.89%；未利用地、耕地和建設用地參與其他土地利用轉變率最高，轉變量分別占2000年未利用地、耕地和建設用地的348.08%，69.73%和18.13%。太湖流域土地利用轉變過程中，耕地主要轉變?yōu)榻ㄔO用地、水域和林地；建設用地主要轉變?yōu)楦兀凰蚝土值刂饕D變?yōu)楦亍?/p>

表3 1980—1990年太湖流域土地利用轉移矩陣

表4 1990—2000年太湖流域土地利用轉移矩陣

表5 2000—2010年太湖流域土地利用轉移矩陣

表6 2010—2018年太湖流域土地利用轉移矩陣

2.2 不同土地利用方式下非點源氮、磷輸出年際變化特征

太湖流域1980年、1990年、2000年和2018年TN負荷量總體呈逐年上升趨勢，從1980年的6.95萬t/a上升至2018年的7.51萬t/a，增幅達到8.06%。耕地是TN負荷的最大貢獻源，但負荷量逐年減少。林地、草地和未利用地TN負荷貢獻極少且年際變化不大。建設用地和水域的TN負荷逐年上升。

圖4表明，太湖流域1980年、1990年、2000年、2018年TP負荷量也呈上升的趨勢，從1980年的0.37萬t/a上升至2018年的0.96萬t/a，增幅達到159.46%。林地、草地、水域和未利用地TP負荷貢獻較少且年際變化不大。耕地是TP負荷的最大貢獻源，2000年以前耕地TP負荷變化不大，2000年以后起伏明顯，耕地TP負荷總體呈現上升趨勢。建設用地TP負荷2000年以前上升緩慢，2000年以后上升明顯。

a 總氮負荷量

2.3 土地利用結構轉變對非點源氮、磷輸出的影響

為了定性定量描述土地利用變化過程對非點源氮、磷輸出的影響，通過相對系數(兩個時期土地利用輸出系數差)與土地利用變化量之間的關系計算得出非點源氮磷過程量。該過程量能真實反映不同土地利用類型變化給非點源氮磷內部輸出所帶來的影響[32]。

圖5表示某一土地利用類型發(fā)生轉變以及其他土地利用類型轉變?yōu)樵撏恋乩妙愋瓦^程中，該土地利用類型上所發(fā)生的非點源氮磷負荷的輸出凈量。正軸表示起到增加非點源氮磷的輸出凈量，負軸表示起到削減非點源氮磷的輸出凈量。

a 耕地TN輸出凈量

4個時期TN負荷凈量總量分別為2 365.26，2 057.31，3 664.25和2 121.92 t，對TN負荷凈量的貢獻較大的地類均是耕地和建設用地，貢獻較小的地類均是草地和未利用地。其中耕地貢獻占比最高值為1980—1990年間的47.08%，建設用地貢獻占比最高值為2000—2010年間的35.88%。

4個時期TP負荷凈量總量分別為89.19，204.45，1 513.62和1 113.91 t，其中1980—1990年間對TN負荷凈量的貢獻較大的地類為耕地和水域，而其余3個時期均是耕地和建設用地，耕地貢獻占比最高值為2000—2010年間的47.95%，建設用地貢獻占比最高值為2010—2018年間的40.95%。4個時期對TP負荷凈量貢獻較小的地類均是草地和未利用地。

不同土地利用轉變對TN和TP負荷凈量輸出總量增長或消減的貢獻能力不同。1980—1990年建設用地轉變?yōu)楦貢r對TN的表現為增加，耕地轉變?yōu)榱值貢r對TN的表現為消減；水域轉變?yōu)楦貢r對TP的表現為增加，耕地轉變?yōu)榱值貢r對TP的表現為消減。1990—2000年建設用地轉變?yōu)楦貢r對TN的表現為增加，耕地轉變?yōu)榱值貢r對TN的表現為消減；建設用地轉變?yōu)楦貢r對TP的表現為增加，耕地轉變?yōu)榱值貢r對TP的表現為消減。2000—2010年建設用地轉變?yōu)楦貢r對TN的表現為增加，建設用地轉變?yōu)榱值貢r對TN的表現為消減；建設用地轉變?yōu)楦貢r對TP的表現為增加，耕地轉變?yōu)榱值貢r對TP的表現為消減。2010—2018年建設用地轉變?yōu)楦貢r對TN的表現為增加，耕地轉變?yōu)榱值貢r對TN的表現為消減；耕地轉變?yōu)榻ㄔO用地時對TP的表現為增加，建設用地轉變?yōu)楦貢r對TP的表現為消減。對于TN而言，建設用地轉變?yōu)楦貢r表現為增加，耕地轉變?yōu)榱值亍⒔ㄔO用地轉變?yōu)榱值貢r表現為削減；對于TP而言，1980—1990年、1990—2000年、2000—2010年間建設用地轉變?yōu)楦?、水域轉變?yōu)楦貢r表現為增加，耕地轉變?yōu)榱值貢r表現為削減，而2010—2018年間耕地轉變?yōu)榻ㄔO用地時表現為增加，建設用地轉變?yōu)榱值貢r表現為削減。

3 結論

本文以太湖流域1980年、1990年、2000年、2010年、2018年土地利用數據為基礎，并通過輸出系數模型分析太湖流域土地利用變化對非點源氮磷輸出的影響，得出以下結論。

(1)耕地和建設用地是太湖流域近40年間面積變化最劇烈的土地利用類型；土地利用方式以耕地向建設用地、水域和林地轉變?yōu)橹鳌?/p>

(2)太湖流域TN和TP負荷量總體逐年上升，TN從1980年的6.95萬t到2018年的7.51萬t，TP從1980年的0.37萬t到2018年的0.96萬t；耕地是TN和TP負荷的最大貢獻源，5個時期平均占比分別達到67.12%和65.05%。

(3)不同土地利用方式轉變對TN和TP負荷凈量輸出總量增長或消減的貢獻能力不同。對于TN而言，建設用地向耕地轉變時表現為增加，耕地向林地、建設用地向林地轉變時表現為削減；對于TP而言，1980—1990年、1990—2000年、2000—2010年間建設用地向耕地轉變、水域向耕地轉變時表現為增加，耕地向林地轉變時表現為削減，而2010—2018年間耕地向建設用地轉變時表現為增加，建設用地向林地轉變時表現為削減。