沮河流域景觀格局對水體氮磷輸出的影響

摘要：探究沮河土地利用與景觀格局對氮磷輸出的影響，為當(dāng)?shù)厮Y源和生態(tài)系統(tǒng)保護(hù)以及沮河流域長期可持續(xù)發(fā)展提供數(shù)據(jù)支撐。2023年2月（枯水期）和7月（豐水期）開展沮河流域水質(zhì)調(diào)查，運(yùn)用遙感影像數(shù)據(jù)解譯沮河流域土地利用現(xiàn)狀的空間分布以及計(jì)算流域景觀格局指數(shù)，采用相關(guān)分析和冗余分析方法，從時(shí)間和空間2個(gè)維度，探討沮河子流域景觀格局與出口處氮、磷輸出關(guān)系。結(jié)果表明：（1）磷污染最大的影響因子是磷礦，枯水期氮磷污染主要來源于果園、農(nóng)田和城鎮(zhèn)；（2）近圓形狀指數(shù)面積加權(quán)平均值對流域氮磷污染物的截留有著正向作用；（3）枯水期，礦場和水體的歐氏鄰近距離面積加權(quán)平均值，森林回旋半徑面積加權(quán)平均值、草地對比度加權(quán)邊緣密度，農(nóng)田邊緣對比指數(shù)平均值是影響沮河水體中氮磷濃度變化的主要因子；在豐水期，農(nóng)田和森林邊緣對比指數(shù)平均值、農(nóng)田和礦場的對比度加權(quán)邊緣密度、森林近圓形狀指數(shù)面積加權(quán)平均值、礦場歐氏鄰近距離面積加權(quán)平均值、森林鄰近指數(shù)面積加權(quán)平均值是影響沮河水體中氮磷濃度變化的主要因子。研究結(jié)果對控制沮河水體氮磷濃度的變化有一定指導(dǎo)意義。

關(guān)鍵詞：景觀格局；土地利用；氮磷輸出；沮河流域

中圖分類號：P901；X143" " " " 文獻(xiàn)標(biāo)志碼：A" " " " 文章編號：1674-3075（2025）02-0062-10

景觀格局包括景觀要素的組成和空間配置，通過影響生態(tài)過程對流域產(chǎn)生作用（Turner，1989）。不同的土地利用方式作為各種自然活動和社會活動的外在客觀表現(xiàn)，對流域水質(zhì)有著重要的影響（Chen et al，2008; Liu et al，2020）。不同的景觀空間配置會使污染物的發(fā)生、遷移和轉(zhuǎn)化等過程產(chǎn)生差異，從而影響河流水質(zhì)（Xu et al，2021）。目前，景觀格局對河流水質(zhì)的影響已經(jīng)成為國內(nèi)外河流生態(tài)領(lǐng)域關(guān)注的熱點(diǎn)方向之一。流域土地利用和景觀格局對氮磷輸出有著較大影響（Carey et al，2011；袁軼男等，2020；祁蘭蘭等，2021；姚靜，2022），如農(nóng)業(yè)、工業(yè)和城市化等對水體富營養(yǎng)化的影響不可忽視（Savage et al，2010；宋夢廷等，2022）。農(nóng)業(yè)活動（例如施肥、灌溉和農(nóng)作物種植等）是氮磷輸出的主要來源之一（沈濤等，2007；Chen et al，2008；Cao et al，2014）。此外，城市化和工業(yè)化（如工業(yè)廢水、生活污水的排放，以及交通用地的增加等）也會對流域水質(zhì)產(chǎn)生影響，（Teng，2011；Bartarya et al，2012；黃益平；Sarker et al，2021等，2021）。

對于不同空間尺度下流域景觀格局對氮磷輸出的影響仍存在爭議（Zhou et al，2012；蔡宏等，2015；Shen et al，2015；黃益平等，2021）。一些研究表明，在集水區(qū)和緩沖區(qū)等不同空間下，景觀格局對氮磷輸出的影響存在差異（黃益平等，2021；Liu et al，2022）。蔡宏等（2015）發(fā)現(xiàn)緩沖區(qū)尺度的土地利用方式對水質(zhì)的解釋度與子流域尺度表現(xiàn)出一致的規(guī)律。這些不一致的結(jié)論可能是不同研究采用的空間尺度、土地利用方式、水質(zhì)指標(biāo)以及數(shù)據(jù)采集時(shí)間等因素存在差異所導(dǎo)致。除了空間尺度，景觀格局在不同季節(jié)下對水體氮磷輸出影響也有所差異。彭勃等（2019）運(yùn)用GIS技術(shù)及SPSS統(tǒng)計(jì)分析方法探討了枯水期與豐水期土地景觀格局對河流水質(zhì)的影響，結(jié)果顯示枯水期較豐水期影響更為顯著。呂樂婷等（2021）在東江流域采用相關(guān)分析和冗余分析發(fā)現(xiàn)枯水期河流水質(zhì)與景觀格局的相關(guān)性強(qiáng)于豐水期。因此，研究不同季節(jié)下流域景觀格局對氮磷輸出的影響十分必要。

沮河是沮漳河的主脈，為長江中游北岸一條支流，水系不太發(fā)育，兩岸支流短促，全長226 km，河道平均比降0.189%，集水面積 3 382 km2。近年來，隨著社會經(jīng)濟(jì)的發(fā)展和城市化進(jìn)程的加速，沮漳河流域的生態(tài)環(huán)境受到嚴(yán)重影響（陳崇德等，2008；黃永文等，2013）。為此，本文以沮漳河的主脈沮河為對象，探索土地利用與景觀格局對氮磷輸出的影響，分析土地利用和景觀格局對沮河流域水質(zhì)的影響，為沮河流域水環(huán)境保護(hù)、土地利用結(jié)構(gòu)改善和景觀格局優(yōu)化提供一定的科學(xué)依據(jù)和借鑒。

1" "材料與方法

1.1" "區(qū)域概況

研究區(qū)域內(nèi)整體地勢特點(diǎn)為西高東低，西部以山地、丘陵為主，東部以平原、丘陵為主；土地利用類型主要為林地和農(nóng)田；流域?qū)賮啛釒Ъ撅L(fēng)性濕潤氣候，四季分明，雨熱資源豐富，年平均降水量在1 100 mm左右，降雨主要集中在6、7月份。研究區(qū)自峽口水庫以上部分為上游，礦產(chǎn)資源豐富，經(jīng)濟(jì)林、用材林種植面積較廣；峽口水庫以下至百草壩水庫部分為中游，以水稻、油菜、茶葉、柑橘等經(jīng)濟(jì)作物為主；自百草壩水庫以下部分則為下游，水資源優(yōu)勢極大，其水產(chǎn)品及養(yǎng)殖業(yè)也有較大發(fā)展，下游地區(qū)城鎮(zhèn)化、工業(yè)化擴(kuò)展迅速，下游生態(tài)環(huán)境和水質(zhì)不斷下降（黃永文等，2013）。

1.2" "流域土地利用類型確定

本文所用土地利用現(xiàn)狀數(shù)據(jù)來源于中國資源衛(wèi)星應(yīng)用中心（https：//data.cresda.cn/）2022年分辨率為16 m且云量接近于0的環(huán)境二號遙感影像。使用ENVI5.3進(jìn)行正射校正、輻射定標(biāo)、大氣校正后為確保影像解譯精度，在ArcMap10.7中選取50個(gè)控制點(diǎn)進(jìn)行地理配準(zhǔn)。配準(zhǔn)后的結(jié)果在ENVI5.3中進(jìn)行監(jiān)督分類將土地利用劃分為農(nóng)田、水體、林地、草地、城鎮(zhèn)、果園、礦場7種利用類型（圖1）。在Google Earth pro中選取驗(yàn)證樣本與分類結(jié)果進(jìn)行精度驗(yàn)證，精度評價(jià)結(jié)果Kappa系數(shù)為0.85，符合分類精度要求。

研究區(qū)域土地利用類型主要以森林與農(nóng)田為主，研究區(qū)域上游地區(qū)平均海拔高于中下游上百米，此區(qū)域主要為山地，土地利用類型以森林與農(nóng)田為主且水體分布較少，土地利用的類型較單一。由于地勢原因，農(nóng)田主要以坡地和梯田為主，斑塊破碎度較高。自峽口水庫以下土地利用類型豐富，破碎度較高，水體分布較好。研究區(qū)域中游的地勢呈兩邊高中間低，河流在此匯集，從中流過農(nóng)田、果園、草地、城鎮(zhèn)。研究區(qū)域下游，自向家草壩以下地勢開始逐漸平緩，屬于長江中游平原，地勢平緩適宜居住及大面積種植，土地利用類型以農(nóng)田、城鎮(zhèn)、草地為主（圖1）。

1.3" "采樣點(diǎn)分布與氮磷測定

本研究采自地理空間數(shù)據(jù)云網(wǎng)站（https：//www.gscloud.cn/）中30 m分辨率的數(shù)字高程模型（DEM）。通過SWAT插件提取水系及子流域劃分，結(jié)合劃分結(jié)果中的子流域出水口和實(shí)地情況確定10個(gè)采樣點(diǎn)以及與采樣點(diǎn)序號對應(yīng)的10個(gè)子流域（圖2）。1號采樣區(qū)域以森林和水體為主；2號采樣點(diǎn)位于峽口柑橘種植區(qū)域的最下游，此流域的果園面積在所有子流域中占比最大；3號流域內(nèi)有一處地下磷礦廠及一座磷礦加工廠；4號流域緊挨3號流域的磷礦加工廠；5號、8號流域內(nèi)土地利用類型較全，作為正常控制斷面；6號流域位于遠(yuǎn)安縣城下游，此流域土地利用類型主要以城鎮(zhèn)為主；7號采樣點(diǎn)位于向家草壩水庫，該流域同樣以森林與水體為主；9號采樣點(diǎn)位于當(dāng)陽市郊區(qū)，該子流域主要匯集了農(nóng)田污染物；10號采樣點(diǎn)處于研究區(qū)域末端，作為最后一個(gè)控制斷面。其中，3～6號采樣點(diǎn)位于遠(yuǎn)安縣境內(nèi)，7～10號采樣點(diǎn)位于當(dāng)陽市境內(nèi)。

2023年2月（枯水期）和7月（豐水期）進(jìn)行實(shí)地采樣，在河心和距岸邊2 m處0.1 m水深取500 mL水體，混合后放入500 mL采樣瓶中，每個(gè)樣點(diǎn)三瓶，低溫避光保存。在實(shí)驗(yàn)室采用鉬酸銨分光光度法（GB 11893—89）測定總磷（TP）濃度，總氮（TN）濃度測定則采用堿性過硫酸鉀消解紫外分光光度法（GB 11894—89）。

1.4" "景觀指數(shù)提取

將分類好的土地利用類型圖導(dǎo)入FRAGSTAT4.2軟件，選擇指標(biāo)、設(shè)定景觀指數(shù)參數(shù)、輸入類型描述表后運(yùn)行模型，將運(yùn)行結(jié)果與豐水期和枯水期的TN、TP數(shù)據(jù)進(jìn)行相關(guān)性分析。最終分別選取與TN、TP數(shù)據(jù)呈正負(fù)相關(guān)的前3位景觀指數(shù)，近圓形狀指數(shù)面積加權(quán)平均值（CIRCLE AM）、對比度加權(quán)邊緣密度（CWED）、邊緣對比指數(shù)平均值（ECON MN）、鄰近指數(shù)面積加權(quán)平均值（PROX AM）、歐氏鄰近距離面積加權(quán)平均值（ENN AM），回旋半徑面積加權(quán)平均值（GYRATE AM）等6個(gè)指標(biāo)（表1）。

1.5" "分析方法

本研究采用Pearson相關(guān)分析和冗余分析（Redundancy Analysis，RDA）探討不同土地利用景觀格局對氮磷輸出的影響。Pearson相關(guān)分析確定變量之間的關(guān)系，用相關(guān)系數(shù)R反映變量之間的相關(guān)程度。冗余分析（RDA）既可以揭示單個(gè)景觀變量對河流水質(zhì)變化的貢獻(xiàn)率，又可以用二維排序圖直觀展示景觀指數(shù)與水質(zhì)參數(shù)之間的關(guān)系。消除趨勢對應(yīng)分析（DCA）是一種處理長生態(tài)梯度數(shù)據(jù)的分析方法，在運(yùn)用RDA模型前先通過DCA分析，根據(jù)分析結(jié)果顯示的排序軸大小選擇模型，數(shù)值小于4.0時(shí)則采用RDA模型。

2" "結(jié)果與分析

2.1" "氮磷輸出的時(shí)空差異

從圖3中可以看出，不同的時(shí)間和空間的氮磷輸出量存在差異。從時(shí)間上來看，豐水期的TN、TP濃度明顯高于枯水期含量。同時(shí)，兩地豐水期的降雨量高于枯水期（圖4）。豐水期TN、TP的最高濃度分別達(dá)到了19.03 mg/L和0.61 mg/L。TN、TP在枯水期的最高濃度只有6.93 mg/L和0.17 mg/L（圖3）。

從空間上來看，擁有磷礦場占比最大的3號流域輸出的TP含量在豐水期和枯水期都達(dá)到最高值（圖3，表2），分別為0.61 mg/L和0.17 mg/L。TP濃度其次的是含有包含大范圍工礦企業(yè)、城鎮(zhèn)面積占比第2的6號流域，其在豐水期和枯水期TP均排第2高值（表2）。且6號流域的TN濃度在豐水期達(dá)到最高值19.03 mg/L，其余流域的濃度差異性較小，均在17、18 mg/L左右?？菟诘腡N差值相對較小，最高值出現(xiàn)在果園面積最大的2號流域。枯水期最低的TP濃度只有0.01 mg/L，出現(xiàn)在擁有大型水庫、大面積森林及城鎮(zhèn)面積比較低的1號流域。

2.2" "豐水期流域景觀格局與氮磷冗余分析

此次研究將氮磷數(shù)據(jù)與兩個(gè)時(shí)期樣本進(jìn)行消除趨勢對應(yīng)分析，豐水期結(jié)果顯示排序軸中最長梯度值為0.9（小于3.0），所以采用RDA分析，結(jié)果如圖5所示。

RDA分析結(jié)果顯示第一軸解釋率為28.7%，景觀指數(shù)對水質(zhì)指標(biāo)的總解釋率達(dá)到了88.82%。ECON MN、CWED、PROX AM、ENN AM與 TN呈正相關(guān)，GYRATE AM和CIRCLE AM與TN呈負(fù)相關(guān)。ECON MN、PROX AM、ENN AM、GYRATE AM與TP呈正相關(guān)，CWED和CIRCLE AM與TP呈負(fù)相關(guān)（圖5）?？傮w來看CIRCLE AM和ENN AM對水質(zhì)指標(biāo)的影響較其他指數(shù)顯著，其他4個(gè)指數(shù)對水質(zhì)指標(biāo)的影響程度相接近。

2.3" "枯水期流域景觀格局與氮磷冗余分析

枯水期指數(shù)與氮磷數(shù)據(jù)的消除趨勢對應(yīng)分析結(jié)果顯示，排序軸中最長梯度值為0.9（小于3.0），采用RDA分析。

RDA分析結(jié)果顯示第一軸解釋率為42.5%，景觀指數(shù)對水質(zhì)指標(biāo)的總解釋率達(dá)到了89.57%。GYRATE AM、PROX AM與TN呈正相關(guān)，ENN AM、CWED、ECON MN與TN呈負(fù)相關(guān)；ENN AM、PROX AM、GYRATE AM與TP呈正相關(guān)，CWED、CIRCLE AM與TP呈負(fù)相關(guān)（圖6）?？傮w來看GYRATE AM、CWED對水質(zhì)指標(biāo)的影響較其他指數(shù)顯著，CIRCLE AM與水質(zhì)指標(biāo)的影響程度較弱。

2.4" "豐水期流域景觀配置與氮磷輸出關(guān)系

對豐水期水質(zhì)指標(biāo)與景觀配置進(jìn)行Pearson相關(guān)分析，結(jié)果如圖7、表3所示：森林的邊緣對比指數(shù)平均指數(shù)（ECON MN）與TN呈顯著正相關(guān)（Plt;0.05），農(nóng)田的邊緣對比指數(shù)平均值（ECON MN）和礦場的歐氏鄰近距離面積加權(quán)平均值（ENN AM）與TP呈顯著正相關(guān)（Plt;0.01），農(nóng)田和礦場的對比度加權(quán)邊緣密度（CWED）對TP呈顯著正相關(guān)（Plt;0.05），森林的鄰近指數(shù)面積加權(quán)平均值（PROX AM）與TP呈顯著正相關(guān)（Plt;0.05），森林的近圓形狀面積加權(quán)平均值（CIRCLE AM）與TP呈顯著負(fù)相關(guān)（Plt;0.05）。

2.5" "枯水期流域景觀配置與氮磷輸出關(guān)系

枯水期水質(zhì)指標(biāo)與景觀配置的Pearson相關(guān)分析結(jié)果如圖7、表3所示：農(nóng)田的邊緣對比指數(shù)平均值（ECON MN）與TP呈顯著正相關(guān)（Plt;0.01）（圖7），草地的對比度加權(quán)邊緣密度（CWED）對TN呈顯著負(fù)相關(guān)（Plt;0.05），森林的回旋半徑面積加權(quán)平均值（GYRATE AM）與TN呈顯著正向相關(guān)（Plt;0.05），水體的歐氏鄰近距離面積加權(quán)平均值與TP呈顯著正相關(guān)（Plt;0.05），礦場的歐氏鄰近距離面積加權(quán)平均值（ENN AM）與TP呈顯著正相關(guān)（Plt;0.01）。

3" "討論

3.1" "氮磷輸出差異性

豐水期氮磷主要來自3、4、6號子流域，3號子流域內(nèi)存在一個(gè)地下磷礦場和一個(gè)露天磷礦加工場，開采的磷礦以及加工廠內(nèi)大量堆積的磷礦石是該流域重要污染源。遠(yuǎn)安縣6、7月的降雨量是1、2月降雨量的20倍之多（圖4），降雨又是土壤養(yǎng)分流失的主要因素之一（Braskerud，2002；張亞麗等，2004），氮磷污染物便隨之流向下游。4號主要是3號子流域攜帶下來的污染物。6號子流域的城鎮(zhèn)面積占比較大還有少量礦場和較多工礦企業(yè)，在降雨情況下，城鎮(zhèn)不透水面流淌堆積的污染物失去地表植物截留和土壤吸附被沖刷至河流并流向下游（陳利頂?shù)龋?003；Fedorko et al，2005；羅璇等，2010；趙鵬等，2012）。

枯水期流域內(nèi)降雨量低于豐水期，同時(shí)氮磷輸出量也低于豐水期（圖3、4），3、6號子流域在枯水期依然受礦場的影響導(dǎo)致磷素輸出量高于其他流域。9、10號子流域的城鎮(zhèn)和農(nóng)田面積占比為流域最大，同時(shí)2個(gè)流域的森林面積占比是流域內(nèi)最少，枯水期降雨量少于豐水期數(shù)倍，1、2月因春節(jié)的原因城鎮(zhèn)居民區(qū)的活動量較其他月份高，且兩個(gè)流域中農(nóng)田面積占比較大，正值春耕時(shí)期，人們大量翻耕施加底肥增加了氮磷的流失（羅璇等，2010；李明濤等，2013），表明9、10號流域的氮磷污染大概率主要來自城鎮(zhèn)和農(nóng)田（羅璇等，2010；趙鵬等，2012；李明濤等，2013）。2個(gè)時(shí)期氮素輸出量差值在10 mg左右，枯水期的TN主要來自擁有果園面積最大的2號子流域，此結(jié)果與韓黎陽等（2014）在蘭陵溪流域的結(jié)論相似。2個(gè)時(shí)期的氮素差異大與氮素輸出有關(guān)，地表徑流中的氮素主要以水溶態(tài)輸出為主（曾立雄，2010），豐水期的TN、TP的濃度高于枯水期，且TN的流失量高于TP（韓黎陽等，2014）。

3.2" "流域景觀格局與氮磷輸出關(guān)系

本研究結(jié)果中顯示CIRCLE AM與2個(gè)時(shí)期的水質(zhì)結(jié)果都呈現(xiàn)負(fù)相關(guān)性，研究區(qū)域斑塊形狀越規(guī)則，與圓形偏離程度越小，其指數(shù)就越小。越接近圓形的斑塊，其邊界長度越短，與基質(zhì)的相互作用越小，且與相鄰的斑塊間的交流較少，有利于控制上游攜帶和本斑塊產(chǎn)生的污染物。同時(shí)圓形斑塊在自然資源保護(hù)方面有著最高效率。CWED能夠衡量一個(gè)斑塊與其他鄰近地區(qū)的對比程度，同時(shí)反映出邊緣密度和邊緣對比，本研究結(jié)果表明對TP都呈現(xiàn)負(fù)相關(guān)性。而ECON MN和CWED在豐水期和枯水期對TN的影響結(jié)果截然相反，豐水期降雨量充足使斑塊間的邊緣效應(yīng)得到更好地發(fā)揮，豐水期氮素輸出以地表徑流為主，枯水期降雨量成倍減少，氮素輸出途徑由地表徑流轉(zhuǎn)變?yōu)橐曰鳛橹鳎ɡ钗某龋?018），斑塊間的邊緣效應(yīng)受到限制，污染物遷移緩慢。PROX AM與豐水期的水質(zhì)指標(biāo)和枯水期的TP呈正相關(guān)，PROX AM量化了斑塊與其他鄰近斑塊間的空間環(huán)境關(guān)系，數(shù)值越小意味著斑塊的高度隔離，相較于其他指數(shù)對氮磷輸出的解釋度較低。ENN AM已被廣泛用于量化斑塊距離，是基于距離各自斑塊最近的2個(gè)斑塊中心之間的距離，斑塊間距離越近污染物的遷移速度也越快。ENN AM對豐水期的解釋度高于枯水期，基于TP主要以泥沙結(jié)合態(tài)經(jīng)地表徑流流失這一結(jié)論（曾立雄，2010），而TN在枯水期因降雨量減少，其輸出途徑由地表徑流轉(zhuǎn)變?yōu)榛餍问捷敵觯ɡ钗某龋?018）。GYRATE AM是對斑塊范圍的測量，斑塊越大回旋半徑越大，面積不變的同時(shí)斑塊的范圍變大，回旋半徑也越大，反映斑塊面積的同時(shí)也反映了斑塊的延伸程度。面積較大斑塊能夠保護(hù)水質(zhì)、發(fā)育河流網(wǎng)絡(luò)、緩沖能力強(qiáng)。研究結(jié)果表明，在枯水期對水質(zhì)指標(biāo)呈現(xiàn)正相關(guān)性，豐水期降雨量增加，地表徑流量和徑流總量都迅速增加，氮磷流失量也隨之增加（呂喚春，2002）。豐水期流域內(nèi)的水稻、玉米等主要農(nóng)作物正值開花結(jié)穗時(shí)期，經(jīng)過斑塊大面積土壤和作物根系的攔截與吸收，污染物的遷移受到了很好的控制。

3.3" "流域景觀配置與水質(zhì)指標(biāo)輸出關(guān)系

農(nóng)田的ECON MN和CWED值在2個(gè)時(shí)期只與枯水期TN呈負(fù)相關(guān)，枯水期大部分農(nóng)田處于休耕期，豐水期農(nóng)田大量施肥以及地面徑流增加，污染物攔截率下降，污染濃度較枯水期高。礦場的ECON MN和CWED也只與枯水期的TN指標(biāo)呈負(fù)相關(guān)，研究區(qū)內(nèi)的礦場主要分為地下磷礦、磷礦石加工廠。豐水期降雨充足，會有大量的徑流和滲漏對TN、TP進(jìn)行遷移，枯水期降雨量減少后地下礦井上層周邊的植被會對經(jīng)基流輸出的TN進(jìn)行攔截、吸附。草地的CWED與枯水期的水質(zhì)指標(biāo)和豐水期的TP呈負(fù)相關(guān)，草地依然是典型的污染物攔截土地利用類型與較多學(xué)者研究結(jié)果相似（李明濤等，2013；王晶萍等，2016）。森林的CIRCLE AM與水質(zhì)指標(biāo)呈負(fù)相關(guān)，森林對水質(zhì)的凈化功能在多處被證實(shí)（李明濤等，2013；焦勝等，2014；吉冬青等，2015；楊婭楠等，2016），較圓的斑塊形狀與相鄰和遠(yuǎn)距離基質(zhì)間的交流最少，斑塊的面積及形狀越規(guī)則時(shí)指數(shù)值越低，對水質(zhì)的凈化作用越突出（黃金良等，2011）。森林的GYRATE AM在枯水期與TN呈顯著正相關(guān)，GYRATE AM測量了斑塊面積的同時(shí)也測量了斑塊的延伸范圍，在森林的面積大小對水質(zhì)優(yōu)化有正向的作用時(shí)，也反映出斑塊的延伸程度對改善水質(zhì)有反向作用（黃金良等，2011）。礦場的ENN AM對TP呈顯著正相關(guān)，磷礦是環(huán)境中磷污染的主要因素，礦場周邊的土地利用類型會對TP的輸出進(jìn)行一定量的攔截和吸附，當(dāng)周圍都是同一土地利用類型時(shí)會加劇TP的聚集，造成更大的污染，所以磷礦之間距離的遠(yuǎn)近能較好地控制磷TP的聚集和輸出。森林的PROX AM對水質(zhì)指標(biāo)呈正相關(guān)，PROX AM指數(shù)越大斑塊的隔離度越低，森林斑塊與其他土地利用類型斑塊的交流就越豐富，物質(zhì)傳輸就更快捷，森林斑塊的隔離度對水質(zhì)的凈化作用有所減弱。水體在流域中呈上疏下密，整體地勢由山地到丘陵再到平原進(jìn)行一個(gè)過渡，上游河道水體分布較少，中下游水庫、小溪流及池塘分布逐漸增加，水體斑塊數(shù)量、密度增加（圖1）。城鎮(zhèn)依水而建，農(nóng)田、果園靠水灌溉，生活污水、工業(yè)廢水和灌溉用水通過河道、小溪、灌溉渠流向下游，污染物依附于水體向下輸出，以至于水體ENN AM與水質(zhì)指標(biāo)呈正相關(guān)。

4" "結(jié)論

本研究采用相關(guān)分析和冗余分析，對沮河流域景觀格局、景觀配置與氮、磷輸出關(guān)系進(jìn)行了不同時(shí)間和空間的研究，結(jié)論如下：

（1）磷礦是TP污染的最直接、影響力最大的因素。城鎮(zhèn)和其工礦企業(yè)的面積及數(shù)量是豐水期TP的主要貢獻(xiàn)類型?？菟赥P主要來源于農(nóng)田和城鎮(zhèn)， TN主要來源果園，豐水期的TN、TP濃度高于枯水期，且TN的流失量高于TP。

（2）從流域景觀格局來看，CIRCLE AM與水質(zhì)指標(biāo)呈負(fù)相關(guān)，證明斑塊的形狀越規(guī)則對流域氮磷污染物的截留越有效，2種邊緣指標(biāo)ECON MN和CWED在2個(gè)時(shí)期結(jié)合降雨情況都能很好反映不同的氮磷輸出情況。ENN AM和GYRATE AM都能從斑塊面積對水質(zhì)做出解釋，且豐水期解釋度高于枯水期。

（3）森林CIRCLE AM、GYRATE AM、PROX AM與水質(zhì)指標(biāo)表現(xiàn)出良好的相關(guān)性，對水質(zhì)的影響有較好的解釋，改善流域水質(zhì)可以從增加森林斑塊的面積，規(guī)范森林斑塊形狀，加強(qiáng)森林斑塊隔離性入手。同時(shí)也要增加草地斑塊的邊緣密度，加強(qiáng)對水質(zhì)的優(yōu)化作用。礦場的ENN AM對TP污染物的輸出表現(xiàn)最顯著，多個(gè)礦場之間的距離會起到一個(gè)連鎖反應(yīng)，對臨近斑塊及整體流域造成極大的影響。

參考文獻(xiàn)

陳利頂， 丘君， 張淑榮， 等， 2003. 復(fù)雜景觀中營養(yǎng)型非點(diǎn)源污染物時(shí)空變異特征分析[J]. 環(huán)境科學(xué)， 24（3）：85-90.

CHEN L D， QIU J， ZHANG S R， et al， 2003. Tempo-spatial variation of non-point source pollutants in a complex landscape[J]. Environmental Science， 24（3）：85-90.

陳崇德， 方修銀， 尚元， 等， 2008. 沮漳河流域水庫群聯(lián)合調(diào)度初探[J]. 中國防汛抗旱， 18（6）：47-50.

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

CAI H， HE Z W， AN Y L， et al， 2015. Response relationship between land use and water quality in Chishui river basin based on rs and GIS[J]. Resources and Environment in the Yangtze Basin， 24（2）：286-291.

黃金良， 李青生， 洪華生， 等， 2011. 九龍江流域土地利用/景觀格局-水質(zhì)的初步關(guān)聯(lián)分析[J]. 環(huán)境科學(xué)， 32（1）：64-72.

HUANG J L， LI Q S， HONG H S， et al， 2011. Preliminary study on linking land use amp; landscape pattern and water quality in the Jiulong River watershed[J]. Environmental Science， 32（1）：64-72.

黃永文， 陳安， 平川， 等， 2013. 宜昌市沮漳河流域當(dāng)陽段水環(huán)境現(xiàn)狀調(diào)查與評價(jià)[J]. 環(huán)境保護(hù)與循環(huán)經(jīng)濟(jì)， 33（8）：48-51.

韓黎陽， 黃志霖， 肖文發(fā)， 等， 2014. 三峽庫區(qū)蘭陵溪小流域土地利用及景觀格局對氮磷輸出的影響[J]. 環(huán)境科學(xué)， 35（3）：1091-1097.

HAN L Y， HUANG Z L， XIAO W F， et al， 2014. Effects of land use and landscape pattern on nitrogen and phosphorus exports in lanlingxi watershed of the Three Gorges Reservoir area， China[J]. Environmental Science， 35（3）：1091-1097.

黃益平， 王鵬， 徐啟渝， 等， 2021. 袁河流域土地利用方式對河流水體碳、氮、磷的影響[J]. 環(huán)境科學(xué)研究， 34（9）：2132-2142.

HUANG Y P， WANG P， XU Q Y， et al， 2021. Influence of land use on carbon， nitrogen and phosphorus in water of Yuan River Basin[J]. Research of Environmental Sciences， 34（9）：2132-2142.

焦勝， 楊娜， 彭楷， 等， 2014. 溈水流域土地景觀格局對河流水質(zhì)的影響[J]. 地理研究， 33（12）：2263-2274.

JIAO S， YANG N， PENG K， et al， 2014. The effects of land-use and landscape pattern on water quality in Weihe river watershed[J]. Geographical Research， 33（12）：2263-2274.

吉冬青， 文雅， 魏建兵， 等， 2015. 流溪河流域景觀空間特征與河流水質(zhì)的關(guān)聯(lián)分析[J]. 生態(tài)學(xué)報(bào)， 35（2）：246-253.

JI D Q， WEN Y， WEI J B， et al， 2015. Relationships between landscape spatial characteristics and surface water quality in the Liu Xi River watershed[J]. Acta Ecologica Sinica， 35（2）：246-253.

呂喚春， 2002. 千島湖流域農(nóng)業(yè)非點(diǎn)源污染及其生態(tài)郊應(yīng)的研究[D]. 杭州：浙江大學(xué).

羅璇， 史志華， 尹煒， 等， 2010. 小流域土地利用結(jié)構(gòu)對氮素輸出的影響[J].環(huán)境科學(xué)， 31（1）：58-62.

LUO X， SHIZ H， YIN W， et al， 2010. Effect of land use structure on nitrogen output in small watershed [J]. Environmental Science， 31（1）：58-62.

李明濤， 王曉燕， 劉文竹， 2013. 潮河流域景觀格局與非點(diǎn)源污染負(fù)荷關(guān)系研究[J]. 環(huán)境科學(xué)學(xué)報(bào)， 33（8）：2296-2306.

LI M T， WANG X Y， LIU W Z， 2013. Relationship between landscape pattern and non-point source pollution loads in the Chaohe River Watershed[J]. Acta Scientiae Circumstantiae， 33（8）：2296-2306.

李文超， 雷秋良， 翟麗梅， 等， 2018. 流域氮素主要輸出途徑及變化特征[J]. 環(huán)境科學(xué)， 39（12）：5375-5382.

LI W C， LEI Q L， ZHAI L M， et al， 2018. Main nitrogen output routes and change characteristics in watershed [J]. Environmental Science， 39（12）：5375-5382.

呂樂婷， 高曉琴， 劉琦， 等， 2021. 東江流域景觀格局對氮、磷輸出的影響[J]. 生態(tài)學(xué)報(bào)， 41（5）：1758-1765.

Lü L T， GAO X Q， LIU Q， et al， 2021. Influence of landscape pattern on nitrogen and phosphorus output in the Dongjiang River Basin[J]. Acta Ecologica Sinica， 41（5）：1758-1765.

彭勃， 付永勝， 趙君鳳， 等， 2019. 小流域土地利用景觀格局對水質(zhì)的影響[J]. 生態(tài)科學(xué)， 38（3）：90-99.

PENG B， FU Y S， ZHAO J F， et al， 2019. Effect of land use landscape pattern on water quality in small watershed[J]. Ecological Science， 38（3）：90-99.

祁蘭蘭， 王金亮， 葉輝， 等， 2021. 滇中“三湖流域” 土地利用景觀格局與水質(zhì)變化關(guān)系研究[J]. 水土保持研究， 28（6）：199-208.

QI L L， WANG J L， YE H， et al， 2021. Study on the land use landscape pattern and water quality change in the ‘three lakes’ basin in central Yunnan[J]. Research of Soil and Water Conservation， 28（6）：199-208.

沈濤， 劉良云， 馬金峰， 等， 2007. 基于L-THIA模型的密云水庫地區(qū)非點(diǎn)源污染空間分布特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)， 23（5）：62-68.

SHEN T， LIU L Y， MA J F， et al， 2007. Spatial distribution characteristics of non-point pollution of Miyun reservoir areas based on L-THIA model[J]. Transactions of the Chinese Society of Agricultural Engineering， 23（5）：62-68.

宋夢廷， 劉睿， 陳芳， 2022. 流域生態(tài)富營養(yǎng)化水體中氮磷污染修復(fù)植物篩選及效能研究進(jìn)展[J]. 環(huán)境保護(hù)與循環(huán)經(jīng)濟(jì)， 42（5）：14-20.

楊婭楠， 王金亮， 陳光杰， 等， 2016. 撫仙湖流域土地利用格局與水質(zhì)變化關(guān)系[J]. 國土資源遙感， 28（1）：159-165.

YANG Y N， WANG J L， CHEN G J， et al， 2016. Relationship between land use pattern and water quality change in Fuxian Lake basin [J]. Land and Resources Remote Sensing， 28（1）：159-165.

袁軼男， 金云峰， 聶曉嘉， 等， 2020. 基于Fragstats4的上海市城市景觀格局指數(shù)動態(tài)研究[J]. 山東農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版）， 51（6）：1157-1162.

YUAN Y N， JIN Y F， NIE X J， et al， 2020. Study on dynamic urban landscape pattern indexes in Shanghai City based on Fragstats4[J]. Journal of Shandong Agricultural University （Natural Science Edition）， 51（6）：1157-1162.

姚靜， 2022. 渭河流域景觀格局變化對徑流的影響研究[D]. 鄭州：鄭州大學(xué).

王晶萍， 李兆富， 劉紅玉， 等， 2016. 中田河流域景觀異質(zhì)性對水體總氮濃度影響研究[J]. 環(huán)境科學(xué)， 37（2）：527-533.

WANG J P， LI Z F， LIU H Y， et al， 2016. Influence of landscape heterogeneity on total nitrogen concentration in Zhongtian River watershed[J]. Environmental Science， 37（2）：527-533.

張亞麗， 張興昌， 邵明安， 等， 2004. 降雨強(qiáng)度對黃土坡面礦質(zhì)氮素流失的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)， 20（3）：55-58.

ZHANG Y L， ZHANG X C， SHAO M A， et al， 2004. Impact of rainfall intensity on soil mineral nitrogen loss by runoff on loess slope[J]. Transactions of the Chinese Society of Agricultural Engineering， 20（3）：55-58.

曾立雄， 2010. 三峽庫區(qū)蘭陵溪小流域養(yǎng)分的分布、遷移與控制研究[D]. 北京：中國林業(yè)科學(xué)研究院.

趙鵬， 夏北成， 秦建橋， 等， 2012. 流域景觀格局與河流水質(zhì)的多變量相關(guān)分析[J]. 生態(tài)學(xué)報(bào)， 32（8）：2331-2341.

ZHAO P， XIA B C， QIN J Q， et al， 2012. Multivariate correlation analysis between landscape pattern and water quality[J]. Acta Ecologica Sinica， 32（8）：2331-2341.

CHEN M， CHEN J， SUN F， 2008. Agricultural phosphorus flow and its environmental impacts in China[J]. Science of the Total Environment， 405（1-3）：140-152.

CAO D， CAO W Z， FANG J， et al， 2014. Nitrogen and phosphorus losses from agricultural systems in China： a meta-analysis[J]. Marine Pollution Bulletin， 85（2）：727-732.

BRASKERUD B C， 2002. Factors affecting nitrogen retention in small constructed wetlands treating agricultural non-point source pollution[J]. Ecological Engineering， 18（3）：351-370.

BAHUKHANDI D K， 2012. Impact assessment of urbanization and industrialization on surface and groundwater quality[J]. Journal of Engineering Design amp; Technology， 1（1）：11-22.

CAREY R O， MIGLIACCIO K W， LI Y C， et al， 2011. Land use disturbance indicators and water quality variability in the Biscayne Bay Watershed， Florida[J]. Ecological Indicators， 11（5）：1093-1104.

FEDORKO E J， PONTIUS R G JR， ALDRICH S P， et al， 2004. Spatial distribution of land type in regression models of pollutant loading[J]. The Biological Bulletin， 207（2）：173.

LIU J， LIU X L， WANG Y， et al， 2020. Landscape pattern at the class level regulates the stream water nitrogen and phosphorus levels in a Chinese subtropical agricultural catchment[J]. Agriculture， Ecosystems amp; Environment， 295：106897.

LIU H Y， MENG C， WANG Y， et al， 2022. Multi-spatial scale effects of multidimensional landscape pattern on stream water nitrogen pollution in a subtropical agricultural watershed[J]. Journal of Environmental Management， 321：115962.

SHEN Z Y， HOU X S， LI W， et al， 2015. Impact of landscape pattern at multiple spatial scales on water quality： a case study in a typical urbanised watershed in China[J]. Ecological Indicators， 48：417-427.

SAVAGE C， LEAVITT P R， ELMGREN R， 2010. Effects of land use， urbanization， and climate variability on coastal eutrophication in the Baltic Sea[J]. Limnology and Oceanography， 55（3）：1033-1046.

SARKER B， KEYA K N， MAHIR F I， et al， 2021. Surface and ground water pollution： causes and effects of urbanization and industrialization in south Asia[J]. Scientific Review， （73）：32-41.

TURNER M， 1989. Landscape ecology： the effect of pattern on process[J]. Annual Review of Ecology and Systematics， 20：171-197.

TENG Y G， YANG J， ZUO R， et al， 2011. Impact of urbanization and industrialization upon surface water quality： a pilot study of Panzhihua mining town[J]. Journal of Earth Science， 22（5）：658-668.

XU Q Y， WANG P， SHU W， et al， 2021. Influence of landscape structures on river water quality at multiple spatial scales： a case study of the Yuan river watershed， China[J]. Ecological Indicators， 121：107226.

ZHOU T， WU J G， PENG S L， 2012. Assessing the effects of landscape pattern on river water quality at multiple scales： a case study of the Dongjiang River watershed， China[J]. Ecological Indicators， 23：166-175.

（責(zé)任編輯" "鄭金秀" "崔莎莎）

Influence of Landscape on Nitrogen and Phosphorus Exports from Water Bodies in Juhe River Basin

XIANG Qian， JIANG Shunyao， XIONG Qinxue， ZHU Jianqiang

（College of Agriculture， Changjiang University， Jingzhou" "434025， P.R. China）

Abstract：The impact of landscapes on river water quality has become a research hotspot in river ecology both domestically and internationally. In this study， the Juhe River basin was selected as a case study， and we explored the impact of land use and landscape on nitrogen and phosphorus outputs in the basin， aiming to provide data support for local water resource management， ecosystem protection， land use structure improvement and landscape optimization. The spatial distribution of land uses and sub-basin divisions within Juhe River basin was interpreted based on remote sensing images from 2022， and six land use types were delineated， including farmland， surface water， forest， grass land， cities and towns， orchards and mines. A water quality survey in the Juhe River basin was conducted in February （dry season） and July （wet season） of 2023， at 10 river outlets representing 10 sub-basins. The spatiotemporal relationship between the landscape of sub-basins and the nitrogen and phosphorus output at river outlets was investigated using Pearson correlation and redundancy analysis. Results show that： （1） Phosphate mining was the most important factor influencing phosphorus pollution. During the dry season， the primary sources of nitrogen and phosphorus were orchards， farmlands， and urban areas. （2） The area-weighted average of the circularity index had a positive effect on the retention of nitrogen and phosphorus in the basin. （3） During the dry season， key factors affecting nitrogen and phosphorus concentrations in Juhe River included the area-weighted average of the Euclidean proximity distance from mining sites and water bodies， the area-weighted average of forest patch radius， the area-weighted average of grassland contrast， and the area-weighted average of farmland edge contrast index. During the wet season， the primary driving factors included the area-weighted average of farmland and forest edge contrast index， the area-weighted average of the farmland and mining contrast edge density， the area-weighted average of the forest circularity index， the area-weighted average of Euclidean proximity distance from mining sites， and the area-weighted average of the forest proximity index. This study provides valuable insights for controlling nitrogen and phosphorus concentrations in Juhe River.

Key words：landscape; land use; nitrogen and phosphorus outputs; Juhe River basin

基金項(xiàng)目：國家自然基金-區(qū)域創(chuàng)新發(fā)展聯(lián)合基金項(xiàng)目（U21A2039）。

作者簡介：向前，2000年生，男，碩士研究生，專業(yè)方向?yàn)檗r(nóng)業(yè)遙感。E-mail：2297176415@qq.com

通信作者：蔣舜堯。E-mail：759130021@qq.com