近年来中国城市化进程仍旧在不断推进, 预计在2030年中国城镇化水平将达到70.99%^[1-2]。城市化进程不但改变了原本自然或半自然下垫面的景观类型和布局, 而且改变了自然生态系统水文和地球化学循环过程^[3-6]。中国城市水环境污染形势严峻, 90%以上的城市河流受到不同程度污染, 自然过程已经基本丧失了对城市水环境质量的调控能力, 人类活动已经成为水环境及水质变化的主导要素, 城市水环境的保护和治理引起了广泛关注^[7]。人类活动通过改变景观结构和空间格局, 改变了水文过程和径流路径, 从而影响了污染物的发生位置、迁移路径和转化过程, 进而对河流水环境产生深刻影响^[8-9]。因此, 从景观生态学的角度出发, 探索景观内不同土地利用及其格局对河流水环境的影响, 可为通过有目的地调节人类活动来调控景观空间格局、改善其服务功能提供重要的理论依据^[10-13]。

大量研究表明城市景观格局对于河流水环境存在不同程度的影响^[14-17]。如Wang等^[18]在上海城区、郊区和农村不同城市化梯度上应用综合污染指标比较了地表水质的变化, 表明地表水质的空间格局取决于城市化水平；Tu^[19]在不同城市化梯度上选择了129个子流域进行分析表明土地利用对水质的解释能力随着城市化程度的不同而变化, 如农业用地和居住用地比例与污染物浓度在城市化较低的区域呈正相关, 而在高度城市化区域却呈负相关；Shen等^[20]发现景观组成指标可解释46.9%的水质变化, 而如果加入景观构型指标, 则可以解释水质变化的60.3%, 同时类型水平上的景观指数对水质的预测比景观水平上的指数更为有效；Hasenmuellerd等^[21]在农村-城市的城乡梯度带上发现洪水事件中城市河流的水质特征变异程度更大, 可达农村河流水质变异幅度的5倍, 对环境质量有着较大的负作用；Deng^[22]发现水质与河网的结构和连通性具有显著的相关性, 与以城市用地为主的区域与耕地为主的区域相比, 在很大程度上减弱了这种相关性, 耕地与城市用地的季节性作用则相反, 认为河网结构和连通性有利于河流水质改善。

总体上, 不同的景观指数对不同的水质指标在不同的缓冲区研究尺度上有着不同的解释能力。但是, 由于景观的复杂性和异质性, 一些共性规律亟待深入发掘。景观格局和水质关系的研究一直是景观生态学的一个热点和难点^[23-25]。不同景观类型和景观格局的变化对水质的影响存在差异性^[26]。

本文以京杭运河的起点杭州段及其相通相连河网水系区域作为研究区, 河流总趋势为南北向^[27], 南依钱塘江, 北至杭州德清边界, 以水质监测点为中心建立不同宽度的缓冲区, 通过计算景观格局指数, 并与不同的水质指标进行相关性分析、冗余分析等, 定量分析不同空间尺度下景观格局与河流水质的关系, 以探讨流域景观空间分布格局对水质状况的影响程度, 为建立更为合理的景观结构, 改善运河水系河流水环境提供依据。

1 研究区概况

杭州作为长江三角洲中心城市之一, 总面积16850 km², 建成区面积559.2 km², 常住人口1036万人, 城镇人口813.26万人, 城镇化率78.5%^[28]。京杭大运河是世界上里程最长、工程最大的运河, 也是最古老的运河之一, 南起杭州(余杭), 北至北京(涿郡), 贯通海河、黄河、淮河、长江和钱塘江五大水系, 全长约1794 km^[29]。杭州是京杭大运河的南端起点, 京杭大运河从钱塘江配水自南往北穿越杭州城区, 与之相通的河道纵横, 水网发达, 水文过程复杂多变^[27]。运河杭州段(图 1)是杭州主城区水位最低的地表河流, 两侧支流众多, 是市区各类河水、地表径流和众多支流的主要受纳水体, 城区与运河相连的河道多达33条(除钱塘江外)^[30], 每年进入运河的水量达到4.5×10⁸ m³, 对运河水质产生了巨大的负面影响^[31-32]。另外, 运河杭州段从中心城区流经城乡交错带, 景观类型复杂, 人类活动强烈, 对该区域河流水环境具有重要影响。

2 研究方法 2.1 主要数据采集及处理

(1) 遥感数据解译

以2017年高分二号(GF-2)影像为基础信息源, 依据京杭运河杭州段城市土地利用实地调查结果, 将景观类型划分为10个类型(表 1)；利用GF-2遥感影像, 通过eCognition面向对象的分类法对遥感影像进行解译；在此基础上, 为提高解译精度和分类的准确度, 消除时间差, 结合2019年谷歌实时地图和野外实地考察, 对分类结果进行修改和校正, 使得数据匹配。通过野外实地采集的GPS定点信息进行验证, 满足数据分析要求。

通过河网提取、河流分割等步骤, 并结合《浙江省河流手册》^[27]中的流域和河流水系分布图和杭州市行政边界, 确定研究区边界范围(图 1)。

(2) 河流水质监测

选择子流域界线相对清晰的支流作为监测河流, 在汇入上级河道前且不受上一级河道影响的位置进行水质采样, 共设置了72个具有代表性的水质监测点(图 1), 基本实现了监测点在研究区的均匀分布。分别在2019年3月、7月和11月进行了3次水样采集, 采样选择当月中旬连续晴天且采样前三天无降水的时段, 以消除降水对监测结果的影响；样品采集时保证采样位置处于河流中央且流速较大处, 以保证样品的代表性；水样储存于预先处理过的500 mL聚乙烯瓶中。

综合考虑京杭运河杭州段水环境主要污染因子, 本研究选取了氨氮(NH₄⁺-N)、硝氮(NO₃^--N)、总氮(TN)和总磷(TP)4个水质指标, 利用间断化学分析仪(BluVision^TM, Skalar)进行水质指标测定。

2.2 主要研究方法

(1) 缓冲区分析

由于研究区河渠交织密布, 水网发达, 河流人为调控程度高, 水文过程复杂, 河渠水流互通, 汇水单元界线不明确, 因此, 本研究以水质监测点为中心建立圆形缓冲区, 基于缓冲区分析景观格局对河流水质的影响。考虑到研究区位于杭州市的城乡过渡带, 景观异质性程度极高, 在缓冲半径划分时采用了200 m和500 m两种间隔距离, 共设置了8个不同半径的缓冲区, 缓冲半径依次为100、300、500、1000、1500、2000、2500、3000 m(图 2)。

(2) 景观格局指数计算

选择的景观格局指数包括斑块密度(PD)、最大斑块指数(LPI)、边缘密度(ED)、景观形状指数(LSI)、周长—面积分维数(PAFRAC)、景观蔓延度指数(CONTAG)、散布与并列指数(IJI)、景观分割度(DIVISION)、香农多样性指数(SHDI)和聚集度(AI), 以及各景观类型的面积百分比。通过Fragstats 4.2分别计算不同半径缓冲区内的格局指数。

(3) 统计分析

本文采用了Pearson′s相关分析法对景观指标和水质指标进行分析, 以探讨京杭运河杭州段不同尺度缓冲区上的景观类型及其空间格局与水质的相关性。一些研究已经将Pearson′s相关分析法作为衡量任何两个变量之间联系强度的标准^[33-34]。为了确定景观格局与水质指标之间的关系, 本文使用冗余分析对景观格局进行分析。为了定量研究京杭运河杭州段水质变化和景观格局之间的关系, 因此使用多元逐步回归分析法来探讨京杭运河杭州段对不同缓冲区内景观格局和水质的响应关系。

本研究采用SPSS 25.0实现Pearson′s相关分析和多元逐步回归分析, 冗余分析则通过CANOCO 5.0软件实现。

3 结果与分析 3.1 研究区景观格局分析

京杭运河杭州段景观格局如图 3所示, 研究区总面积1509 km², 其中农业用地和林地是2种主要的景观类型, 分别占研究区总面积的25.21%和24.87%, 是由于研究区北部区域属于城乡过渡带, 农业种植用地分布集中, 而西部属于著名的西湖、灵隐寺和西溪湿地风景区, 人类开发程度较低, 林地成大面积连片分布；居民用地、道路和商业与工业用地等建设用地, 分别占研究区总面积的18.74%、10.05%和9.24%, 主要是位于研究区中南部钱塘江北岸的城市建成区, 也是研究区中人类活动强度最高的区域；另外, 研究区河网密布、水系发达, 水体面积占到了7.74%。

人类活动不仅增加了景观斑块的种类, 而且使得斑块的形状更趋于多样化、复杂化, 从而加剧了景观的破碎化程度^[35]。在整个研究区尺度上, 边缘密度为129.69 m/hm², 周长—面积分维数为1.48, 景观形状指数为127.46, 充分体现了区域城乡交错带的典型特征, 斑块形状复杂多变, 斑块破碎化程度高；散布与并列指数为78.03, 聚集度为67.59, 表明研究区斑块分布较为分散, 聚集程度较低；香农多样性指数为1.83, 说明该区域景观类型复杂多样, 景观异质性高, 人类活动强度高低变化明显。

3.2 不同半径水质监测点缓冲区景观格局变化

图 4为不同半径缓冲区的景观类型组成比例, 可以看出不同缓冲区内景观组成差异明显。例如, 农业用地所占比例在100 m缓冲区内最低, 在300 m缓冲区内最高, 说明农村近距离沿河两岸的农田分布模式最为广泛；道路所占比例却在100 m缓冲区内最高, 充分表明了道路与河流堤岸并排相连分布的特征；居民用地比例在500 m缓冲区内最高, 表明居民用地在河流周边地区分布较广泛。

景观格局指数随着缓冲半径的变化发生明显变化(图 5)。斑块密度(PD)与边缘密度(ED)具有相同变化趋势, 在1000 m以内的缓冲区, 均随着缓冲半径的增大而降低, 且其中在100 m缓冲区中远高于其他缓冲区, 但当缓冲半径超过1000 m后仅有小幅波动, 基本维持稳定, 这表明城市河流离岸越近, 受人类活动的调控和干扰就越强烈, 景观破碎化程度越高, 特别是在近岸100 m的区域内尤为剧烈。周长—面积分维数(PAFRAC)和香农多样性指数(SHDI)与PD和ED的变化趋势正好相反, 当缓冲半径小于1000 m时呈增加趋势, 高于1000 m后基本稳定, 这也是由于越靠近河流两侧的近岸区域受人类调控越强烈, 使得斑块形状规则简单, 但景观类型更加丰富多变。散布与并列指数(IJI)在100 m缓冲区最高, 达到了76.4%, 表明景观类型的混合程度在邻接河岸的100 m区域内最高, 而在其他缓冲区则相对稳定(68%左右)。聚集度指数(AI)随着缓冲半径的增大而呈逐步降低趋势, 这主要是河流水体和沿河两侧道路增强了缓冲区内的景观聚集特征。

从景观格局指数在不同半径缓冲区内的变化趋势来看, 基于河流监测点的100 m和1000 m的缓冲半径是两个关键尺度, 在100 m缓冲区内人类活动的干扰最为强烈, 而这种干扰随着缓冲半径增大到1000 m后基本进入稳定状态, 缓冲距离对景观格局特征的影响不明显。

3.3 河流水质统计特征分析

京杭运河杭州段72个监测点的水质指标的描述性统计特征如表 2所示。氨氮含量在11月份的变幅最大(0.05—8.97mg/L), 变异系数达到了0.96, 从平均值来看, 11月最低(1.32 mg/L), 3月和7月较为接近, 分别为1.71 mg/L和1.77 mg/L；硝氮的平均含量在7月最低(1.14 mg/L), 在3月和11月的含量较为接近, 分别为1.76 mg/L和1.68 mg/L；总氮含量的平均值在3月(7.57 mg/L)远高于7月(3.58 mg/L)和11月(3.45 mg/L)；总磷含量呈逐步降低趋势。总体来看, 研究区河流的氨氮和总氮含量整体处于较高水平, 水体污染较为严重。因此, 京杭运河杭州段水质污染主要属于氮源污染。

3.4 景观格局对河流水质的影响分析 3.4.1 景观类型与河流水质的相关性分析

通过对72个监测点不同半径缓冲区内的景观类型比例与水质指标的相关性分析(表 3), 可以看出不同的景观类型对河流水质在不同的季节有着不同的影响。氨氮在汛期(7月)与居民用地和公共设施用地等建设用地的相关性最显著, 同时农业用地也对氨氮有着显著影响；而在非汛期(3月和11月), 氨氮与景观类型的相关性较低；氨氮与景观类型的关系呈现出了较强的季节性变化。硝氮与景观类型的关系最为密切, 其中居民用地在3个时段都与硝氮存在显著相关性；农业用地对硝氮的贡献主要表现在非汛期, 林地在11月份降雨量较低的情况下对硝氮有显著的消减作用。总氮含量受商业与工业用地的影响在3个时段都较为明显, 同时在7月和11月受农业用地的影响显著, 特别是7月份与公共设施用地的相关性极为显著；这也表明了工业和商业生产是总氮的主要来源, 同时农业生产也对总氮有重要影响。总磷与景观类型的相关性整体较弱, 仅3月份与草地和居民用地呈现相关性。京杭运河杭州段河流水质的污染源主要来自建设用地和农业用地。

3.4.2 景观格局与河流水质的冗余分析

景观格局指数随缓冲半径的变化趋势表明, 当缓冲半径大于1000 m时, 缓冲区的景观格局指数趋于稳定, 不再随缓冲半径的增大而发生较大变化, 因此, 本研究以1500 m缓冲区为重点, 基于冗余分析(RDA)量化景观格局指数对河流水质的解释能力排序(图 6), 图中各个景观格局指标箭头与水质指标箭头之间夹角小于90度, 两者关系为正相关, 大于90度为负相关, 等于90度则不相关；两者之间锐角越小或钝角越大, 则正向或负向相关性越强^[36]。结果表明, 散布与并列指数(IJI)与硝氮和总氮在3个月份都呈极强的正相关, 表明景观类型的混合程度越高, 河流硝氮和总氮含量也越高；周长—面积分维数(PAFRAC)、景观形状指数(LSI)和边缘密度(ED)与硝氮和总氮在非汛期(3月和11月)呈负相关, 在汛期(7月)却呈明显正相关, 说明相对复杂的斑块形状在非汛期有利于消减河流的硝氮和总氮污染, 但在汛期增大了污染物经降水冲刷进入河流的可能性；聚集度指数(AI)对水质的影响趋势与形状指数相似, 斑块的集中连片分布在非汛期增加了氮源污染, 但在汛期有助于氮源污染削弱；香农多样性指数(SHDI)在3个时段都与总氮呈明显正相关, 说明景观类型越复杂, 对河流总氮的影响越大；斑块密度(PD)与硝氮在非汛期呈显著负相关, 在汛期成极显著正相关, 城市景观的破碎化程度在非汛期加剧了河流硝氮污染, 但在汛期可以减少大面积连片斑块被降水冲刷从而有助于减缓硝氮污染。总磷含量在3个时段与景观格局指数的相关性均较低。

3.4.3 景观格局与河流水质的多元逐步回归分析

(1) 景观类型与水质之间的关系

表 4为水质指标和景观类型占比之间的多元逐步回归分析结果, 可以看出各个回归模型都达到显著水平(P < 0.05), 景观类型中商业与工业用地、居民用地、水体是主要影响水质指标的类型。在1500 m缓冲区, 不同月份的水质受到不同景观的影响。3月, 硝氮受到了商业与工业用地和农业用地的影响, 随着二者的增加硝氮浓度增加；总氮浓度是随着商业与工业用地的增加而增加, 随着水体增加而减少, 说明了城市的排水系统对水质指标有一定影响；总磷受到了居民用地和草地的影响。7月, 氨氮随着居民用地增加而增加, 居民生活对氨氮影响明显；硝氮随着公共设施用地、草地和湿地增加而增加, 居民用地的减少而减少；总氮随着商业与工业用地和公共设施用地增加而增加。11月, 氨氮随着水体的增加而减少；硝氮随着湿地和居民用地的增加而减少, 与商业与工业用地则相反；总氮随着农业用地和商业与工业用地的增加而增加。

(2) 景观格局指数与水质的关系

逐步多元回归分析能够将回归方程中那些不重要的自变量因子排除, 从而建立最佳模型。在1500 m的缓冲区尺度上, 景观格局指数与水质的逐步多元回归分析结果如表 5所示。3月, 硝态氮受到景观指数IJI的影响, 随着IJI的增加而增加；总氮随着PAFRAC的增加而减少。7月, 氨氮随着DIVISION的增加而增加, 随着IJI的增加而减少；总氮随着PAFRAC、IJI、PD、CONTAG和SHDI的增加而增加。11月, 氨氮随SHDI的增加而增加；硝氮随着CONTAG和IJI的增加而增加；总氮随着LPI的增加而减少, 随着IJI的增加而增加。

4 讨论

水质对于维持健康的河流生态系统的结构和功能至关重要, 自然过程和人为活动的影响使得河流水质和环境健康发生变化。京杭运河杭州段总面积1509 km², 城市景观和农村景观镶嵌分布, 北部以农村景观为主, 中南部为钱塘江北岸城市景观, 西部则为风景名胜区, 农业用地和林地分别占到了研究区总面积1/4, 居民用地、商业与工业用地等建设用地占40%, 水体面积占到了7.74%, 河网水系发达。本研究表明, 京杭运河杭州段水质污染类型主要是氮源污染, 污染源主要来自建设用地和农业用地。总氮与商业与工业用地之间存在着强烈的相关关系, 这表明建设用地是京杭运河杭州段中氮素的主要来源。有研究表明, 农业用地和城镇建设用地是河流水化学污染的重要来源^[37]。城市化进程突出表现在大面积的不透水硬表面取代了自然植被为主的透水软表面, 不透水地面增加、绿地和水体景观减少、景观破碎和离散化加剧, 由此引起了城市景观格局变化, 导致降雨径流过程发生了质的变化, 一定程度上造成了城市水体质量下降^[38-39], 使得城市面临着愈加交织复杂的水环境危机^[38]。同时, 城市排水系统干扰了自然排水过程, 比如污水管网渗漏、地表水及地下水渗入、污水收集处理后异地排放等^[40-42], 对河流水质影响严重。城市景观格局对河流氮磷污染物的发生、迁移和转化产生重要影响, 人类活动干扰加剧了河流氮、磷污染, 在很大程度上增加了河流污染的复杂性、随机性和空间分异性^[43]。本研究从逐步多元回归分析可以看出, IJI是影响水质最重要的景观格局指数, 同时, 京杭运河1000 m缓冲区尺度的景观格局对水质影响最为显著, 因此, 在加强京杭运河杭州段水资源的保护与管理时, 需要重点对此缓冲区尺度的景观格局进行优化与调整。

景观空间格局对进入河流污染物的吸收和降解过程的影响存在差异性^[26], 不同土地利用类型对于降雨径流过程的影响也有所不同^[39]。城市景观格局的水环境效应存在时间尺度性, 即在夏季丰水期、春秋季平水期和冬季枯水期对水体污染的影响程度存在差异, 例如在珠江、太湖和南京28处湿地等水体中, 丰水期总氮(TN)、总磷(TP)、化学需氧量(COD)、叶绿素浓度最高, 水环境污染程度最高^[14]。在本研究中京杭运河杭州段河流水质与PD、ED、LSI、PAFRAC和AI等景观格局指数的关系, 在汛期和非汛期差异较大甚至相反, 表明降水和径流等水文过程的关键性作用, 而景观格局对河流水质的影响在很大程度上是通过对径流的空间再分配等水文过程作用于河流水环境。因此, 针对京杭运河杭州段的河流水质污染防治, 需要结合运河两侧景观格局, 因地制宜通过格局优化和空间配置强化城市景观对污染物的消减和拦截等功能。

5 结论

本文以京杭运河起点杭州段河网水道为研究对象, 设置了72个河网水质监测点, 以水质监测点为中心划分圆形缓冲区, 在分别对该区域景观格局和河流水质监测分析的基础上, 研究景观类型和景观格局特征对河流水质在不同时期的影响特征, 结论如下：

(1) 在研究区整体尺度上, 景观破碎化程度高, 斑块形状复杂多变, 景观类型多样, 异质性高, 城乡交错带特征典型, 人类活动强度高低分化明显。

(2) 在缓冲区尺度上, 景观格局指数随缓冲半径内的变化特征明显, 100 m缓冲区内的格局指数远高于或低于其他缓冲区, 而当缓冲半径达到1000 m以后, 格局指数则相对稳定, 不再随缓冲半径的增加而发生明显变化, 体现了人类活动沿河流两侧向外扩展的特征。

(3) 京杭运河杭州段水系的氮源污染物含量整体处于较高水平, 氨氮、硝氮和总氮的年均含量分别为1.60、1.53、4.87 mg/L, 总磷含量处于低水平, 年均值为0.32 mg/L。

(4) 景观类型与河流水质的相关性在不同半径的缓冲区内表现较为一致, 硝氮受景观类型的影响最为显著, 氨氮仅在汛期与建设用地和农业用地的相关性较高, 总氮与景观类型相关性的季节性差异最明显, 而总磷与景观类型的相关性整体较弱。

(5) 城市景观格局指数与河流水质参数的相关性较为明显, 景观格局指数对不同的水质指标在不同的时期具有较好的解释能力, 且同一景观指数在汛期和非汛期对水环境的影响作用差异较大甚至相反, 景观格局在河网密集型城市区域对河流水环境影响的整体利弊还需进一步权衡量化。