姜 暢，劉鴻雁，*，陳 竹，饒 程 ，李政道，楊通銓，吳 攀

（1.貴州大學(xué)資源與環(huán)境工程學(xué)院，貴陽(yáng) 550025；2.貴州大學(xué)農(nóng)學(xué)院，貴陽(yáng) 550025；3.貴陽(yáng)市兩湖一庫(kù)管理局，貴州 清鎮(zhèn) 551400）

土地利用變化涉及區(qū)域人口、資源、環(huán)境與發(fā)展方面的核心問(wèn)題。研究表明土地利用變化是影響水質(zhì)的重要原因[1]，其對(duì)流域面源污染的影響是導(dǎo)致水環(huán)境惡化的重要因素之一[2-4]。紅楓湖是貴陽(yáng)市重要的飲用水源保護(hù)地，流域土地利用變化與水環(huán)境安全密切相關(guān)。從1990年起，土地利用變化便成為全世界資源與環(huán)境領(lǐng)域的研究熱點(diǎn)[5-6]，眾多學(xué)者利用高光譜遙感、微波遙感、高分辨率影像和多源信息復(fù)合、遙感圖像自動(dòng)識(shí)別與分類、3S集成等技術(shù)[7]，針對(duì)不同時(shí)空土地利用的變化、特征、趨勢(shì)開(kāi)展了大范圍的研究工作。

在土地利用對(duì)水質(zhì)影響的研究中，蔡宏等[8]、Nash等[9]、曹芳芳等[10]或通過(guò)子流域土地利用類型為分析樣本，或以監(jiān)測(cè)斷面為中心建立緩沖區(qū)等手段從空間上研究不同距離土地利用類型對(duì)水質(zhì)的相關(guān)性。郭青海等[11]、楊峰等[3]利用多元逐步回歸模型分析了土地利用類型與流域湖泊水質(zhì)間的相關(guān)關(guān)系，從景觀生態(tài)學(xué)的角度闡述了不同土地利用類型對(duì)營(yíng)養(yǎng)元素輸出的源匯作用，對(duì)貢獻(xiàn)率高的土地利用類型提出了面源污染的分類控制措施。孫金華等[12]通過(guò)子流域歸類分析了土地利用空間格局差異對(duì)水質(zhì)的影響，表明土地利用比例結(jié)構(gòu)和分布格局對(duì)水質(zhì)亦具有重要影響。廖建榮[13]采用冗余分析（RDA）方法揭示土地利用類型與河流水質(zhì)的空間依賴性。因此，通過(guò)量化流域土地利用類型的面積，再與水質(zhì)進(jìn)行關(guān)聯(lián)分析是研究土地利用類型與污染物質(zhì)內(nèi)在聯(lián)系的重要手段[14]。紅楓湖是貴陽(yáng)市的“三大水缸”之一，水質(zhì)質(zhì)量直接關(guān)乎區(qū)域用水安全。近年來(lái)已有較多關(guān)于紅楓湖流域的研究，但流域土地利用變化對(duì)水質(zhì)影響的研究仍鮮有報(bào)道。為明確影響水質(zhì)營(yíng)養(yǎng)指標(biāo)的主要土地利用類型，本研究選取紅楓湖準(zhǔn)保護(hù)區(qū)流域?yàn)檠芯繀^(qū)，以7個(gè)水質(zhì)監(jiān)測(cè)斷面為中心建立緩沖區(qū)，采用Spearman秩相關(guān)和景觀生態(tài)學(xué)理論探究紅楓湖流域15年間土地利用變化與水質(zhì)的響應(yīng)。結(jié)合研究區(qū)實(shí)際情況，文章從空間尺度變化的角度出發(fā)，以更大的緩沖半徑進(jìn)行了多個(gè)連續(xù)尺度的研究，分析不同距離的土地利用類型對(duì)中心點(diǎn)水質(zhì)的影響，使分析結(jié)果更為可靠。該結(jié)果可為流域面源污染的控制提供依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

紅楓湖流域地處貴州高原中部，長(zhǎng)江二級(jí)支流貓?zhí)拥纳嫌螀^(qū)，是貴州省最大的高原人工湖泊之一，也是喀斯特高原深谷型湖泊。該流域地處亞熱帶濕潤(rùn)季風(fēng)區(qū)，海拔1250 m，雨量豐沛，多年平均降水量達(dá)1200 mm[15]。流域位于東經(jīng) 106°19′～106°28′，北緯 26°26′～26°35′，流域面積 1596 km2[16]，湖域面積 57.2 km2，總庫(kù)容6.01億m3，湖內(nèi)小島零星散布，多已開(kāi)發(fā)為旅游景點(diǎn)。經(jīng)過(guò)60多年的發(fā)展，現(xiàn)已發(fā)展成為具有飲用水源、工農(nóng)業(yè)用水、旅游、發(fā)電等多功能的水庫(kù)，供水量占貴陽(yáng)市主城區(qū)用水量的70%[16]。流域內(nèi)碳酸鹽巖在垂直和水平分布上均占有較大的比例，巖溶發(fā)育強(qiáng)烈，巖溶地貌成為流域內(nèi)地貌類型的主要特征。從2003年到2007年紅楓湖水質(zhì)急劇惡化[16]，從Ⅲ類降到劣Ⅴ類，2008年開(kāi)始經(jīng)過(guò)6年的集中治理后水質(zhì)逐步好轉(zhuǎn)。

1.2 數(shù)據(jù)獲取

本研究采用條帶號(hào)與行編號(hào)分別為127、42的Landsat TM遙感影像圖，結(jié)合貴陽(yáng)市兩湖一庫(kù)管理局提供的地形邊界資料得到研究區(qū)流域遙感影像圖，具體包含2001年Landsat4-5 TM衛(wèi)星數(shù)據(jù)（時(shí)相4月13日，分辨率30 m）、2008年Landsat4-5 TM衛(wèi)星數(shù)據(jù)（時(shí)相4月8日，分辨率30 m）、2016年Landsat 8 OLI_TRIS衛(wèi)星數(shù)據(jù)（時(shí)相2月10日，分辨率30 m）共三期遙感影像作為基礎(chǔ)數(shù)據(jù)，數(shù)據(jù)來(lái)源于中科院遙感所（http：//ids.ceode.ac.cn/）與地理空間數(shù)據(jù)云（http：//www.gscloud.cn/），研究區(qū)域選擇為飲用水源保護(hù)區(qū)，含一級(jí)、二級(jí)保護(hù)區(qū)和準(zhǔn)保護(hù)區(qū)。水質(zhì)監(jiān)測(cè)點(diǎn)選取大壩、腰洞、花魚(yú)洞、后午、偏山寨、三岔河、焦家橋等7個(gè)監(jiān)測(cè)斷面，時(shí)間分別為2001、2008年和2016年，數(shù)據(jù)來(lái)源于貴陽(yáng)市兩湖一庫(kù)管理局環(huán)境監(jiān)測(cè)站。選取總磷（TP）、氨氮（NH3-N）、化學(xué)需氧量（COD）、溶解氧（DO）4個(gè)指標(biāo)進(jìn)行分析，監(jiān)測(cè)點(diǎn)分布見(jiàn)圖1。

1.3 土地利用類型及水質(zhì)劃分

參照土地利用類型的分類方法[17-19]，可將土地劃分為6個(gè)一類地和25個(gè)二類地。在分辨率30 m精度的圖形中，二類地?zé)o法得到有效區(qū)分，為了控制樣本的分類精度，結(jié)合分類方法和參考相關(guān)文獻(xiàn)，最終確定土地利用類型為耕地、林地、草地、建設(shè)用地、未利用地、水域。對(duì)遙感影像圖進(jìn)行輻射定標(biāo)、大氣校正和裁剪等預(yù)處理后進(jìn)行分類。通過(guò)對(duì)多光譜不同波段的對(duì)比研究后，最終對(duì)2001年和2008年的遙感影像采用543波段進(jìn)行RGB合成，對(duì)2016年的遙感影像采用654波段進(jìn)行RGB合成。在ENVI 5.2的支持下，采用支持向量機(jī)的監(jiān)督分類方法，結(jié)合Google Earth和實(shí)地調(diào)查對(duì)分類結(jié)果進(jìn)行人工修正，對(duì)結(jié)果進(jìn)行分類處理后，利用ArcGIS 10.2生成滿足高精度要求的土地利用類型圖。通過(guò)ENVI 5.2中Region of Interest Tool計(jì)算樣本的可分離系數(shù)（Separability），同時(shí)通過(guò)Confusion Matrix工具對(duì)2001、2008年和2016年的影像處理成果執(zhí)行精度評(píng)估?？煞蛛x系數(shù)的標(biāo)準(zhǔn)是大于1.8屬于合格樣本，超過(guò)1.9說(shuō)明樣本之間可分離性好，小于1.8則需要編輯樣本或者重新選擇樣本，若小于1，考慮將兩類樣本合成一類樣本?？傮w分類精度（OA）達(dá)到90%則為優(yōu)。Kappa指數(shù)可分為5個(gè)等級(jí)：0～0.20為幾乎不一致、0.21～0.40為較低一致、0.41～0.60 為中度一致、0.61～0.80 為高度一致、0.81～1為幾乎完全一致。根據(jù)各年份水質(zhì)狀況的不同對(duì)水域進(jìn)行劃分，其中Ⅰ、Ⅱ類水劃為一類，水質(zhì)為優(yōu)，Ⅲ類水劃為一類，水質(zhì)為較好，Ⅳ類水劃為一類，水質(zhì)為較差，Ⅴ、劣Ⅴ類劃為一類，水質(zhì)為差。

1.4 土地利用變化分析

土地利用變化幅度是從面積上考察土地利用類型的變化程度，重點(diǎn)描述不同土地利用類型在不同時(shí)間內(nèi)總量上的差異[20]。文章對(duì)土地利用類型的面積進(jìn)行動(dòng)態(tài)監(jiān)測(cè)，利用轉(zhuǎn)移矩陣[5-6]生成2001—2008年、2008—2016年的土地利用轉(zhuǎn)移矩陣，對(duì)紅楓湖流域15年間土地利用類型進(jìn)行轉(zhuǎn)換幅度和轉(zhuǎn)換方向分析。

1.5 研究尺度確定

建立緩沖區(qū)有利于從空間上分析不同距離的土地利用類型對(duì)中心點(diǎn)水質(zhì)的影響，該方法易于量化不同土地利用類型對(duì)水質(zhì)的影響范圍和影響程度[8]。文章以7個(gè)監(jiān)測(cè)斷面為圓心，準(zhǔn)保護(hù)區(qū)為邊界，以不同尺度量化分析土地利用類型與水質(zhì)的關(guān)系。緩沖半徑上，由于不同于城市邊界范圍緩沖區(qū)域的劃分，因此選用1、1.5、2 km和2.5 km尺度略大的緩沖半徑，在ArcGIS 10.2下生成緩沖區(qū)域（圖1），計(jì)算各土地利用類型的面積。

圖1 監(jiān)測(cè)斷面點(diǎn)及緩沖區(qū)域分布圖Figure1 Picture of monitoring points and buffer zone

1.6 土地利用與水質(zhì)的關(guān)聯(lián)分析

本文以7個(gè)監(jiān)測(cè)斷面為樣本，通過(guò)ArcGIS 10.2生成大壩、腰洞、花魚(yú)洞、后午、偏山寨、三岔河、焦家橋監(jiān)測(cè)斷面的圓形緩沖區(qū)（圖1），以1、1.5、2 km和2.5 km為緩沖尺度，分別將2001、2008年和2016年不同緩沖范圍內(nèi)土地利用類型的面積（林地、耕地、草地、建設(shè)用地、未利用地）和對(duì)應(yīng)水質(zhì)數(shù)據(jù)（TP、NH3-N、COD、DO）進(jìn)行關(guān)聯(lián)，利用 SPSS 17.0 進(jìn)行 Spearman相關(guān)性分析，探究15年間土地利用變化與水質(zhì)的響應(yīng)。根據(jù)景觀生態(tài)學(xué)“源匯”理論，對(duì)不同土地利用類型與水質(zhì)的相關(guān)性進(jìn)行“源匯”分析。不同地類對(duì)面源污染起到了正向推動(dòng)作用即為“源”景觀，相反，對(duì)面源污染起到了負(fù)向滯緩作用即為“匯”景觀[2，21]。

2 結(jié)果與討論

2.1 土地利用變化分析

2.1.1 土地利用分類精度評(píng)估

利用ENVI 5.2生成可分離系數(shù)（Separability）、總體分類精度（OA）、Kappa指數(shù)三個(gè)指標(biāo)。通過(guò)對(duì)紅楓湖流域2001、2008年和2016年的遙感影像分類進(jìn)行精度評(píng)估（表1）可知，分類精度整體較好，符合解譯要求。

2.1.2 土地利用類型分布

利用 ENVI 5.2與 ArcGIS 10.2生成 2001、2008年和2016年的土地利用類型及水質(zhì)變化圖（圖2）。利用Arc GIS 10.2的空間數(shù)據(jù)分析功能，對(duì)各個(gè)土地利用類型的面積進(jìn)行量化統(tǒng)計(jì)（表2）。15年間紅楓湖流域的總體特征：耕地和未利用地面積呈現(xiàn)逐年下降趨勢(shì)，林地面積先上升后下降，草地、建設(shè)用地面積持續(xù)上升，水域面積在小范圍內(nèi)波動(dòng)，變化比例較小。流域土地利用類型以林地、耕地為主，2016年分別占總面積的33.70%和29.67%。從2001年到2016年，林地面積減少了10.68 km2，減幅達(dá)7.30%，建設(shè)用地增加了11.41 km2，增幅最高，達(dá)36.11%，這也與經(jīng)濟(jì)發(fā)展和城鎮(zhèn)化建設(shè)的推進(jìn)相符合。

表1 土地利用類型分類精度Table1 Classification accuracy of land-use type

圖2 2001、2008年和2016年紅楓湖流域土地利用類型及水質(zhì)變化圖Figure2 Picture of land-use type and water quality change in Hongfeng basin in 2001,2008 and 2016

2.1.3 土地利用變化趨勢(shì)

分別對(duì)2001、2008年和2016年的土地利用類型圖進(jìn)行動(dòng)態(tài)監(jiān)測(cè)操作，利用ENVI 5.2的轉(zhuǎn)移矩陣分析得到2001—2008年、2008—2016年土地利用轉(zhuǎn)移矩陣（表3）。從2001年到2016年耕地持續(xù)減少，累計(jì)減少面積達(dá)8.61 km2。從2001年到2008年，耕地主要流出方向是建設(shè)用地和林地，兩者面積在7年間均有所增加。除耕地的貢獻(xiàn)外，未利用地也主要流入建設(shè)用地和林地，在2008年，未利用地流入建設(shè)用地和林地的面積共計(jì)1.34 km2，占未利用地整體變化量的68%。從2008年到2016年耕地的主要流出方向?yàn)榻ㄔO(shè)用地和草地。從2001年到2008年，林地面積有所增加，但從2001年到2016年，總體呈現(xiàn)下降趨勢(shì)，2016年林地面積下降了12.26 km2（表3），下降總量與其他地類相比最多，轉(zhuǎn)出的林地主要流入了建設(shè)用地和草地。綜上，耕地、林地和未利用地是建設(shè)用地和草地面積增加的主要來(lái)源，說(shuō)明15年間流域土地利用變化主要朝建設(shè)用地和草地增加。

表2 2001、2008年和2016年研究區(qū)土地利用類型Table2 Land-use type of the research area in 2001,2008 and 2016

表3 2001—2008年和2008—2016年研究區(qū)土地利用轉(zhuǎn)移矩陣（km2）Table3 The transfer matrix of land-use in research area from 2001 to 2008 and 2008 to 2016（km2）

2.2 土地利用類型與水質(zhì)的相關(guān)性分析

2.2.1 水質(zhì)變化趨勢(shì)

2001年紅楓湖水質(zhì)總體為Ⅲ類，水質(zhì)情況較好。2008年水質(zhì)污染嚴(yán)重，總體達(dá)Ⅳ類，局部呈現(xiàn)Ⅴ類、劣Ⅴ類水平，工業(yè)廢水、生活污水的排放以及農(nóng)業(yè)面源污染可能是導(dǎo)致水質(zhì)惡化的主要原因。2008年啟動(dòng)了污染治理工程，經(jīng)過(guò)治理，2016年水質(zhì)總體達(dá)Ⅱ類，局部呈現(xiàn)Ⅰ類水平（圖2），15年間流域水質(zhì)經(jīng)歷了由較好變差，由差轉(zhuǎn)為優(yōu)的過(guò)程。將土地利用類型（林地、耕地、草地、建設(shè)用地、未利用地）與水質(zhì)（TP、NH3-N、COD、DO）進(jìn)行 Spearman秩相關(guān)分析（表 4）。由相關(guān)性分析可以看出，建設(shè)用地、耕地與水質(zhì)營(yíng)養(yǎng)指標(biāo)普遍呈正相關(guān)關(guān)系，有些達(dá)到顯著水平，林地呈明顯負(fù)相關(guān)關(guān)系，草地對(duì)水質(zhì)的影響較小，未利用地對(duì)水質(zhì)的影響較為波動(dòng)，正負(fù)相關(guān)性都有存在，影響方向不明確。

表4 2001、2008年和2016年土地利用類型與水質(zhì)的相關(guān)性分析Table4 Correlation analysis between land-use type and water quality parameters in 2001,2008 and 2016

2.2.2 建設(shè)用地、耕地與水質(zhì)的相關(guān)性

建設(shè)用地與水質(zhì)指標(biāo)的相關(guān)性表現(xiàn)顯著，與TP、NH3-N、COD均在一般正相關(guān)和極顯著正相關(guān)之間，相關(guān)系數(shù)分別可達(dá)0.818、0.883、0.571，與DO均在一般負(fù)相關(guān)和顯著負(fù)相關(guān)之間，相關(guān)系數(shù)最高達(dá)-0.865。以不同尺度比較建設(shè)用地與TP、NH3-N的相關(guān)系數(shù)。當(dāng)尺度為1、1.5 km和2 km時(shí)，2001、2008年和2016年相關(guān)系數(shù)大小的順序總體為2016年＜2008年＜2001年，尺度為2.5 km時(shí)，2008年＜2016年＜2001年，且2001年不同尺度上建設(shè)用地的相關(guān)系數(shù)較其他土地利用類型相比最高，說(shuō)明建設(shè)用地與水質(zhì)營(yíng)養(yǎng)指標(biāo)在2001年的相關(guān)性最為顯著，表明建設(shè)用地所代表的工業(yè)和生活污染源對(duì)水中氮、磷元素的引入有重要貢獻(xiàn)作用。耕地與TP、NH3-N、COD均在一般正相關(guān)和極顯著正相關(guān)之間，相關(guān)系數(shù)分別可達(dá)0.757、0.857、0.893，與DO主要呈顯著負(fù)相關(guān)，相關(guān)系數(shù)最高達(dá)-0.847。以不同尺度比較耕地與TP、NH3-N的相關(guān)系數(shù)。當(dāng)尺度為1.5、2 km和2.5 km時(shí)，2001、2008年和2016年相關(guān)系數(shù)大小的順序總體為2001年＜2008年＜2016年，尺度為1 km時(shí)，2008年＜2001年＜2016年，且2016年不同尺度上耕地的相關(guān)系數(shù)較其他土地利用類型相比最高。結(jié)合緩沖區(qū)面積統(tǒng)計(jì)發(fā)現(xiàn)，耕地面積比例越高，對(duì)水質(zhì)的影響越顯著。說(shuō)明在2016年耕地所代表的農(nóng)業(yè)面源污染成為紅楓湖保護(hù)區(qū)內(nèi)的主要污染源，這一結(jié)果也與紅楓湖開(kāi)展五大治理工程后，點(diǎn)源污染得到有效控制相吻合。2008年不同尺度上土地利用類型與水質(zhì)的相關(guān)性規(guī)律不明顯，紅楓湖的水質(zhì)較差可能是生活污染源、工業(yè)污染源和農(nóng)業(yè)面源共同作用的結(jié)果。

2.2.3 林地、草地、未利用地與水質(zhì)的相關(guān)性

緩沖區(qū)域?yàn)? km時(shí)，林地與水質(zhì)的相關(guān)性并不明顯（表4）。隨著緩沖區(qū)域的擴(kuò)大，林地面積比例增加，與水中TP、NH3-N、COD濃度呈明顯負(fù)相關(guān)關(guān)系，具體表現(xiàn)在2008年和2016年。在2 km尺度的緩沖區(qū)域上，其相關(guān)系數(shù)可達(dá)-0.643、-0.750、-0.571。這是由于一方面林地屬于透水性下墊面，對(duì)氮磷元素滯留、吸收可降低面源污染對(duì)受納水體的影響[22]；另一方面，林地面積比例增加導(dǎo)致耕地、建設(shè)用地面積比例減少，使污染物輸出也相應(yīng)減少，這也與相關(guān)研究結(jié)果[23-26]一致。草地與各水質(zhì)營(yíng)養(yǎng)指標(biāo)的相關(guān)系數(shù)較小，且不同尺度上相關(guān)系數(shù)變化未見(jiàn)明顯規(guī)律，說(shuō)明草地對(duì)流域水質(zhì)影響較小。未利用地對(duì)水質(zhì)指標(biāo)的相關(guān)性較為波動(dòng)，例如與TP和NH3-N，兩者既有一般正相關(guān)，也存在一般負(fù)相關(guān)，沒(méi)有明確的變化規(guī)律，相關(guān)性指向不明確。未利用地屬于透水性下墊面，對(duì)徑流的截留作用較強(qiáng)。廖建榮[13]指出，DO與未利用地呈現(xiàn)正相關(guān)關(guān)系。結(jié)合表4發(fā)現(xiàn)，2001年和2008年的不同緩沖尺度上，未利用地與DO呈現(xiàn)負(fù)相關(guān)關(guān)系，2016年的不同緩沖尺度上，未利用地與DO呈現(xiàn)正相關(guān)關(guān)系。這可能是因?yàn)樵缒觊g紅楓湖流域周邊正處于開(kāi)發(fā)建設(shè)時(shí)期，未利用地較多，長(zhǎng)期堆放建筑垃圾，因此雨源型面源污染更為突出，導(dǎo)致了未利用地與DO相關(guān)性波動(dòng)的結(jié)果。

2.2.4 源匯效應(yīng)分析

為比較土地利用類型對(duì)水質(zhì)營(yíng)養(yǎng)指標(biāo)的影響，根據(jù)表4相關(guān)系數(shù)的正負(fù)分別將土地利用劃為“源”、“匯”兩大類。結(jié)果顯示建設(shè)用地、耕地與水體中TP、NH3-N、COD普遍呈正相關(guān)關(guān)系，有些達(dá)到顯著水平，且以耕地的相關(guān)系數(shù)為最高，與DO呈現(xiàn)明顯負(fù)相關(guān)關(guān)系。2008年和2016年，林地與TP、NH3-N、COD呈現(xiàn)明顯負(fù)相關(guān)關(guān)系，但未達(dá)到顯著水平。因此，耕地和建設(shè)用地對(duì)氮磷等污染負(fù)荷起“源”的作用，且以耕地為主，林地對(duì)氮磷等產(chǎn)出則起到“匯”的效應(yīng)。

3 結(jié)論

（1）紅楓湖流域保護(hù)區(qū)內(nèi)土地利用類型主要以林地、耕地為主，占到了總面積的60%以上。從2001年到2016年，林地和耕地面積呈現(xiàn)下降趨勢(shì)，其中林地減少了10.68km2，減幅達(dá)7.30%，耕地減少了8.61 km2，減幅達(dá)5.30%。建設(shè)用地和草地面積持續(xù)增加，其中建設(shè)用地增幅最大，達(dá)36.11%，面積增加了11.41 km2，草地增加了9.65 km2，增幅達(dá)14.62%，兩者面積的增加主要來(lái)源于耕地、林地和未利用地的流入。

（2）2001年水質(zhì)總體為Ⅲ類，建設(shè)用地與TP、NH3-N、COD呈現(xiàn)正相關(guān)關(guān)系，有些達(dá)到顯著水平，表明在沒(méi)有進(jìn)行環(huán)境治理之前，污染主要來(lái)自于工業(yè)和生活污染源；2008年不同尺度上土地利用類型與水質(zhì)的相關(guān)性規(guī)律不明顯，水質(zhì)總體為Ⅳ類，紅楓湖水質(zhì)較差可能是生活污染源、工業(yè)污染源和農(nóng)業(yè)面源共同作用的結(jié)果；2016年耕地與 TP、NH3-N、COD和DO的相關(guān)系數(shù)分別達(dá)0.757*、0.750、0.786*和-0.847*，說(shuō)明在工業(yè)污染源和生活污染源得到治理后，農(nóng)業(yè)面源污染成為紅楓湖流域的主要污染源，但其排放系數(shù)較小，水質(zhì)總體為Ⅱ類，沒(méi)有對(duì)水質(zhì)產(chǎn)生較大影響。

（3）紅楓湖流域保護(hù)區(qū)內(nèi)建設(shè)用地、耕地與TP、NH3-N、COD普遍呈現(xiàn)正相關(guān)關(guān)系，與DO呈現(xiàn)明顯負(fù)相關(guān)關(guān)系，表現(xiàn)為污染“源”的作用，林地與TP、NH3-N、COD呈現(xiàn)明顯負(fù)相關(guān)關(guān)系，表現(xiàn)出“匯”的效應(yīng)。

[1]夏 叡,李云梅,王 橋,等.基于遙感的無(wú)錫市土地利用與過(guò)境水質(zhì)響應(yīng)關(guān)系的研究[J].地理科學(xué),2010,30（1）：129-133.

XIA Rui,LI Yun-mei,WANG Qiao,et al.Response relationship between land-use and transit water quality in Wuxi City based on remote sensing[J].Scientia Geographica Sinica,2010,30（1）：129-133.

[2]孟曉云,于興修,泮雪芹.云蒙湖流域土地利用變化對(duì)非點(diǎn)源氮污染負(fù)荷的影響[J].環(huán)境科學(xué),2012,33（6）：1789-1793.

MENG Xiao-yun,YU Xing-xiu,PAN Xue-qin.Impact of the land-use change on the non-point source nitrogen load in Yunmeng Lake watershed[J].Environmental Sciences,2012,33（6）：1789-1793.

[3]楊 峰,王鵬舉,楊珊珊,等.城郊農(nóng)業(yè)區(qū)小流域土地利用結(jié)構(gòu)對(duì)氮素輸出的影響[J].環(huán)境科學(xué),2012,33（8）：2652-2658.

YANG Feng,WANG Peng-ju,YANG Shan-shan,et al.Influence of land use structure on nitrogen output in the watershed of suburban agriculture regions[J].Environmental Sciences,2012,33（8）：2652-2658.

[4]李兆富,楊桂山,李恒鵬.西苕溪流域土地利用對(duì)氮素輸出影響研究[J].環(huán)境科學(xué),2006,27（3）：498-502.

LI Zhao-fu,YANG Gui-shan,LI Heng-peng.Effects of land use on nitrogen export in Xitiaoxi watershed[J].Environmental Sciences,2006,27（3）：498-502.

[5]于興修,楊桂山.中國(guó)土地利用/覆被變化研究的現(xiàn)狀與問(wèn)題[J].地理科學(xué)進(jìn)展,2002,21（1）：51-57.

YU Xing-xiu,YANG Gui-shan.The advances and problems of land use and land cover change research in China[J].Progress in Geography,2002,21（1）：51-57.

[6]許俐俐,胡寶清,嚴(yán)志強(qiáng),等.喀斯特土地利用變化及其區(qū)域生態(tài)環(huán)境效應(yīng)：以廣西都安瑤族自治縣為例[J].地域研究與開(kāi)發(fā),2005,24（6）：66-71.

XU Li-li,HU Bao-qing,YAN Zhi-qiang,et al.Karst land use change and their regional eco-environmental effect：A case study of Du′an Yao Autonomous County of Guangxi[J].Areal Research and Development,2005,24（6）：66-71.

[7]彭 建,王仰麟,張 源,等.滇西北生態(tài)脆弱區(qū)土地利用變化及其生態(tài)效應(yīng)：以云南省永勝縣為例[J].地理學(xué)報(bào),2004,59（4）：629-638.

PENG Jian,WANG Yang-lin,ZHANG Yuan,et al.Land use change and its ecological effect in the ecotone of northwest of Yunnan Province,China：A case study of Yongsheng County[J].Acta Geographica Sinica,2004,59（4）：629-638.

[8]蔡 宏,何政偉,安艷玲,等.基于遙感和GIS的赤水河水質(zhì)對(duì)流域土地利用的響應(yīng)研究[J].長(zhǎng)江流域資源與環(huán)境,2015,24（2）：287-291.

CAI Hong,HE Zheng-wei,AN Yan-ling,et al.Response relationship between land use and water quality in Chishui River basin on RS and GIS[J].Resources and Environment in the Yangtze Basin,2015,24（2）：287-291.

[9]Nash M S,Heggem D T,Ebert D,et al.Multi-scale landscape factors influencing stream water quality in the state of oregon[J].Environmental Monitoring and Assessment,2009,156（1/2/3/4）：343-360.

[10]曹芳芳,李 雪,王 東,等.新安江流域土地利用結(jié)構(gòu)對(duì)水質(zhì)的影響[J].環(huán)境科學(xué),2013,34（7）：2582-2587.

CAO Fang-fang,LI Xue,WANG Dong,et al.Effects of land use structure on water quality in Xin′anjiang River[J].Environmental Sciences,2013,34（7）：2582-2587.

[11]郭青海,馬克明,楊 柳.城市非點(diǎn)源污染的主要來(lái)源及分類控制對(duì)策[J].環(huán)境科學(xué),2006,27（11）：2170-2175.

GUO Qing-hai,MA Ke-ming,YANG Liu.Main sources of urban nonpoint source pollution and control measures for classified catchments[J].Environmental Sciences,2006,27（11）：2170-2175.

[12]孫金華,曹曉峰,黃 藝.滇池流域土地利用對(duì)入湖河流水質(zhì)的影響[J].中國(guó)環(huán)境科學(xué),2011,31（12）：2052-2057.

SUN Jin-hua,CAO Xiao-feng,HUANG Yi.Effect of land use on inflow rivers water quality in lake Dianchi watershed[J].China Environmental Science,2011,31（12）：2052-2057.

[13]廖建榮.山地城市土地利用與河流水質(zhì)響應(yīng)關(guān)系的研究：以重慶市龍景湖流域?yàn)槔齕D].重慶：重慶大學(xué),2016.

LIAO Jian-rong.Research on the response relationship between river water quality and land-using in mountainous city：Taking L-J Lake watershed for example[D].Chongqing：Chongqing University,2016.

[14]鮑全盛,王華東.我國(guó)水環(huán)境非點(diǎn)源污染研究與展望[J].地理科學(xué),1996,16（1）：66-72.

BAO Quan-sheng,WANG Hua-dong.The research and prospect on non-point source pollution on water environment in China[J].Scientia Geographica Sinica,1996,16（1）：66-72.

[15]肖致強(qiáng),安艷玲.2002—2009年紅楓湖水污染趨勢(shì)分析[J].環(huán)境科學(xué)導(dǎo)刊,2012,31（3）：29-34.

XIAO Zhi-qiang,AN Yan-ling.An analysis of water environment changes of lake Hongfeng（2002—2009）[J].Environmental Science Survey,2012,31（3）：29-34.

[16]楊通銓,劉鴻雁,喻陽(yáng)華.紅楓湖水質(zhì)變化趨勢(shì)及原因分析[J].長(zhǎng)江流域資源與環(huán)境,2014,23（Z1）：96-101.

YANG Tong-quan,LIU Hong-yan,YU Yang-hua.Variation trend of water quality and its causing effect of the Hongfeng Lake[J].Resources and Environment in the Yangtze Basin,2014,23（Z1）：96-101.

[17]劉紀(jì)遠(yuǎn).中國(guó)資源環(huán)境遙感宏觀調(diào)查與動(dòng)態(tài)研究[M].北京：中國(guó)科學(xué)技術(shù)出版社,1996.

LIU Ji-yuan.Macro-scale survey and dynamic study of natural resources and environment of China by remote sensing[M].Beijing：China Science and Technology Press,1996.

[18]劉紀(jì)遠(yuǎn),匡文慧,張?jiān)鱿?等.20世紀(jì)80年代末以來(lái)中國(guó)土地利用變化的基本特征與空間格局[J].地理學(xué)報(bào),2014,69（1）：3-14.

LIU Ji-yuan,KUANG Wen-hui,ZHANG Zeng-xiang,et al.Spatiotemporal characteristics,patterns and causes of land use changes in China since the late 1980s[J].Acta Geographica Sinica,2014,69（1）：3-14.

[19]Liu J Y,Zhang Z X,Xu X L,et al.Spatial patterns and driving forces of land use change in China during the early 21st Century[J].Journal of Geographical Sciences,2010,20（4）：483-494.

[20]陳 江,馬松梅,劉 琳,等.博樂(lè)墾區(qū)土地利用/覆被變化及預(yù)測(cè)研究[J].水土保持研究,2015,22（3）：44-49.

CHEN Jiang,MA Song-mei,LIU Lin,et al.Study on land use/cover change and the prediction of Bole reclamation area[J].Research of Soil and Water Conservation,2015,22（3）：44-49.

[21]陳利頂,傅伯杰,趙文武.“源”“匯”景觀理論及其生態(tài)學(xué)意義[J].生態(tài)學(xué)報(bào),2006,26（5）：1444-1449.

CHEN Li-ding,FU Bo-jie,ZHAO Wen-wu.Source-sink landscape theory and its ecological significance[J].Acta Ecologica Sinica,2006,26（5）：1444-1449.

[22]楊 柳,馬克明,郭青海,等.城市化對(duì)水體非點(diǎn)源污染的影響[J].環(huán)境科學(xué),2004,25（6）：32-39.

YANG Liu,MA Ke-ming,GUO Qing-hai,et al.Impacts of the urbanization on waters non-point source pollution[J].Environmental Sciences,2004,25（6）：32-39.

[23]Sliva L,Williams D D.Buffer zone versus whole catchment approaches to studying land use impact on river water quality[J].Water Research,2001,35（14）：3462-3472.

[24]帥 紅,夏北成.廣佛區(qū)域土地利用結(jié)構(gòu)對(duì)非點(diǎn)源污染的影響[J].熱帶地理,2006,26（3）：230-233.

SHUAI Hong,XIA Bei-cheng.Effects of land use structure on nonpoint source pollution in the areas of Guangzhou-Foshan[J].Tropical Geography,2006,26（3）：230-233.

[25]Moreno J L,Navarro C,Delas H J.Abiotic ecotypes in south-central Spanish rivers：Reference conditions and pollution[J].Environmental Pollution,2006,143：388-396.

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

FU Bo-jie,MA Ke-ming,ZHOU Hua-feng,et al.Impacts of land use structure on soil nutrient distribution in loess hilly region[J].Chinese Science Bulletin,1998,43（22）：2444-2447.