1 引言(Introduction)

流域水质是多尺度环境因子的综合反映, 污染物进入河流水体前在地表迁移、转化过程中会受到地理环境(植被、土壤)、气候环境和人类活动(城市化、工业化等)等因子的共同作用(Wu et al., 2018).这些因子反映了景观格局在空间上的组成和结构, 因此, 景观格局影响着流域内生态和水文过程, 是造成河流水质空间分异的重要原因(井云清等, 2017; Mainali et al., 2018).

随着城市化发展和人口增长, 景观格局从河岸带到流域再到全球等不同尺度上都随之表现出不同的空间分布特征(Ding et al., 2016), 而水环境对景观格局变化的响应尤为敏感(Finlayson et al., 2018).已有国内外较多研究表明, 景观格局对水环境的效应从局部小尺度到流域等不同空间尺度上不尽相同(Ding et al., 2016; De Mello et al., 2018).这种由空间尺度差异引起的景观与水质关系的不同也是目前研究面临的一个挑战(李艳利等, 2015).但之前的研究结果对具体识别流域中景观格局与河流水质变化相关性最强的空间尺度具有较大的争议(Shi et al., 2017).Meneses等(2015)和Ding等(2016)分别在葡萄牙和中国广东某流域的研究中发现, 景观格局在流域尺度上对水质变化的影响程度大于河岸带缓冲区尺度.而Hille等(2018)和Xu等(2019)分别在丹麦515个河流流域和中国贵州乌江河流域的研究结果表明, 河岸带缓冲区尺度可以更好地保护河流水质.这些不一致的研究结果主要由于影响水质的每个流域所具备的独特性所导致(Sun et al., 2014).但至今的研究中仍缺乏各种河岸带缓冲区宽度变化对水质动态影响的定量分析.河岸带直接作用于水体, 是敏感的生态交错区, 同时又浓缩了整个流域生态环境的特征(Zhang et al., 2018; Ruan et al., 2019), 且单一空间尺度往往无法全面解释景观格局与河流水质之间的关系.因此, 划分不同空间尺度河岸带缓冲区, 识别流域内景观格局与河流水质变化相关性最强的空间尺度, 对于流域水环境管理和景观生态学在流域水环境中的应用具有十分重要的意义.

随着河岸带空间尺度的变化, 景观组成和景观结构对水质变化的影响机制和程度也具有差异性.Clément等(2017)在加拿大约600个河流流域案例的研究中表明, 水质对不同景观组成类型面积变化(如城市用地、农用地、林地等)的响应程度存在差异性.Xie等(2018)在流经中国香港与广东的深圳河流域的研究中发现, 景观组成变量与溶解氧、总磷、悬浮物等水质参数的相关性随空间尺度变化而发生改变.Lee等(2009)在韩国某流域中的研究结果表明, 城市区域斑块破碎化程度越低对水环境质量改善效果越明显, 而林地分散分布的区域其水质退化的现象非常显著.Zhang等(2019)在中国三峡水库流域的研究则反映了景观结构中斑块密度和最大斑块指数对总氮、硝态氮、总磷等水质参数的影响作用相反.通过前人的研究发现, 可能由于景观类型在空间上的组成和结构不同, 其对水体中营养物浓度的效应和影响机制存在一定的差异性.然而目前对景观组成、结构与水质变化之间具体关系的研究还多为定性研究, 多数仅限于不同景观类型对水质的影响或聚焦在较大的空间尺度上, 而将景观格局的空间组成和结构作用量化的研究相对较少; 同时也缺乏在不同空间尺度上识别影响流域水质变化的主要景观因子及各景观格局指数产生的具体贡献率的研究(Shi et al., 2017).因此, 使用景观格局指数方法区分景观组成和结构变量, 定量分析其与水质变化的关系, 对于明确区域土地利用/覆被类型、景观要素空间最佳分布及其最优配置关系有着重要的参考价值.

随着国家"长江大保护"战略的提出, 汉江作为长江的第一大支流, 其水环境保护意义重大.襄阳市作为汉江中游的最大城市, 景观格局变化剧烈, 人口密集, 地表水污染趋势加重, 对汉江水环境的压力仍在持续增加.同时, 汉江作为襄阳市城区的唯一水源地和纳污水体, 其水环境状况将在很大程度上受城区发展与扩张的影响.因此, 研究汉江水质变化与襄阳市城区景观格局的响应关系对于汉江流域水环境保护和长江经济带的可持续发展具有现实意义.因此, 本文以汉江流域襄阳城区段为研究对象, 在河岸带尺度上, 用包含景观组成和结构信息的景观格局指数表征景观格局状况, 通过遥感(RS)、GIS技术、冗余分析等模型与方法对景观格局与河流水质变化之间的关系进行研究.主要包括3个方面:1识别景观格局对流域水质变化影响最大的空间尺度; 2区分景观组成和景观结构对流域水质的影响程度大小; 3定量分析不同缓冲区宽度上景观格局对水质变化影响效应的差异性.

2 数据来源与研究方法(Data sources and methods) 2.1 研究区域

汉江是长江最大的支流, 长1577 km, 流域面积约1.59×10⁵ km².襄阳市位于汉江中游(110°45'~113°43'E, 31°14'~31°37'N), 全市总面积约为19774 km², 为北亚热带季风气候, 年平均气温约为14 ℃, 年平均降水量为800~1200 mm(图 1).自2004年以来, 襄阳市城市扩张迅速, 城市化率从37.8%上升到2017年的59.7%.汉江流域襄阳城区段总面积约3557 km², 包括从白家湾到余家湖的监测断面, 其主河道长度约49 km.小清河和唐白河是汉江流域襄阳城区段内的两条支流, 分别于清河口和张湾汇入汉江.但据湖北省和襄阳市水资源公报, 该区域内人均水资源量仅约为947 m³, 远远低于湖北省人均水资源量的2052 m³, 水资源状况非常严峻, 同时该段城区河流是襄阳城区内约220万居民重要饮用水源地, 并具有灌溉、养殖和运输功能.

2.2 水质数据监测与分析

根据汉江流域襄阳城区段干流和支流的特点, 在汉江流入、流出城区的进口和出口, 以及支流汇入干流的汇入点等区域设置8个监测断面(图 1), 其中4个位于干流(S1~S4), 3个位于支流小清河(S5~S7), 1个位于支流唐白河(S8).选取水温(T)、pH、溶解氧(DO)、高锰酸盐指数(COD_Mn)、总氮(TN)、总磷(TP)、氨氮(NH₄⁺-N)等7项指标来分析河流水环境质量.其中, 2017年1—6月的水质监测数据来源于襄阳市环境监测站在每个断面的逐月监测分析, 2017年7—12月的水质监测数据来源于研究团队每月一次的现场采样和实验室分析; 水质采样与分析测试采用国家标准分析方法(赵晓等, 2018), 水质评价标准依据国家《地表水环境质量标准》(GB 3838—2002).

2.3 缓冲区空间尺度选取

结合汉江流域襄阳城区段河岸带的实际情况和研究目的, 应用ArcGIS软件中的空间分析功能, 分别对汉江的8个水质监测断面建立8个空间梯度的河岸带缓冲区, 以监测断面为基点, 平行河道分别向河流上游延伸至研究区域边界, 以垂直河道100、200、300、400、500、1000、1500、2000 m的宽度为规律在河流两岸建立缓冲区, 由此获得64个作用区.每个缓冲区范围内的景观类型均在ArcGIS软件界面划分提取、叠加和分类统计.同时, 为便于景观格局指数的统计分析及从研究精度和复杂度上考虑, 在各个尺度缓冲区上矢量数据转换成栅格数据时, 设置栅格大小为10 m×10 m, 河流缓冲区划分见图 1.

2.4 景观格局指数选取

研究区域2017年的景观格局数据是以2009年第二次全国土地调查数据为基础, 并以2017年6月Landsat TM影像为数据源, 通过监督分类解译进行更新.分类系统参考《土地利用现状分类》(GB/T21010—2007)中方法, 并结合当地的实际情况, 将研究区域景观类型划分为8种:城市用地(URB)、农用地(AGR)、林地(FOR)、草地(GRA)、道路(ROA)、水域(WAT)、园地(GAR)、未利用地(UNU).

景观格局指数能够高度浓缩景观格局信息, 并反映其空间组成和结构等方面特征, 其中, 景观组成变化表现为不同景观组成类型数量或均匀度的变化, 而景观结构变化则表现为空间格局复杂程度的变化(赵鹏等, 2012; Shen et al., 2015).根据研究目的, 结合汉江流域襄阳城区段的土地利用及社会经济发展的实际情况, 从表征景观优势度、破碎度和多样性等景观意义出发, 优先选择经文献研究可能对水质变化有较大影响的景观组成和景观结构变量.其中, 景观组成变量包括斑块面积比(PLAND)和均匀度指数(SHEI), 景观结构变量包括斑块密度(PD)、最大斑块指数(LPI)、斑块数量(ND)、聚集度(AI)、香农多样性指数(SHDI)、景观形状指数(LSI)(Palmer, 2004; 刘怡娜等, 2019).同时, 在类型水平上分析PD、LPI、PLAND、ND、AI等5个景观格局指数, 在景观水平上分析SHDI、SHEI、LSI等3个景观格局指数, 所有景观格局指数均由基于栅格数据的Fragstats 4.2软件计算获得.

2.5 数据统计与分析

使用SPSS 22.0统计各监测断面每月的水质指标数据, 并将每个断面所有月份数据做平均分析, 得到各断面的水质指标浓度年均值.以影响水质变化的景观类型和景观格局指数等景观因子为研究对象, 采用冗余分析(RDA)研究景观格局指数与水质参数年均值之间的关系.RDA是一种约束性直接梯度排序方法, 可用二维排序图直观地定量分析景观格局指数与水质参数之间的关系.选择RDA之前, 需采用CANOCO 4.5软件对各监测断面的水质参数进行降趋势对应分析(DCA), 结果显示, 4个坐标轴在所有缓冲区上的最长梯度值都小于3.因此, 适用于线性RDA模型来研究多空间尺度下景观格局指数与水质变化的关系(Shen et al., 2015), 该模型可进一步定量揭示单一景观变量对湖泊水质变化的解释贡献率.在RDA结果图中, 当两个变量的箭头方向相同时, 表示存在正相关关系, 反之亦然; 两个箭头之间的夹角角度与它们的相关程度成反比(Ding et al., 2016).

3 结果(Results) 3.1 不同景观类型面积的空间变化特征

不同宽度河岸缓冲区上景观类型面积比例的平均值分析结果见图 2.各景观类型在100~2000 m的河岸带缓冲区内都呈现了明显的空间变化特征.农业用地、水域、林地、城市用地这4种景观类型的面积比例在100~2000 m的空间尺度内占据了约95%以上.其中, 农业用地占比随缓冲区宽度增加呈逐步上升的空间变化特征, 其比例从100 m缓冲区的54.3%逐渐上升到2000 m缓冲区的64.5%, 并且占比最高, 是流域内优势景观类型.城市用地主要分布于汉江襄阳城区段的中游和下游之间, 其面积比例总体上随缓冲区宽度增加而上升, 并在1000 m缓冲区上达到最大值16.92%.林业用地主要分布于汉江干流的右侧, 其面积比例随缓冲区宽度增加而下降, 从100 m缓冲区的9.6%下降到2000 m缓冲区的6.9%.其他景观类型中, 水域、草地的面积比例在空间尺度上呈现随缓冲区宽度增加而逐步下降的特征, 分别从100 m缓冲区的28%、0.83%下降到2000 m缓冲区的8%、0.51%;园地、未利用地、道路面积比例呈现相反的逐步增加趋势, 分别从100 m缓冲区的0.30%、0.62%、0.58%上升到2000 m缓冲区的0.66%、1.06%、1.27%.以上各景观类型的变化趋势表明林业和水域等自然生态用地面积在空间尺度上逐渐减少, 而造成城市空气环境、水环境污染压力增加的工业、交通等城市用地和带来面源污染的农业用地的总面积逐渐增长, 综合反映了研究区域内城市河流水环境面临的巨大压力.

3.2 水质参数空间变化特征

通过对S1~S8样点断面水质参数年均浓度的对比分析(表 1), 结果表明, DO、COD_Mn、TP、NH₄⁺-N均不超过地表水Ⅲ类.但支流小清河、唐白河汇入干流的S7、S8断面中COD_Mn、TP、NH₄⁺-N、TN的浓度均高于汉江干流, 其中, TN浓度最高, 均达到了地表水Ⅴ类.反映了汉江干流的水质在流域内较好, 支流对干流的水质退化有一定的影响.同时, 汉江干流出城区域的S4断面水质浓度高于进城区域的S1断面, 表明随着河流的流向, 从上游到中游、下游, 污染物不断地富集于汉江干流水体中, 造成河流水质污染逐步加重, 尤其是TP、NH₄⁺-N浓度增加幅度较大, 反映了城市区域的点源(包括工业废水、生活污水等)和面源污染对汉江的水质产生了较大的影响.

3.3 景观格局对水质影响的空间尺度性特征

景观格局指数与水质指标的RDA分析结果如表 2所示.结果表明, 河岸带缓冲区景观格局指数与水质指标有着不同程度的空间相关性, 揭示了景观格局对水质变化效应的空间尺度性特征.表 2显示了在100 ~2000 m宽度缓冲区上最主要的景观格局指数对水质影响的解释率变化情况.在所有空间尺度内景观格局指数对水质空间变化的解释率均达到50%以上, 并且前两个轴占所有轴总解释量的92%以上, 数值较大, 能够很好地反映分析结果.总体上, 景观格局指数对水质的解释量随缓冲区宽度的增加呈现先升后降的特征.在100 m缓冲区尺度上, 景观格局指数对水质变化的累计解释量为70.5%;当缓冲区宽度增加到300 m时, 总解释率增加到78.5%, 达到了所有河岸带缓冲区中的最大值; 但当缓冲区宽度从400 m增加到2000 m时, 景观格局对水质的解释能力出现下降趋势, 在2000 m缓冲区中总解释量最低, 仅为50.7%.因此, 在本研究中的河流流域内, 300 m宽度的河岸缓冲区被确定为流域内景观格局在不同空间尺度上对水质变化解释率最大的区域, 也是研究区域河流水质管理的有效缓冲区.

3.4 景观组成和结构对水质变化的影响程度

RDA分析结果可以很好地识别影响水质空间分异的主要景观因子, 以及从整体上区分景观组成变量和景观结构变量对水质空间变化影响的大小.表 2清晰地表明研究区内不同景观类型对水质的影响存在差异性.景观结构变量PD_FOR在100 m、1000~2000 m范围内是对河流水质变化影响最明显的变量, 解释能力达到了46.5%~56.5%;PD_URB在200~300 m范围内是对河流水质影响最显著的变量, 解释能力分别达到了52.7%和57.2%;PD_GRA和LPI_FOR分别在400 m和500 m缓冲区上是对河流水质变化影响最明显的变量, 对水质的解释量分别达到了56.2%和52.9%.而景观组成变量仅PLAND_GAR在400 m尺度上对水质有一定的影响, 表明研究区域内景观组成变量对水质的影响程度小于景观结构变量.同时上述结果也表明, 在河流流域内随着缓冲区宽度的增加, 林业用地对流域水质的影响比较明显, 而城市用地的影响作用在减小.反映了城市用地在较小尺度(小于300 m)的河岸带缓冲区宽度内对水质影响最显著, 而林业用地在尺度相对较大的缓冲区宽度内(大于500 m)是影响水质的主要因子.同时以上结果也表明, 景观格局在类型水平上对河流水质的影响程度更为明显, 在景观水平上仅LSI在300 m宽度缓冲区上对水质有一定解释量(21.3%), 与水质参数的相关关系较小.

3.5 不同缓冲区宽度景观格局对水质变化的影响

图 3揭示了不同缓冲区宽度上景观格局指数与水质参数之间的相关关系.PD_URB和PD_FOR在100~300 m、500~2000 m内基本与所有水质参数(COD_Mn、TN、NH₄⁺-N、TP)呈显著的负相关关系, 与DO为正相关关系, 表明林业用地和建设用地的斑块密度越小, 对水质参数的上升有显著作用.LPI_GRA在200 m缓冲区宽度上与COD_Mn呈正相关, 与DO呈显著负相关, 但PD_GRA在400 m尺度上却与COD_Mn呈负相关关系, 揭示了斑块面积越大的草地对水质的负面影响越大, 但斑块密度越大却对水质的净化作用更显著.PLAND_GAR和LPI_GAR在400 m和500 m尺度上与TN、NH₄⁺-N、TP呈正相关, 表明园地的面积比例和最大斑块面积与水质的恶化呈正相关.景观水平上的景观格局指数仅LSI在300 m尺度上对水质有一定解释量(21.3%), 但与水质参数的相关关系较小.

4 讨论(Discussion) 4.1 景观格局对水质效应的最有效空间尺度识别

生态系统及其格局在空间上多尺度存在, 景观格局和地表水环境质量之间的关系受不同空间尺度上各种因子的影响, 具有空间异质性(Margriter et al., 2014; Zhang et al., 2019).虽然河岸带缓冲区对水质影响的最大宽度范围仍然存在争议, 但之前研究表明每一条河流均有一个临界缓冲区, 即存在流域内景观格局和水质变化之间相关关系最大的河岸带缓冲区(Sun et al., 2014).Tran等(2010)在纽约哈德逊河流域的研究中发现, 在200 m空间尺度上景观格局可以最大程度的解释水质变化, 但随观测空间尺度的增加, 景观格局对水质变化的影响呈逐渐下降趋势.王小平等(2017)通过对中国新疆艾比湖流域的案例研究, 发现景观格局在4 km缓冲区上对流域水质变化的解释能力最强.Shen等(2015)发现景观格局指数在300 m宽度缓冲区可以达到所有空间尺度上对水质变异的最大解释率.Dai等(2017)在中国抚仙湖流域中发现200 m河岸带是保护地表水环境质量的最适宜空间尺度.在本研究区域内, 景观格局(包括景观组成和结构)在不同河岸带缓冲区尺度上对河流水质指标的解释贡献率差异较大, 但300 m河岸带宽度对水质变化的影响程度最大.因此, 该宽度的河岸缓冲区是作为景观格局与水质关联分析的最有效空间尺度, 是河岸缓冲区对河流水环境影响最为直接和迅速的区域.该空间宽度的识别为确定河流水质改善的有效范围, 实现研究区域内景观格局的合理配置及土地利用的合理优化与管理提供了重要的科学依据.

4.2 景观组成和结构对水质影响程度的差异性分析

不同景观类型的组成(面积等)和景观结构(斑块密度、聚集度、多样性指数等)的变化对进入水体污染物的吸收和降解过程存在差异性, 是影响河流水质的重要因素(Shi et al., 2017).之前的研究中发现, 不同流域面积内景观结构变量与景观组成变量对水质的影响程度不一致(Campagnaro et al., 2017), 在大于250 km²的流域中, 景观组成对水质的影响程度大于景观结构(Clément et al., 2017).在中国太湖流域的研究中发现景观结构与水质的相关关系更加显著(Xiao et al., 2016), 而在河南伊河流域却发现景观组成与水质的相关性更加明显(Bian et al., 2019).在本研究区域内水质参数与斑块密度、最大斑块指数、景观形状指数等景观结构变量有显著的相关关系, 但景观组成变量仅PLAND_GAR在400 m缓冲区上对水质有一定的影响, 解释率仅为18.2%, 影响程度明显小于景观结构变量PD_GRA.表明不同景观类型的密度、最大斑块面积、形状指数等对水质的影响作用比景观类型的面积比例更显著, 反映了景观结构变量在流域中对河流水质的影响程度在不同空间尺度中均大于景观组成变量.

4.3 不同空间尺度景观格局对水质的差异性效应分析

由于人类活动干扰及景观结构异质性对污染负荷影响程度的差异性, 使得景观类型和水质之间必然存在着关联性(夏品华等, 2016; Duarte et al., 2018).本研究分析结果表明, 城市用地、林业用地、园地的景观结构变化是影响河流水体中TN、COD_Mn、TP、NH₄⁺-N等污染因子浓度变化的重要因素.在之前较多研究中发现, 农业用地(如水田、旱地)与流域内的水质之间具有显著相关关系, 并得出农业用地是流域内水质面源污染的重要来源, 如化肥、农药的过量使用等, 特别是在城市化率不高的流域表现尤为显著(Lee et al., 2009; Yang et al., 2012).在农业用地比例最高的太湖流域研究中发现, 林地、城市用地、农业用地对进入水体的TN、TP、COD_Cr等污染物的输出能力依次上升(Li, 2006).但在同样以农业用地占主导类型的本研究区域的结果却相反.农业用地与河流水质之间没有发现明显的相关关系, 这可能是由于与城市区域的废水排放相比, 流域内农业活动产生的污染物相对较少.这一结果与之前的研究类似(Tu, 2011), 即当一个地区的城市化率超过农业用地占土地利用总面积的百分比时, 农业活动产生的污染不是主要污染来源, 特别是在快速城市化流域.流域内园地面积虽然占比极低, 分布较分散, 且在各缓冲区域内分布不均衡, 但在本研究中PLAND_GAR和LPI_GAR与大部分水质参数呈正相关关系, 揭示分布面积较大的园地对河流水质产生了不利影响, 是流经城市区域的河流流域内重要的面源来源之一(Wang et al., 2019).

城市用地一般对流域内的水质退化起着重要作用(Shen et al., 2015), 尤其是小尺度上的城市用地对水质的影响更为明显(Shi et al., 2017).研究区内城市建设用地面积占比不高, 但表征其景观破碎度的斑块密度在整个河岸带缓冲区中对水质产生了极大的影响, 并在200 m和300 m宽度缓冲区中对水质的影响最大.表明在汉江流域襄阳段的河岸带中存在大量不透水面的城市区域, 会导致水质急剧恶化, 同时由人类活动产生的大量工业废水和生活污水在分散的城市区域会加剧水质退化程度.反映了随着流域内人口、经济和城市化水平的提高, 在分散的城市区域内其污水收集和处理设施没有跟上城市和人口的扩张速度与规模, 从而对河流水质造成了严重的压力(Ding et al., 2016; Liu et al., 2018).

通常情况下, 林地和草地被认为是有污染物截留功能的"汇"景观, 能够减少面源污染对水质的负面影响(Ouyang et al., 2010; Xiao et al., 2016).但本研究区域内PD_GRA与LPI_GRA对水质变化产生了不同的影响, 可能是由于密集分布的自然草地对进入河流的污染物有较好的净化作用.而在靠近河流区域, 大面积的人工草地被用作城市绿化用地, 在其种植和养护过程中大量施用化肥和农药, 最终导致未被完全利用的污染物进入河流造成水质下降(Xiao et al., 2016).因此, 汉江流域襄阳段河岸带中的人工草地和自然草地因种植方式的差异, 其对水环境的影响作用相反.与此相反, 在全球不同流域的研究案例中发现林业用地与水质在相对较大的空间尺度上相关关系更为显著(Woli et al., 2004; Zhang et al., 2019); 且城市化流域内林地比例较低的区域其水体污染物浓度较高(Bernhardt et al., 2007; Ruan et al., 2019).在本研究中, 林业用地的斑块密度和最大斑块面积指数在100 m、500~2000 m上是影响流域河流水质退化的主要景观格局指数.说明林地景观的斑块面积和密度越大, 且林地整体连通性和聚集程度越高, 对水质污染物的固定和截留作用越显著.其原因主要是由于林地冠层下土壤层具有涵养水源、保持水土、滞留部分水体污染物的作用, 对地表径流具有良好的净化效果(郭玉静等, 2018).因此, 流域内大面积分布、聚集且连接良好的林业用地对进入河流的污染物具有更强的遏制作用, 以及在临水河岸带(小于100 m)中分布的林业用地是维持河流水质良好的重要因素.

综合以上分析得出, 为了防止流域内河流水环境污染趋势加重, 对河流水质改善应重点强调与保护300 m宽度河岸缓冲区范围; 在城市区域需加强污水收集与处理设施的建设, 以减少污染物排入水体; 同时在靠近水体区域, 需要降低园地和人工草地种植强度, 增加林地和自然草地的覆盖密度和面积, 从而达到对入河污染物截留的最佳效果.

5 结论(Conclusions)

1) 不同宽度缓冲区对河流水质指标的效应具有显著的空间尺度性特征, 300 m缓冲区是流域内景观格局对水质变化作用最强的河岸带宽度.

2) 景观结构变量在不同空间尺度上对河流水质指标的影响程度均大于景观组成变量, 同时景观格局在类型水平上对水质的影响比景观水平显著.

3) 城市用地在较小尺度上(不超过300 m)是对河流水质影响最大的景观类型, 林地在相对较大的空间尺度及临水区域作用最显著.

4) 城市区域分布越集中对水质退化影响越严重; 林地分布越密集、面积越大对流域内进入河流水体污染物的吸收和降解作用越强; 受流域内人类活动方式的影响, 农业用地与水质之间没有发现显著的相关性; 受施肥、农药等因素影响, 大面积分布的园地和人工草地对水体污染物浓度上升有一定的贡献, 而自然草地呈现相反的作用.