項 頌,龐 燕,*,竇嘉順,呂興菊,薛力強(qiáng),儲昭升

1 中國環(huán)境科學(xué)研究院湖泊生態(tài)環(huán)境創(chuàng)新基地,北京 100012 2 大理洱海湖泊研究中心, 大理 671000

人類的不合理活動導(dǎo)致土地利用發(fā)生改變,土地利用的改變又反過來影響著人類生活的諸多領(lǐng)域,而土地利用變化對地表水環(huán)境的影響是業(yè)界關(guān)注的重點之一[1- 4]。入湖河流既是外源匯入湖泊的重要通道,也是湖泊水資源的重要補(bǔ)給源,其水質(zhì)水量變化直接影響湖泊水環(huán)境質(zhì)量[5-6]。因此,進(jìn)行土地利用與入湖河流關(guān)系研究對流域水環(huán)境保護(hù)及土地開發(fā)利用意義重大。

土地利用作為人類活動的產(chǎn)物,通過改變流域的水文循環(huán)、水土流失、養(yǎng)分遷移轉(zhuǎn)化等生態(tài)過程,進(jìn)而對河流生態(tài)環(huán)境產(chǎn)生影響[7]。入湖河流作為生態(tài)景觀的重要組分,其與土地利用的關(guān)系受空間尺度和時間變化影響[8-9]。針對兩者關(guān)系的時空差異性研究目前已開展了很多,例如Sliva和Williams[10]開展了全流域和100m緩沖區(qū)尺度下加拿大Ontario省南部3個流域水質(zhì)與土地利用關(guān)系的對比研究,Tudesque等[11]對5種不同空間尺度下法國Adour-Garonne流域土地利用和水質(zhì)、藻類的關(guān)系進(jìn)行了研究,Ding等[12]開展了東江流域多空間尺度下土地利用格局與水質(zhì)的關(guān)系研究,他們的研究均表明相對于河岸緩沖帶尺度,土地利用類型對水質(zhì)的影響在流域尺度更顯著；而Shen、Zhou等[13-14]的研究則得出相反結(jié)論,Shen等在北運(yùn)河流域開展全流域尺度和5種緩沖區(qū)尺度下土地利用與水質(zhì)關(guān)系的研究,結(jié)果表明較全流域尺度,100m緩沖區(qū)尺度下土地利用格局對河流水質(zhì)的影響更顯著,Zhou等對1990—2006年東江流域土地利用與水質(zhì)的關(guān)系進(jìn)行研究,選取全流域、子流域、緩沖區(qū)3種尺度開展兩者關(guān)系的對比研究,結(jié)果表明多數(shù)水質(zhì)指標(biāo)在3種空間尺度下均受土地利用格局的影響,其中子流域尺度下兩者關(guān)系較全流域、緩沖區(qū)尺度強(qiáng)；此外Yu、lvarez-Cabria、Prathumratana等[15-17]對兩者關(guān)系受季節(jié)性如降雨、氣候等因素影響而產(chǎn)生的變化也進(jìn)行了研究,研究表明土地利用對河流水質(zhì)的影響與季節(jié)性變化關(guān)系密切,不同季節(jié)下降雨量、溫度和農(nóng)業(yè)活動的差異對兩者關(guān)系影響較大。研究手段包括傳統(tǒng)分析方法如相關(guān)性分析[18-19]、多元回歸分析[14,20],以及近年發(fā)展起來的地理信息技術(shù)[21]及模型[22-23]。雖然目前涉及這方面的研究較多且已走向深入,但受各流域自身特性、人類干擾強(qiáng)度和數(shù)據(jù)精度的影響,土地利用與水質(zhì)關(guān)系的研究在某些方面仍未有新突破,如土地利用與水質(zhì)關(guān)系的最適或最強(qiáng)空間尺度問題一直未形成統(tǒng)一定論[10,12-13,24]。

洱海作為云南省第二大高原淡水湖泊,是大理州乃至云南省社會經(jīng)濟(jì)可持續(xù)發(fā)展的基礎(chǔ),但隨流域社會經(jīng)濟(jì)發(fā)展和區(qū)域人口增加,流域土地開發(fā)利用加劇,洱海水質(zhì)日益下降。洱海西區(qū)作為洱海入湖河流的核心區(qū),為洱海提供重要清潔水源,提高洱海水質(zhì),但由于西區(qū)地理環(huán)境優(yōu)越、水資源豐富,成為流域人類活動的中心區(qū)域,區(qū)內(nèi)高密度的人類活動和高強(qiáng)度的土地開發(fā)問題十分突出,以非點源污染為主的環(huán)境污染問題日趨加重,影響入湖河流水質(zhì),威脅洱海水環(huán)境,因此如何平衡土地開發(fā)與水環(huán)境保護(hù)是該區(qū)目前正面臨的重要難題,也是改善區(qū)域河流水質(zhì),保障清水入洱海,促進(jìn)流域經(jīng)濟(jì)發(fā)展的關(guān)鍵問題。該文在此背景下,結(jié)合前期的研究[25-26],基于研究區(qū)土地利用特征,采用數(shù)理統(tǒng)計和空間分析方法,開展入湖河流水質(zhì)隨土地利用空間尺度及季節(jié)性效應(yīng)的變化研究,以期為當(dāng)?shù)厮h(huán)境保護(hù)和土地開發(fā)利用提供一定參考。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于云南省大理白族自治州大理市洱海流域西部,99°58′—100°14′E、25°35′—25°56′N之間(圖1)。區(qū)域面積336.9km2,占洱海流域面積的13.1%；氣候特征為典型的低緯亞熱帶高原季風(fēng)氣候,多年年均溫度為15.1℃；干濕季分明,5—10月為雨季,11月—次年4月為旱季,多年年均降雨量為1048mm,雨季降水量占全年降水量的88.7%,年內(nèi)降水高峰月出現(xiàn)在7—8月。

研究區(qū)水系發(fā)達(dá),水資源豐富,18條垂直于洱海的山溪型河流構(gòu)成“蒼山十八溪”水系,是洱海清水補(bǔ)給的重要水源地,年平均地表徑流量2.339億m3,占洱海水資源量的24.7%。蒼山十八溪發(fā)源于蒼山,流經(jīng)人類活動密集的大理壩子,灌溉壩區(qū)萬頃良田,提供大理市居民生活飲用水,最終進(jìn)入洱海,為洱海提供清水補(bǔ)給。作為洱海的主要水源和核心的入湖河道,蒼山十八溪水質(zhì)水量變化直接影響洱海水環(huán)境質(zhì)量[27]。因陽南溪、葶溟溪、雙鴛溪、霞移溪4條河流一年中多數(shù)時間處于斷流狀態(tài),該次研究不包括這4條河流。

1.2 數(shù)據(jù)來源

遙感影像數(shù)據(jù)為2012年資源一號02c衛(wèi)星拍攝的2.36m高精度衛(wèi)片,受區(qū)域人類活動及經(jīng)濟(jì)社會發(fā)展等因素的影響,經(jīng)調(diào)研發(fā)現(xiàn)土地利用變化相對緩慢,該文選取的2012年遙感影像數(shù)據(jù)基本能反映2012年前后一段時間的土地利用結(jié)構(gòu)特征。DEM為中國科學(xué)院計算機(jī)網(wǎng)絡(luò)信息中心ASTER GDEM 30m分辨率高程數(shù)據(jù)。

1.3 分析方法

1.3.1 空間分析

基于遙感影像,參照最新土地分類標(biāo)準(zhǔn)土地利用現(xiàn)狀分類[29],采用非監(jiān)督分類法進(jìn)行土地利用分類,并以洱海流域第二次全國土地調(diào)查結(jié)果為參考修正,分為建設(shè)用地、農(nóng)業(yè)用地、植被、水體、其他用地5類(圖2)。其中植被包括林地和牧草地,雖然林地和牧草地對水質(zhì)的凈化作用不同,但鑒于研究區(qū)牧草地面積極小,林地面積大,將兩者歸為一類；其他用地包括特殊用地、自然保留地和風(fēng)景名勝設(shè)施用地。

圖1 研究區(qū)位置Fig.1 Location of study area

圖2 研究區(qū)土地利用空間分布Fig.2 The spatial distribution of land use in study area

在Arcgis 10.1(ESRI)平臺上,基于流域DEM利用水文分析工具劃分小流域,提取河網(wǎng),并在此基礎(chǔ)上劃分河岸帶。目前河岸帶劃分類型主要有矩形和圓形兩種,對于河流的研究多以矩形為主[30]；研究區(qū)河流為山溪型河流,流域較窄小,結(jié)合區(qū)域土地利用空間分布特征,采用緩沖區(qū)法沿河流分別向兩岸劃分30、60、90m等3種寬度的矩形河岸帶緩沖區(qū)。

1.3.2 統(tǒng)計分析

基于SPSS 19軟件進(jìn)行河流水質(zhì)特征分析。分析之前采用SPSS 19進(jìn)行數(shù)據(jù)異常值剔除和篩選,并進(jìn)行數(shù)據(jù)正態(tài)性檢驗,發(fā)現(xiàn)所有數(shù)據(jù)的雙尾漸進(jìn)概率P值均大于給定的顯著性水平0.05,數(shù)據(jù)服從正態(tài)分布。在此基礎(chǔ)上采用描述性統(tǒng)計分析方法對14條入湖河流水質(zhì)數(shù)據(jù)進(jìn)行統(tǒng)計分析,確定水質(zhì)變化特征。

在Canoco 4.5軟件平臺上進(jìn)行土地利用和入湖河流水質(zhì)的相關(guān)性分析。采用趨勢對應(yīng)分析(DCA)確定排序模型,發(fā)現(xiàn)排序軸最大梯度值小于3,則選用常用的線性排序模型冗余分析(RDA)進(jìn)行兩者相關(guān)性研究[10,13]。

基于SPSS 19軟件,在0.05的顯著性水平下,采用逐步多元回歸分析方法建立入湖河流水質(zhì)對土地利用的響應(yīng)關(guān)系,回歸的進(jìn)入概率和刪除概率分別為0.05和0.1[31]。

2 結(jié)果與分析

2.1 研究區(qū)入湖河流水質(zhì)與土地利用特征

結(jié)合GIS工具,對不同空間尺度下14條入湖河流所在小流域土地利用特征進(jìn)行空間統(tǒng)計分析(圖3),結(jié)果顯示：研究區(qū)土地利用類型以植被、農(nóng)業(yè)用地和建設(shè)用地為主,小流域尺度下植被為優(yōu)勢地類,河岸帶尺度下植被和農(nóng)業(yè)用地占比最高,建設(shè)用地次之。14個小流域各類土地利用面積占比在不同空間尺度下存在差異,就其平均值而言,植被：小流域>河岸帶90m緩沖區(qū)>河岸帶60m緩沖區(qū)>河岸帶30m緩沖區(qū),河岸帶緩沖區(qū)尺度下差異不大,整體而言植被面積占比在河岸帶緩沖區(qū)尺度小于小流域尺度；建設(shè)用地和農(nóng)業(yè)用地：河岸帶60m緩沖區(qū) >河岸帶30m和90m緩沖區(qū)>小流域,兩種地類面積占比均在河岸帶緩沖區(qū)尺度大于流域尺度；水體：隨河岸帶緩沖區(qū)尺度增大其面積占比減小,整體而言水體面積占比在河岸帶緩沖區(qū)尺度大于小流域尺度；其他用地：河岸帶90m緩沖區(qū)>河岸帶30m緩沖區(qū)>河岸帶60m緩沖區(qū) >小流域,整體而言其他用地面積占比在河岸帶緩沖區(qū)尺度大于流域尺度。

圖3 不同空間尺度下土地利用面積比例Fig.3 Land use area proportion under multiple scales

2.2 土地利用與入湖河流水質(zhì)的關(guān)系

通過冗余分析進(jìn)行不同空間尺度下土地利用面積百分比與入湖河流水質(zhì)的關(guān)系研究,檢驗不同空間尺度下土地利用對于水質(zhì)變異程度的解釋能力,從而篩選對水質(zhì)影響最大的空間尺度(圖4和表2)。結(jié)果顯示：土地利用類型面積百分比對河流水質(zhì)影響程度的空間尺度排序為小流域>河岸帶30m緩沖區(qū)>河岸帶60m緩沖區(qū)>河岸帶90m緩沖區(qū),小流域尺度下土地利用面積百分比對河流水質(zhì)總解釋率為66%,是對河流水質(zhì)影響最大的空間尺度；所選取的3種河岸帶緩沖區(qū)尺度下,河岸帶30m緩沖區(qū)尺度下土地利用面積百分比對河流水質(zhì)影響最顯著,解釋率為47.4%。

圖4 土地利用與入湖河流水質(zhì)關(guān)系的冗余分析Fig.4 Redundancy analysis between land use and inflow river water qualitycon：建設(shè)用地,construction land；agr：農(nóng)業(yè)用地,agricultural land；veg：植被,vegetation land；wat：水體,water；oth：其他用地,other land；COD：化學(xué)需氧量,Chemical Oxygen Demand；氨氮,Ammonia Nitrogen；TN：總氮,Total Nitrogen；TP：總磷,Total Phosphorus

空間尺度Spatialscales參數(shù)Parameters第一軸Axis1第二軸Axis2第三軸Axis3第四軸Axis4總解釋率/%Totalexplainedvariance小流域Subcatchment特征值0.5050.1360.019066.00土地利用和入湖河流水質(zhì)相關(guān)性的累積百分率/%76.5097.10100100河岸帶30m緩沖區(qū)特征值0.3920.0720.009047.4030mriparianbufferzone土地利用和入湖河流水質(zhì)相關(guān)性的累積百分率/%82.8098.1099.90100河岸帶60m緩沖區(qū)特征值0.3850.0730.012046.9060mriparianbufferzone土地利用和入湖河流水質(zhì)相關(guān)性的累積百分率/%8297.5099.90100河岸帶90m緩沖區(qū)特征值0.380.0530.0030.00143.7090mriparianbufferzone土地利用和入湖河流水質(zhì)相關(guān)性的累積百分率/%86.9099.0099.70100.00

圖5 不同空間尺度下各土地利用類型解釋率Fig.5 Explained ratio of different land use under multiple scales

各種土地利用類型在不同空間尺度下對河流水質(zhì)的解釋率如圖5所示,結(jié)果顯示：選取的小流域、河岸帶30m緩沖區(qū)、河岸帶60m緩沖區(qū)、河岸帶90m緩沖區(qū)4種空間尺度下,整體而言對河流水質(zhì)解釋率最高的地類為建設(shè)用地和植被,偏蒙特卡羅置換檢驗的P均小于0.05,為顯著影響地類；但小流域尺度下,除建設(shè)用地(P=40.7%)和植被(P=24.8%)對河流水質(zhì)的解釋率較高,其他用地類型對河流水質(zhì)也有一定程度的影響。

2.3 入湖河流水質(zhì)對土地利用的響應(yīng)關(guān)系

表3 入湖河流水質(zhì)對土地利用類型的響應(yīng)關(guān)系

2.4 土地利用與入湖河流水質(zhì)的季節(jié)性關(guān)聯(lián)

圖6 建設(shè)用地與水質(zhì)的季節(jié)性關(guān)聯(lián)分析Fig.6 Linear regression between percentage of construction land area and water quality under different seasons

圖7 植被與水質(zhì)的季節(jié)性關(guān)聯(lián)分析Fig.7 Linear regression between percentage of vegetation land use area and water quality under different seasons

3 討論

3.1 土地利用與入湖河流水質(zhì)關(guān)系的尺度效應(yīng)

土地利用與入湖河流水質(zhì)存在尺度效應(yīng)已成定論,但關(guān)于兩者關(guān)系的最佳空間尺度并未形成統(tǒng)一結(jié)論[12,21,24]。該文采用冗余分析進(jìn)行土地利用與河流水質(zhì)的空間效應(yīng)研究,發(fā)現(xiàn)兩者的關(guān)系隨空間尺度變化不同,小流域尺度下土地利用對水質(zhì)的解釋率最高(達(dá)66%),是對入湖河流水質(zhì)影響最強(qiáng)的空間尺度[10-12,14],而Shen、Zhou等[13-14]的研究則發(fā)現(xiàn)相對于全流域尺度,河岸帶尤其是100m寬度的河岸帶尺度下土地利用與入湖河流水質(zhì)的關(guān)系更顯著。該文出現(xiàn)與此不同的研究結(jié)論是基于研究區(qū)自身的土地利用特征,研究區(qū)小流域尺度下土地利用優(yōu)勢地類為植被,面積占比超過50%,遠(yuǎn)高于其他地類,優(yōu)勢地類突出且單一,其與水質(zhì)的顯著關(guān)系使得總解釋率較高；而河岸帶尺度下,植被、農(nóng)業(yè)用地和建設(shè)用地共同構(gòu)成優(yōu)勢地類,優(yōu)勢地類組分復(fù)雜且各地類對水質(zhì)的影響效果各異,因此疊加作用下其對水質(zhì)的解釋率較小流域尺度低[11,15,22]。類似的結(jié)論出現(xiàn)在Sliva、Tudesque、Ding等[10-12]的研究中,他們也發(fā)現(xiàn)土地利用與入湖河流水質(zhì)的最佳影響尺度為流域尺度。

該次研究中小流域尺度之所以對河流水質(zhì)影響最大,與其高植被覆蓋率密切相關(guān)；河岸帶是阻隔污染物遷移入河的最后屏障,同樣對保護(hù)河流水質(zhì)意義重大[19,32]。目前洱海流域西部部分山溪河流受人類活動影響,生態(tài)功能出現(xiàn)退化或逐漸喪失,因此除保證小流域整體植被覆蓋率不降低外,有必要提高河岸帶周邊的植被覆蓋率、理性約束人類活動。

3.2 入湖河流水質(zhì)對土地利用的響應(yīng)關(guān)系

3.3 土地利用與入湖河流水質(zhì)關(guān)系的季節(jié)效應(yīng)

土地利用與入湖河流水質(zhì)的關(guān)系是非點源研究的重要內(nèi)容之一,與降雨徑流密切相關(guān)[7,15-16]。已有的關(guān)于兩者關(guān)系的季節(jié)性效應(yīng)研究結(jié)論各異,如孫芹芹等[35]對九龍江流域城鎮(zhèn)建設(shè)用地與河流水質(zhì)關(guān)系的季節(jié)性關(guān)聯(lián)研究發(fā)現(xiàn),隨城鎮(zhèn)建設(shè)用地面積增加,河流水質(zhì)在豐水期和枯水期惡化均較快；Li等[19]對漢江流域土地利用驅(qū)動下水質(zhì)的變化規(guī)律研究發(fā)現(xiàn),豐水期土地利用與入湖河流水質(zhì)關(guān)系在緩沖區(qū)尺度下較全流域尺度弱,枯水期則相反；Pratt和Chang[9]的研究發(fā)現(xiàn)河流水質(zhì)雨季與城鎮(zhèn)用地關(guān)系顯著,旱季與高程、地形等表征地理特征的因子關(guān)系強(qiáng)。而該次數(shù)據(jù)量有限,研究發(fā)現(xiàn)土地利用與入湖河流水質(zhì)關(guān)系的存在一定季節(jié)性差異,其中雨季建設(shè)用地和植被與水質(zhì)的回歸關(guān)系較旱季更強(qiáng),這主要是由于土地利用通過改變地表水文循環(huán)、污染遷移轉(zhuǎn)化途徑進(jìn)而影響入湖河流水質(zhì),雨季其對水質(zhì)的作用機(jī)理受徑流沖刷影響更大,更多污染物能被雨季地表徑流沖刷進(jìn)入河流[7,19],因此在土地利用管控時尤其應(yīng)限制雨季建設(shè)用地面積擴(kuò)張,同時保障植被覆蓋率。此外為更好的揭示土地利用與入湖河流水質(zhì)關(guān)系的季節(jié)效應(yīng),后續(xù)研究中應(yīng)盡可能獲取更多的數(shù)據(jù)用于此部分內(nèi)容的深入研究。

4 結(jié)論

1)選取的小流域、河岸帶30m緩沖區(qū)、河岸帶60m緩沖區(qū)、河岸帶90m緩沖區(qū)4種空間尺度下對河流水質(zhì)影響最主要的地類為建設(shè)用地和植被,最強(qiáng)空間尺度為小流域尺度,河岸帶30m緩沖區(qū)次之。最強(qiáng)空間尺度下土地利用對水質(zhì)的解釋率最高為66%,建設(shè)用地和植被對水質(zhì)的解釋率分別為40.7%和24.8%。

[1] Meneses B M, Reis R, Vale M J, Saraiva R. Land use and land cover changes in Zêzere watershed (Portugal)—Water quality implications. Science of the Total Environment, 2015, 527- 528: 439- 447.

[2] Giri S, Qiu Z Y. Understanding the relationship of land uses and water quality in Twenty First Century: A review. Journal of Environmental Management, 2016, 173: 41- 48.

[3] Sajikumar N, Remya R S. Impact of land cover and land use change on runoff characteristics. Journal of Environmental Management, 2015, 161: 460- 468.

[4] Zhou P, Huang J L, Pontius R G Jr, Hong H S. New insight into the correlations between land use and water quality in a coastal watershed of China: does point source pollution weaken it? Science of the Total Environment, 2016, 543: 591- 600.

[5] Zhang S H, Li Y Q, Zhang T X, Peng Y. An integrated environmental decision support system for water pollution control based on TMDL-A case study in the Beiyun River watershed. Journal of Environmental Management, 2015, 156: 31- 40.

[6] Dong C Y, Zhang W G, Ma H L, Feng H, Lu H H, Dong Y, Yu L Z. A magnetic record of heavy metal pollution in the Yangtze River subaqueous delta. Science of the Total Environment, 2014, 476- 477: 368- 377.

[7] Putro B, Kjeldsen T R, Hutchins M G, Miller J. An empirical investigation of climate and land-use effects on water quantity and quality in two urbanising catchments in the southern United Kingdom. Science of the Total Environment, 2016, 548- 549: 164- 172.

[8] Yu S Y, Xu Z X, Wu W, Zuo D P. Effect of land use types on stream water quality under seasonal variation and topographic characteristics in the Wei River basin, China. Ecological Indicators, 2016, 60: 202- 212.

[9] Pratt B, Chang H. Effects of land cover, topography, and built structure on seasonal water quality at multiple spatial scales. Journal of Hazardous Materials, 2012, 209- 210: 48- 58.

[10] Sliva L, Williams D D. Buffer zone versus whole catchment approaches to studying land use impact on river water quality. Water Research, 2001, 35(14): 3462- 3472.

[11] Tudesque L, Tisseuil C, Lek S. Scale-dependent effects of land cover on water physico-chemistry and diatom-based metrics in a major river system, the Adour-Garonne basin (South Western France). Science of the Total Environment, 2014, 466- 467: 47- 55.

[12] Ding J, Jiang Y, Liu Q, Hou Z J, Liao J Y, Fu L, Peng Q Z. Influences of the land use pattern on water quality in low-order streams of the Dongjiang River basin, China: A multi-scale analysis. Science of the Total Environment, 2016, 551- 552: 205- 216.

[13] Shen Z Y, Hou X S, Li W, Anni G, Chen L, Gong Y W. Impact of landscape pattern at multiple spatial scales on water quality: A case study in a typical urbanised watershed in China. Ecological Indicators, 2015, 48: 417- 427.

[14] Zhou T, Wu J G, Peng S L. Assessing the effects of landscape pattern on river water quality at multiple scales: A case study of the Dongjiang River watershed, China. Ecological Indicators, 2012, 23: 166- 175.

[15] Yu X B, Hawley-Howard J, Pitt A L, Wang J J, Baldwin R F, Chow A T. Water quality of small seasonal wetlands in the Piedmont ecoregion, South Carolina, USA: Effects of land use and hydrological connectivity. Water Research, 2015, 73: 98- 108.

[17] Prathumratana L, Sthiannopkao S, Kim K W. The relationship of climatic and hydrological parameters to surface water quality in the lower Mekong River. Environment International, 2008, 34(6): 860- 866.

[18] Tran C P, Bode R W, Smith A J, Kleppel G S. Land-use proximity as a basis for assessing stream water quality in New York State (USA). Ecological Indicators, 2010, 10(3): 727- 733.

[19] Li S Y, Gu S, Tan X, Zhang Q F. Water quality in the upper Han River basin, China: The impacts of land use/land cover in riparian buffer zone. Journal of Hazardous Materials, 2009, 165(1/3): 317- 324.

[20] Seeboonruang U. A statistical assessment of the impact of land uses on surface water quality indexes. Journal of Environmental Management, 2012, 101: 134- 142.

[21] Smart R P, Soulsby C, Cresser M S, Wade A J, Townend J, Billett M F, Langan S. Riparian zone influence on stream water chemistry at different spatial scales: a GIS-based modelling approach, an example for the Dee, NE Scotland. Science of the Total Environment, 2001, 280(1/3): 173- 193.

[22] Chang H. Spatial analysis of water quality trends in the Han River basin, South Korea. Water Research, 2008, 42(13): 3285- 3304.

[23] Tu J, Xia Z G. Examining spatially varying relationships between land use and water quality using geographically weighted regression I: model design and evaluation. Science of the Total Environment, 2008, 407(1): 358- 378.

[24] Qian S S. Ecological threshold and environmental management: a note on statistical methods for detecting thresholds. Ecological Indicators, 2014, 38: 192- 197.

[25] 翟玥, 尚曉, 沈劍, 王欣澤. SWAT模型在洱海流域面源污染評價中的應(yīng)用. 環(huán)境科學(xué)研究, 2012, 25(6): 666- 671.

[26] 龐燕, 項頌, 儲昭升, 薛力強(qiáng), 葉碧碧. 洱海流域農(nóng)業(yè)用地與入湖河流水質(zhì)的關(guān)系研究. 環(huán)境科學(xué), 2015, 36(11): 4005- 4012.

[27] 翟羽佳, 劉春學(xué). 蒼山十八溪流域水資源生態(tài)系統(tǒng)服務(wù)功能價值評估. 中國農(nóng)村水利水電, 2015, (5): 77- 80.

[28] 國家環(huán)境保護(hù)總局, 廢水監(jiān)測分析方法編委會. 水和廢水監(jiān)測分析方法(第四版). 北京: 中國環(huán)境科學(xué)出版社, 2002: 210- 281.

[29] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局, 中國國家標(biāo)準(zhǔn)化管理委員會. GB/T 21010- 2007 土地利用現(xiàn)狀分類. 北京: 中國標(biāo)準(zhǔn)出版社, 2007.

[30] 韓路, 王海珍, 于軍. 河岸帶生態(tài)學(xué)研究進(jìn)展與展望. 生態(tài)環(huán)境學(xué)報, 2013, 22(5): 879- 886.

[31] 張紅兵, 賈來喜, 李潞. SPSS寶典. 北京: 電子工業(yè)出版社, 2007: 242- 244, 277- 282.

[32] Fernández D, Barquín J,lvarez-Cabria M, Peas F J. Land-use coverage as an indicator of riparian quality. Ecological Indicators, 2014, 41: 165- 174.

[33] 曾立雄, 黃志霖, 肖文發(fā), 田耀武. 三峽庫區(qū)不同土地利用類型氮磷流失特征及其對環(huán)境因子的響應(yīng). 環(huán)境科學(xué), 2012, 33(10): 3390- 3396.

[34] 李瀟然, 李陽兵, 邵景安. 非點源污染輸出對土地利用和社會經(jīng)濟(jì)變化響應(yīng)的案例研究. 生態(tài)學(xué)報, 2016, 36(19): 6050- 6061.

[35] 孫芹芹, 黃金良, 洪華生, 馮媛. 九龍江流域城鎮(zhèn)建設(shè)用地與河流水質(zhì)關(guān)系研究. 環(huán)境科學(xué), 2011, 32(10): 2849- 2854.