呂樂婷,高曉琴,劉 琦,江 源,*

1 遼寧師范大學(xué) 地理科學(xué)學(xué)院, 大連 116029 2 北京師范大學(xué) 地理科學(xué)學(xué)部, 北京 100875

景觀格局包括景觀要素的組成和空間配置,通過影響生態(tài)過程對(duì)流域產(chǎn)生作用。景觀組成即不同的土地利用方式,其作為各種自然活動(dòng)和社會(huì)活動(dòng)的外在客觀表現(xiàn),對(duì)流域水質(zhì)有著重要的影響。不同的景觀空間配置會(huì)使污染物的發(fā)生、遷移和轉(zhuǎn)化等過程產(chǎn)生差異,從而影響河流水質(zhì)。景觀特征與河流水質(zhì)之間的關(guān)系非常復(fù)雜,景觀特征之間的相互作用和相互關(guān)系有很大的研究空間[1]。目前,景觀格局對(duì)河流水質(zhì)的影響已經(jīng)成為國內(nèi)外河流生態(tài)領(lǐng)域關(guān)注的熱點(diǎn)方向之一[2]。

隨著地理信息技術(shù)和景觀生態(tài)學(xué)的發(fā)展,景觀-水質(zhì)的研究已經(jīng)從最初的簡(jiǎn)單定性研究發(fā)展為結(jié)合各種水文、統(tǒng)計(jì)模型的定量研究[3]。劉怡娜等[4]研究發(fā)現(xiàn)：長(zhǎng)江流域農(nóng)田和城鎮(zhèn)面積百分比分別與生態(tài)系統(tǒng)水質(zhì)凈化服務(wù)存在顯著對(duì)數(shù)關(guān)系。王杰等[5]采用相關(guān)分析法、冗余分析法研究了丹江上游景觀格局與水質(zhì)的關(guān)系,發(fā)現(xiàn)耕地強(qiáng)烈影響著水質(zhì)變化,林地對(duì)水質(zhì)有顯著的正效應(yīng)。此外,許多學(xué)者[6- 8]認(rèn)同河流水質(zhì)受到空間尺度不同的影響,但是如何對(duì)空間尺度進(jìn)行界定目前還沒有統(tǒng)一的標(biāo)準(zhǔn)。

張勇榮等[9]在對(duì)筑壩河流域的研究中發(fā)現(xiàn)：子流域、河流緩沖區(qū)及監(jiān)測(cè)點(diǎn)圓形緩沖區(qū)三種空間尺度中,河流緩沖區(qū)和子流域尺度相關(guān)規(guī)律基本一致,但與圓形緩沖區(qū)尺度相關(guān)規(guī)律差異較大。方娜等[10]在對(duì)鄱陽湖濕地水質(zhì)的研究中發(fā)現(xiàn):緩沖區(qū)尺度的土地利用方式對(duì)水質(zhì)的解釋度高于小流域尺度。許多研究[3,11]表明,不同空間尺度差異對(duì)景觀-水質(zhì)關(guān)系的表現(xiàn)不同。

除了空間尺度,景觀格局在不同季節(jié)和時(shí)間對(duì)水質(zhì)影響有所差異。Li等[12]建立了我國洞庭湖典型的農(nóng)業(yè)集約小流域景觀特征與水質(zhì)的關(guān)系模型,發(fā)現(xiàn)在夏季和冬季,景觀指標(biāo)對(duì)水質(zhì)的解釋量高于春季和秋季。歐洋等[3]在對(duì)密云水庫上游流域的研究中發(fā)現(xiàn)：雨季前景觀對(duì)水質(zhì)的影響最弱,雨季中景觀對(duì)水質(zhì)的解釋能力最強(qiáng),雨季后景觀對(duì)水質(zhì)的影響減弱。目前,景觀-水質(zhì)研究大多采用空間分析、相關(guān)分析、多元回歸分析、冗余分析和方差分析等方法對(duì)不同季節(jié)、不同空間尺度(集水區(qū)、子流域、緩沖區(qū)等)下流域水質(zhì)與景觀關(guān)系進(jìn)行探究[13- 15]。

東江是珠江支流,流域中下游地段是我國東南沿海地區(qū)快速城市化的代表性區(qū)域之一[16]。近年來,城市化導(dǎo)致流域內(nèi)景觀發(fā)生了劇烈的變化,對(duì)生態(tài)環(huán)境影響深遠(yuǎn)。本文以東江流域?yàn)橹饕芯繉?duì)象,結(jié)合實(shí)地采樣測(cè)得的氮磷數(shù)據(jù),運(yùn)用ArcGIS水文分析、相關(guān)分析和冗余分析等方法,探究景觀組成和格局與流域水質(zhì)之間的關(guān)系,以集水區(qū)和緩沖區(qū)兩種空間尺度為基礎(chǔ),探究景觀格局與水質(zhì)關(guān)系,辨析不同時(shí)間和空間尺度下流域景觀格局對(duì)氮磷輸出的影響,明確景觀格局在影響東江流域水質(zhì)中的主導(dǎo)因子,結(jié)果可為東江流域的科學(xué)治理提供建議。

1 資料與方法

1.1 研究區(qū)概況

東江流域位于珠江三角洲東北端,它是深圳、香港、廣東等一線城市的供水源地。東江發(fā)源于江西省尋烏縣境內(nèi),自東北向西南流入珠江。流經(jīng)龍川、河源、紫金、惠陽、博羅、東莞等縣市,干流總長(zhǎng)563 km,流域總面積35340 km2。流域地勢(shì)北高南低,北部和中部為丘陵山地,南部為三角洲、低洼地和沿海平原(圖1)。流域內(nèi)上游多林地,山地平原多分布在中下游。東江流域?qū)儆趤啛釒Ъ撅L(fēng)氣候,年平均氣溫21℃,年平均降水量約1750 mm,降水季節(jié)分布不均勻,主要集中在4—9月。流域植被以亞熱帶常綠闊葉林為主,農(nóng)作物中北部以一年兩熟或三熟水旱輪作為主,南部以一年三熟糧食作物為主。

1.2 數(shù)據(jù)來源及處理

本研究所用水質(zhì)數(shù)據(jù)來自于2017年3月(平水期)和7月(豐水期)對(duì)70個(gè)水質(zhì)樣點(diǎn)的實(shí)地采集,水樣采集完成后低溫保存并迅速運(yùn)回實(shí)驗(yàn)室進(jìn)行參數(shù)分析。采用紫外分光法(UV2800)對(duì)總氮(WTN,water total nitrogen)進(jìn)行測(cè)定,采用鉬睇抗分光光度法(UNICO2100)對(duì)總磷(WTP,water total phosphorus)進(jìn)行測(cè)定,每個(gè)指標(biāo)測(cè)定3次后取平均值。

土地利用數(shù)據(jù)基于2015年Landsat TM/ETM影像,來源于中國科學(xué)院國際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)(http://www.cnic.cn/zcfw/sjfw/gjkxsjjx/),分辨率為30 m,經(jīng)解譯,將土地利用分為耕地、林地、草地、城鎮(zhèn)、水域、未利用地6種類型(圖2)。本研究所采用的數(shù)字高程模型(Digital Elevation Model,DEM)數(shù)據(jù)來源于美國地質(zhì)調(diào)查局網(wǎng)站(http://glovis.usgs.gov/),分辨率為30 m。水系數(shù)據(jù)采用國家基礎(chǔ)地理信息中心(http://www.ngcc.cn/ngcc/html/1/index.html)提供的1:250000水系矢量圖層。

圖2 東江流域土地利用概況Fig.2 Land use types of the Dongjiang River Basin

本研究采用ArcGIS的水文分析工具對(duì)各采樣點(diǎn)的集水區(qū)進(jìn)行嵌套式劃分,利用buffer工具,提取河岸兩側(cè)寬500 m,長(zhǎng)度覆蓋樣點(diǎn)上游河流范圍的河岸緩沖區(qū)。結(jié)合土地利用數(shù)據(jù),提取各樣點(diǎn)集水區(qū)及緩沖區(qū)的土地利用數(shù)據(jù),導(dǎo)入Fragstats軟件計(jì)算其景觀指數(shù)。

1.3 研究方法

1.3.1景觀指數(shù)的選取

結(jié)合前人對(duì)東江流域的研究以及該流域的實(shí)際情況[17- 22],本文從景觀水平和斑塊類型水平兩個(gè)方面選取了斑塊密度(Patch Density,PD)、景觀分裂指數(shù)(Landscape Division Index,DIVISION)、邊緣密度(Edge Density,ED)、蔓延度指數(shù)(Contagion,CONTAG)、最大斑塊指數(shù)(Largest Patch Index,LPI)、相似鄰近比例(Percentage of Like Adjacencies,PLADJ)、散布與并列指數(shù)(Interspersion & Juxtaposition Index,IJI)、香濃多樣性指數(shù)(Shannon′s Diversity Index,SHDI)等8個(gè)指標(biāo)。

1.3.2分析方法

本研究采用Pearson相關(guān)分析和冗余分析(Redundancy Analysis,RDA)探索景觀格局與氮、磷輸出的關(guān)系。其中相關(guān)分析可以確定變量之間的關(guān)系,常用相關(guān)系數(shù)|r|來反映變量之間相關(guān)關(guān)系的密切程度。冗余分析是一種回歸分析結(jié)合主成分分析的梯度排序分析方法,能夠從統(tǒng)計(jì)學(xué)角度評(píng)價(jià)一組變量與另一組多變量數(shù)據(jù)之間的關(guān)系。它能夠單獨(dú)解釋各個(gè)環(huán)境變量對(duì)水質(zhì)變化的貢獻(xiàn)率,并能夠有效地對(duì)多個(gè)解釋變量進(jìn)行統(tǒng)計(jì)檢驗(yàn)。環(huán)境變量與水質(zhì)參數(shù)之間的夾角大于90 °,呈負(fù)相關(guān)關(guān)系；小于90°,呈正相關(guān)關(guān)系；等于90°,則不存在關(guān)系。同時(shí),線段長(zhǎng)度反映了環(huán)境變量對(duì)水質(zhì)影響程度的強(qiáng)弱[23]。

2 結(jié)果與分析

2.1 景觀格局與水質(zhì)相關(guān)分析結(jié)果

對(duì)兩次采樣水質(zhì)總氮(WTN, water total nitrogen)、總磷(WTP, water total phosphoru)數(shù)據(jù)、景觀格局指數(shù)進(jìn)行Pearson相關(guān)分析。結(jié)果如表1所示。

表1 不同時(shí)空尺度土地利用和景觀格局與氮、磷輸出的相關(guān)分析

在集水區(qū)：LPI、CONTAG與河流氮磷含量都表現(xiàn)出顯著負(fù)相關(guān)；IJI、DIVISION、SHDI與河流氮磷含量表現(xiàn)出顯著正相關(guān)。ED、PD與河流氮磷含量表現(xiàn)出正相關(guān),PLADJ與河流氮磷含量表現(xiàn)出負(fù)相關(guān),但它們都與7月份總磷相關(guān)性不顯著。從土地利用角度來看：林地與河流氮磷含量表現(xiàn)出強(qiáng)負(fù)相關(guān),城鎮(zhèn)和水域與河流氮磷含量呈強(qiáng)正相關(guān)。未利用地在3月份與河流氮磷含量正相關(guān),7月份相關(guān)性不顯著。耕地與河流氮磷含量正相關(guān)、草地與河流氮磷含量負(fù)相關(guān),但相關(guān)性都很微弱。

在緩沖區(qū)：土地利用對(duì)河流水質(zhì)的影響比景觀格局強(qiáng),3月份的影響比7月份顯著。CONTAG與河流氮磷含量顯著負(fù)相關(guān),IJI與河流氮磷含量顯著正相關(guān)；PD、ED和SHDI等指標(biāo)與3月份河流總氮含量顯著正相關(guān),PLADJ與3月份河流總氮含量顯著負(fù)相關(guān),但7月份的相關(guān)性都不顯著。LPI、DIVISION與河流水質(zhì)相關(guān)性不顯著。從土地利用角度來看：林地、草地與河流水質(zhì)呈顯著負(fù)相關(guān),城鎮(zhèn)與河流水質(zhì)顯著正相關(guān),水域與河流水質(zhì)表現(xiàn)出正相關(guān)關(guān)系,其中7月份水域與河流總氮含量無關(guān),其余土地利用方式與河流水質(zhì)相關(guān)性微弱。

總體上,景觀組成與格局對(duì)河流總氮含量的影響大于總磷。就相關(guān)系數(shù)|r|來看,大部分集水區(qū)的相關(guān)系數(shù)比緩沖區(qū)大,3月份的相關(guān)系數(shù)比7月份大。

2.2 景觀格局與水質(zhì)冗余分析結(jié)果

本研究首先對(duì)各樣本的水質(zhì)指標(biāo)進(jìn)行除趨勢(shì)對(duì)應(yīng)分析(DCA分析),結(jié)果顯示所有排序軸中的最長(zhǎng)梯度值為0.988(小于3.0),所以采用線性模型(RDA)[23]。樣本與景觀指數(shù)RDA分析結(jié)果如表2所示。

表2 土地利用、景觀格局對(duì)水質(zhì)的解釋量

在集水區(qū)尺度,前兩個(gè)排序軸特征根為0.445和0.016。DIVISION和SHDI對(duì)氮磷含量變化的解釋量較高,分別為0.31和0.30；其次為L(zhǎng)PI、ED、CONTAG、PLADJ、林地和城鎮(zhèn)；耕地、未利用地和草地對(duì)氮磷含量變化的解釋量最小,其值不足0.1。在緩沖區(qū)尺度,前兩個(gè)排序軸的特征根為0.415和0.026。林地和城鎮(zhèn)對(duì)氮磷含量變化的解釋量最高,分別為0.26和0.24；其次為CONTAG、IJI、SHDI、草地和水域,其解釋量在0.1以上；其余指標(biāo)對(duì)氮磷含量變化的解釋量不足0.1?？傮w上來看,流域?qū)Φ缀孔兓慕忉屃扛哂诰彌_區(qū),即流域尺度對(duì)氮磷含量變化的影響更大,這與相關(guān)分析所得出的結(jié)果一致。

圖3 集水區(qū)和緩沖區(qū)景觀格局與水質(zhì)的RDA排序 Fig.3 Redundancy Analysis (RDA) ranking of land use and landscape patterns with water quality in catchment areas and buffer zonePD：斑塊密度,Patch Density；DIVISION：景觀分裂指數(shù),Landscape Division Index；ED：邊緣密度Edge Density；CONTAG：蔓延度指數(shù),Contagion；LPI：最大斑塊指數(shù),Largest Patch Index；PLADJ：相似鄰近比例,Percentage of Like Adjacencies；IJI：散布與并列指數(shù),Interspersion & Juxtaposition Index；SHDI：香濃多樣性指數(shù),Shannon′s Diversity Index

RDA排序結(jié)果如圖3所示。在集水區(qū),IJI、PD、SHDI、DIVISION、ED等景觀指數(shù)與總氮總磷含量呈正相關(guān)關(guān)系,LPI、PLADJ、CONTAG等景觀指數(shù)與總氮總磷含量呈負(fù)相關(guān)關(guān)系；耕地與氮磷含量關(guān)系不大,水域、城鎮(zhèn)與氮磷含量呈正相關(guān)關(guān)系。林地、草地與氮磷含量呈負(fù)相關(guān)關(guān)系。如圖3所示,在緩沖區(qū),IJI、SHDI、PD、DIVISION、ED以及未利用地、水域、城鎮(zhèn)等指標(biāo)與河流總氮總磷含量表現(xiàn)出強(qiáng)正相關(guān)；LPI、PLADJ、CONTAG和林地、草地等指標(biāo)與河流總氮總磷含量表現(xiàn)出強(qiáng)負(fù)相關(guān)；耕地與氮磷含量的相關(guān)性不強(qiáng)。在兩組排序圖中都可以看出總氮的箭頭較長(zhǎng),且3月的箭頭比7月的長(zhǎng),表明土地利用及景觀指數(shù)對(duì)總氮的影響程度較大,且對(duì)3月份比7月份強(qiáng)。

3 討論

3.1 景觀組成與氮磷含量關(guān)系

在不同的時(shí)空尺度,東江流域河流氮磷含量與城鎮(zhèn)面積呈正相關(guān)關(guān)系,城鎮(zhèn)是河流氮磷的主要貢獻(xiàn)源。這與夏品華等[24]、羅璇等[25]等的研究結(jié)論一致。在東江流域,城鎮(zhèn)面積不斷增加且表現(xiàn)出集中化趨勢(shì)[22]。居民區(qū)、道路以及建設(shè)用地等不透水面的增加導(dǎo)致下滲減少,污染物直接進(jìn)入河流；居民的生產(chǎn)、生活活動(dòng)產(chǎn)生大量的污水未經(jīng)處理排入水體；此外,居民區(qū)產(chǎn)生的各類垃圾得不到妥善處理,通過地表徑流沖刷進(jìn)入河流,對(duì)水體造成污染,導(dǎo)致河道氮磷含量升高。

東江流域河流眾多,水域(包括灘涂、水庫、湖泊、基塘和溝渠等)分布廣泛。流域內(nèi)養(yǎng)殖業(yè)數(shù)量多,排污量大。且流域內(nèi)畜禽養(yǎng)殖業(yè)管理粗放,多采用水沖糞工藝,未配套污水處理設(shè)施。當(dāng)某一水體受污染后就會(huì)隨水流擴(kuò)展到周邊河網(wǎng)。此外,農(nóng)村勞動(dòng)力外出打工,造成基塘缺乏管理,水體富營養(yǎng)化[26],也加劇了水體中氮磷含量。因此水域面積與河道氮磷含量呈顯著正相關(guān)關(guān)系。劉旭攏等[27]的研究結(jié)果也證實(shí)了這一點(diǎn)。

在東江流域,河流氮磷含量與耕地之間沒有顯著的相關(guān)性。這與丁佼[17]、查智琴[28]、夏品華[24]等人研究結(jié)論相同。研究區(qū)渠網(wǎng)交錯(cuò),農(nóng)業(yè)活動(dòng)以水田為主。水田作為人工濕地,可以有效截留流域內(nèi)氮磷物質(zhì)[29],凈化水質(zhì)。此外,各種土地利用類型在空間上具有消長(zhǎng)關(guān)系,且研究區(qū)城市化水平高,人口稠密,經(jīng)濟(jì)發(fā)達(dá),城鎮(zhèn)點(diǎn)源污染嚴(yán)重,相比之下,耕地對(duì)河流水質(zhì)的影響不顯著。

東江流域河流氮磷含量與林地、草地呈負(fù)相關(guān)關(guān)系。這與楊婭楠等[30]、王晶萍等[7]學(xué)者等的觀點(diǎn)一致。林地和草地都有涵養(yǎng)水源,保持水土的功效,對(duì)污染物起到吸附、滯留和過濾的作用,可以有效阻滯污染物向河流運(yùn)輸。研究區(qū)位于我國東南部,草木茂盛,林地面積占比約70%左右,多分布在流域中上游且斑塊集中成片,連通性高,有利于對(duì)氮磷等污染物的截留和阻攔。

3.2 景觀格局與氮磷含量關(guān)系

東江流域氮磷污染受景觀整體性的影響。分析結(jié)果顯示,DIVISION和IJI這兩個(gè)指標(biāo)都與河流總氮總磷含量呈正相關(guān)關(guān)系,這與宮殿林等[13]、蔡宏等[31]學(xué)者的研究結(jié)果一致。東江流域上游景觀分裂度低,中下游地區(qū)景觀分裂度高[22]。這是因?yàn)榱饔蛏嫌味嗌降厍鹆?人類活動(dòng)受限；中下游城市擴(kuò)張迅速,經(jīng)濟(jì)以高新技術(shù)產(chǎn)業(yè)為主,居民生產(chǎn)生活對(duì)基礎(chǔ)設(shè)施的需求大,導(dǎo)致對(duì)自然景觀的開發(fā)強(qiáng)度增加,景觀整體性受損,河流氮磷污染也隨之增加。

SHDI高值表明景觀多樣性豐富,ED與PD高值說明景觀破碎程度大。東江流域SHDI、ED與PD皆與河流總氮總磷含量呈正相關(guān)關(guān)系,意味著隨著人類活動(dòng)干擾強(qiáng)度增加,河流氮磷污染加重。這與劉曉君等[32]、普軍偉等[33]、曹燦等[14]學(xué)者的研究結(jié)果一致。東江流域景觀多樣性呈上游小,下游大的趨勢(shì),且高值區(qū)為經(jīng)濟(jì)發(fā)達(dá)的城鎮(zhèn)地區(qū)[22]。這說明隨著快速城鎮(zhèn)化,為滿足需求,原有的自然景觀斑塊被割裂成若干多元小斑塊,景觀多樣性和破碎度增加。例如道路以及河道等基礎(chǔ)設(shè)施的修建；耕地、林地等轉(zhuǎn)移為建設(shè)用地,都會(huì)導(dǎo)致景觀破碎程度增加,連通度下降,河流氮磷污染也隨之加重。

反之,LPI值高意味著景觀整體性好,破碎度低,人類干擾活動(dòng)少。這與吉冬青[34]等學(xué)者的研究結(jié)論一致。東江流域上游多山地丘陵,植被覆蓋度高,整體景觀格局較完整,最大斑塊指數(shù)較高,人類活動(dòng)對(duì)自然環(huán)境的開發(fā)程度有限,氮磷污染輸出少；流域下游為三角洲和沿江平原,地形平坦,城市擴(kuò)張迅速,最大斑塊值相對(duì)較低[22],人為因素對(duì)景觀的影響突出,氮磷污染輸出較多。因此,LPI在不同季節(jié)都與河流氮磷含量呈負(fù)相關(guān)關(guān)系。

CONTAG和PLADJ表示斑塊集聚程度和景觀面積比,這兩個(gè)指標(biāo)在東江流域都與氮磷含量負(fù)相關(guān)。流域上游以林地為主,景觀團(tuán)聚程度較好,斑塊連通性高,有助于對(duì)污染物質(zhì)的截留和過濾。流域下游鄰近粵港澳大灣區(qū),城市擴(kuò)張迅速,對(duì)自然景觀的干擾性大,景觀破碎,異質(zhì)性高,各斑塊分散且連通性差,導(dǎo)致氮磷輸出增加。這與徐建鋒等[35]的研究結(jié)果一致。

3.3 不同時(shí)空尺度下景觀格局與氮磷含量關(guān)系

本研究結(jié)果顯示：東江流域3月份河流氮磷含量與景觀指數(shù)的相關(guān)性比7月份高,即平水期高于豐水期,郭玉靜等[36]的研究也證實(shí)了這一結(jié)論。究其原因可能是由于豐水期降水多,沖刷能力強(qiáng),河流流量和流速增加,污染物質(zhì)停留時(shí)間短,各種植物的吸附過濾作用受到影響,且豐水期流量大,雨水對(duì)污染物質(zhì)起到稀釋作用,使得污染物濃度降低。平水期降雨相對(duì)較少,匯流減少,工業(yè)廢水和生活污水長(zhǎng)期穩(wěn)定排放,其對(duì)氮磷輸出的貢獻(xiàn)凸顯。此時(shí),河流流量、流速相對(duì)穩(wěn)定,污染物在水體中維持時(shí)間較長(zhǎng)。

在空間尺度,東江流域集水區(qū)景觀格局與河流氮磷含量的相關(guān)性更高,即流域尺度的景觀格局對(duì)河流氮磷污染的影響更大。此結(jié)論與張軍等[37]在對(duì)丹江流域的研究結(jié)果一致。另外,氮素受控于不同的土地利用方式上的人類活動(dòng),而磷的影響因素更復(fù)雜多樣,因此總氮含量與景觀指數(shù)的相關(guān)性要高于總磷,趙軍等[38]、Molinero等[39]的研究也得到了相同結(jié)論。

4 結(jié)論

本研究采用相關(guān)分析和冗余分析,對(duì)東江流域景觀格局與氮、磷輸出的關(guān)系進(jìn)行了不同空間尺度和季節(jié)變化的研究,結(jié)論如下：

(1)城鎮(zhèn)、水域?qū)恿鞯孜廴居胸?fù)面影響,其中城鎮(zhèn)是氮、磷輸出的主要來源。林地、草地面積對(duì)河流氮磷含量有正面影響。

(2)景觀破碎度、分離度與多樣性與河流氮磷輸出呈正相關(guān)關(guān)系,景觀聚散性與優(yōu)勢(shì)度與河流氮磷輸出呈負(fù)相關(guān)關(guān)系。

(3)平水期河流氮磷含量與景觀格局的相關(guān)性強(qiáng)于豐水期；集水區(qū)尺度河流氮磷含量與景觀格局相關(guān)性高于緩沖區(qū)尺度；此外,總氮對(duì)景觀格局的響應(yīng)比總磷更敏感。

參考文獻(xiàn)(References):

[1] Lintern A, Webb J A, Ryu D, Liu S, Bende-Michl U, Waters D, Leahy P, Wilson P, Western A W. Key factors influencing differences in stream water quality across space. Wiley Interdisciplinary Reviews Water, 2018, 5(1): e1260.

[2] 劉麗娟, 李小玉, 何興元. 流域尺度上的景觀格局與河流水質(zhì)關(guān)系研究進(jìn)展. 生態(tài)學(xué)報(bào), 2011, 31(19): 5460- 5465.

[3] 歐洋, 王曉燕, 耿潤(rùn)哲. 密云水庫上游流域不同尺度景觀特征對(duì)水質(zhì)的影響. 環(huán)境科學(xué)學(xué)報(bào), 2012, 32(5): 1219- 1226.

[4] 劉怡娜, 孔令橋, 肖燚, 鄭華. 長(zhǎng)江流域景觀格局變化對(duì)生態(tài)系統(tǒng)水質(zhì)凈化服務(wù)的影響. 環(huán)境保護(hù)科學(xué), 2018, 44(4): 6- 13.

[5] 王杰, 李鵬, 高海東, 時(shí)鵬, 張秦嶺, 楊倩楠, 馬勇勇. 丹江上游土地利用/景觀指數(shù)與水質(zhì)關(guān)系初探. 水土保持研究, 2018, 25(6): 383- 389.

[6] 鄧靈稚. 基于GWR模型的景觀格局與水環(huán)境相關(guān)性分析[D]. 重慶: 重慶師范大學(xué), 2018.

[7] 王晶萍, 李兆富, 劉紅玉, 王剛, 辛強(qiáng). 中田河流域景觀異質(zhì)性對(duì)水體總氮濃度影響研究. 環(huán)境科學(xué), 2016, 37(2): 527- 533.

[8] 楊傳鳳. 基于“源—匯”格局的山區(qū)流域面源污染評(píng)價(jià)與分區(qū)管理[D]. 泰安: 山東農(nóng)業(yè)大學(xué), 2017.

[9] 張勇榮, 周忠發(fā), 蔣翼, 尚夢(mèng)佳. 筑壩河流不同尺度土地利用結(jié)構(gòu)對(duì)庫區(qū)水質(zhì)的影響——以平寨水庫為例. 水利水電技術(shù), 2019, 50(4): 138- 145.

[10] 方娜, 劉玲玲, 游清徽, 田娜, 吳燕平, 陽文靜. 不同尺度土地利用方式對(duì)鄱陽湖濕地水質(zhì)的影響. 環(huán)境科學(xué), 2019, 40(12): 5348- 5357.

[11] 李艷利, 李艷粉, 徐宗學(xué), 關(guān)中美. 渾太河上游流域河岸緩沖區(qū)景觀格局對(duì)水質(zhì)的影響. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2015, 31(1): 59- 68.

[12] Li H Q, Liu L M, Ji X. Modeling the relationship between landscape characteristics and water quality in a typical highly intensive agricultural small watershed, Dongting Lake Basin, South Central China. Environmental Monitoring and Assessment, 2015, 187(3): 129.

[13] 宮殿林, 洪曦, 曾冠軍, 王毅, 左雙苗, 劉新亮, 吳金水. 亞熱帶典型農(nóng)業(yè)流域河流水質(zhì)多元線性回歸預(yù)測(cè). 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2017, 33(6): 509- 518.

[14] 曹燦, 張飛, 阿依尼格爾·亞力坤, 朱世丹, 郭苗, 阿麗米熱·塔力甫江, Hsiangte K. 艾比湖區(qū)域景觀格局與河流水質(zhì)關(guān)系探討. 環(huán)境科學(xué), 2018, 39(4): 1568- 1577.

[15] 王小平, 張飛, 李曉航, 曹燦, 郭苗, 陳麗華. 艾比湖區(qū)域景觀格局空間特征與地表水質(zhì)的關(guān)聯(lián)分析. 生態(tài)學(xué)報(bào), 2017, 37(22): 7438- 7452.

[16] 彭秋志, 秦國玲, 呂樂婷, 吳亞玲. 基于MODIS數(shù)據(jù)的東江流域云干擾時(shí)空特征分析. 國土資源遙感, 2018, 30(1): 109- 115.

[17] 丁佼. 東江流域土地利用與地表水質(zhì)關(guān)系的景觀生態(tài)學(xué)研究[D]. 北京: 北京師范大學(xué), 2016.

[18] 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.

[19] 周婷, 彭少麟, 任文韜. 東江河岸緩沖帶景觀格局變化對(duì)水體恢復(fù)的影響. 生態(tài)學(xué)報(bào), 2009, 29(1): 231- 239.

[20] 周婷, 彭少麟, 任文韜. 東江河岸帶各區(qū)段景觀格局對(duì)水體恢復(fù)的影響. 中山大學(xué)學(xué)報(bào): 自然科學(xué)版, 2009, 48(4): 78- 83.

[21] 任文韜, 彭少麟, 周婷, 李巖. 東江流域集水區(qū)城市化差異及其對(duì)景觀格局的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2008, 19(12): 2680- 2686.

[22] 呂樂婷, 張杰, 彭秋志, 任斐鵬, 江源. 東江流域景觀格局演變分析及變化預(yù)測(cè). 生態(tài)學(xué)報(bào), 2019, 39(18): 6850- 6859.

[24] 夏品華, 孔祥量, 喻理飛. 草海濕地小流域土地利用與景觀格局對(duì)氮、磷輸出的影響. 環(huán)境科學(xué)學(xué)報(bào), 2016, 36(8): 2983- 2989.

[25] 羅璇, 史志華, 尹煒, 陳佳, 李璐, 伍世良. 小流域土地利用結(jié)構(gòu)對(duì)氮素輸出的影響. 環(huán)境科學(xué), 2010, 31(1): 58- 62.

[26] 王文靜, 韓寶龍, 鄭華, 歐陽志云. 粵港澳大灣區(qū)生態(tài)系統(tǒng)格局變化與模擬. 生態(tài)學(xué)報(bào), 2020, 40(10): 3364- 3374.

[27] 劉旭攏, 鄧孺孺, 秦雁, 梁業(yè)恒, 熊龍海. 東江流域地表水功能區(qū)水質(zhì)對(duì)土地利用的響應(yīng). 熱帶地理, 2016, 36(2): 296- 302.

[28] 查智琴. 哈尼梯田濕地水體氮磷時(shí)空盈虧及其關(guān)鍵匯景觀識(shí)別研究[D]. 昆明: 云南師范大學(xué), 2019.

[29] 查智琴, 角媛梅, 劉志林, 劉澄靜, 高璇, 丁銀平, 趙冬梅. 哈尼梯田濕地景觀水體富營養(yǎng)化及截留效應(yīng)評(píng)價(jià). 生態(tài)學(xué)雜志, 2018, 37(11): 3413- 3421.

[30] 楊婭楠, 王金亮, 陳光杰, 習(xí)曉環(huán), 王成. 撫仙湖流域土地利用格局與水質(zhì)變化關(guān)系. 國土資源遙感, 2016, 28(1): 159- 165.

[31] 蔡宏, 林國敏, 康文華. 赤水河流域中上游坡地景觀特征對(duì)河流水質(zhì)的影響. 地理研究, 2018, 37(4): 704- 716.

[32] 劉曉君, 李占斌, 李鵬, 張鐵鋼, 徐國策, 高海東. 基于土地利用/覆被變化的流域景觀格局與水沙響應(yīng)關(guān)系研究. 生態(tài)學(xué)報(bào), 2016, 36(18): 5691- 5700.

[33] 普軍偉, 趙筱青, 顧澤賢, 盧飛飛, 鄭田甜, 高翔宇, 陳俊旭, 易琦. 云南高原杞麓湖流域的景觀格局與水質(zhì)變化. 水生態(tài)學(xué)雜志, 2018, 39(5): 13- 21.

[34] 吉冬青, 文雅, 魏建兵, 吳志峰, 劉慶, 程炯. 流溪河流域景觀空間特征與河流水質(zhì)的關(guān)聯(lián)分析. 生態(tài)學(xué)報(bào), 2015, 35(2): 246- 253.

[35] 徐建鋒, 尹煒, 閆峰陵, 史志華. 農(nóng)業(yè)源頭流域景觀異質(zhì)性與溪流水質(zhì)耦合關(guān)系. 中國環(huán)境科學(xué), 2016, 36(10): 3193- 3200.

[36] 郭玉靜, 王妍, 劉云根, 鄭毅, 張超, 侯磊. 普者黑巖溶湖泊濕地湖濱帶景觀格局演變對(duì)水質(zhì)的影響. 生態(tài)學(xué)報(bào), 2018, 38(5): 1711- 1721.

[37] 張軍. 丹江流域植被格局演變及其與水質(zhì)響應(yīng)關(guān)系研究[D]. 西安: 西安理工大學(xué), 2017.

[38] 趙軍, 楊凱, 邰俊, 單福征. 區(qū)域景觀格局與地表水環(huán)境質(zhì)量關(guān)系研究進(jìn)展. 生態(tài)學(xué)報(bào), 2011, 31(11): 3180- 3189.

[39] Molinero J, Burke R A. Relations between land use and stream nutrient concentration for small watersheds in the Georgia piedmont//Proceedings of the 2003 Georgia Water Resources Conference. Georgia: Georgia Water Resources Insitute, 2003.