水资源是人类生存和发展不可替代的重要资源，而水质对区域生态环境、农业发展等有重要影响[1],然而随着城市化进程加快、人口数量增加，水质恶化等问题突出[2]。景观格局受人类活动和自然环境的影响，景观格局也影响着水体内生物活动和水文过程，进而影响水质。景观格局对非点源污染物的产生、转移和变换过程起着重要作用[3-5]。因此，开展景观格局指数与水质相关性分析对于土地利用管理和水资源保护具有重要意义[6]。

景观格局指数与水质相关性研究是生态学研究的热点之一[7-9],研究主要分为两方面：一是探究水质与景观格局指数之间的相关性，二是利用景观格局指数分析不同土地利用类型对水质的影响[10-12]。目前关于不同空间尺度景观格局指数与水质相关性进行了大量研究，如：方娜等[13]研究鄱阳湖湿地水质发现，小流域尺度土地利用类型与水质的相关性低于缓冲区尺度的；曹灿等[14]研究发现艾比湖区300 m缓冲区内景观格局指数与水质的相关性最强。

密云水库作为北京市重要地表饮用水水源地，对保障北京市水资源调度配置和供水安全发挥着不可替代的作用。目前，不少学者开展了密云水库水质研究[15],但关于密云水库不同空间尺度下景观格局指数与水质相关性的研究较少。因此，本文分析2015年、2018年、2020年密云水库上游水质空间变化特征，并分析不同空间尺度下景观格局指数与水质的关系，精准识别水源保护的关键区域，以期为密云水库上游水源保护的精细化管控提供支撑。

1、材料与方法

1.1 研究区概况

密云水库上游流域面积1.58万km2,其中北京市所占面积为0.35万km2,涉及密云区、延庆区、怀柔区；河北省所占面积为1.23万km2,涉及张家口市沽源县、赤城县，承德市兴隆县、滦平县、丰宁满族自治县。流域地势西北高、东南低，以山地为主，其西部为白河流域，东部为潮河流域，主要有潮河、白河、天河、黑河、汤河5大支流。流域属温带季风性气候区，土壤类型主要为褐土和棕壤土，降水集中在6—8月，多年平均降水量为494 mm[16]。

1.2 数据来源及处理

研究所采用水质数据来源于实测数据，其中监测点主要位于河流与山区沟道交汇处，能够反映产汇流特点。共设17个监测点，其中白河6个、黑河2个、汤河2个、潮河7个。水质监测指标包括溶解氧(DO)、总氮(TN)、总磷(TP)、氨氮(NH

-N)浓度及高锰酸盐指数(CODMn)。每年8月进行水样采集，其中溶解氧采用《水质 溶解氧的测定 碘量法》(GB 7489—1987)中方法测定，总氮采用《水质 总氮的测定 碱性过硫酸钾消解紫外分光光度法》(HJ 636—2012)中方法测定，总磷采用《水质 总磷的测定 钼酸铵分光光度法》(GB/T 11893—1989)中方法测定，氨氮采用《水质 氨氮的测定纳氏试剂分光光度法》(HJ 535—2009)中方法测定，高锰酸盐指数采用《水质 高锰酸盐指数的测定》(GB/T 11892—1989)中方法测定。土地利用数据来源于中国30 m土地覆被数据集，土地利用类型分为耕地、林地、灌木、草地、水域、未利用地、不透水面7种[17]。

运用ArcGIS建立17个水质监测点的缓冲区，缓冲区半径选取100、200、500、1 000、1 500、2 000、2 500 m, 通过与土地利用数据叠加，获取不同尺度缓冲区的土地利用类型，并利用Fragstats4.2软件计算不同尺度缓冲区的景观格局指数，通过相关性分析和冗余分析探究密云水库上游不同空间尺度景观格局指数与水质的关系。

1.3 研究方法

1)景观格局指数。

根据研究流域实际情况，从景观丰富度、景观破碎度、景观聚集度及景观多样性等角度选取10个景观格局指数[6],见表1。利用ArcGIS和Fragstats4.2软件计算各景观格局指数。

表1 景观格局指数

2)统计分析方法。

利用SPSS软件对不同尺度缓冲区景观格局指数与水质指标进行相关性分析。为了更直观显示景观格局指数对水质的影响，利用Canoco5软件进行冗余分析(RDA),进一步揭示水质与景观格局指数的关系。

2、结果与分析

2.1 不同尺度缓冲区土地利用类型

密云水库上游不同尺度缓冲区土地利用类型面积计算结果表明，耕地、林地、草地面积占比较大；缓冲区半径为100、200 m时，主要土地利用类型为耕地和不透水面；缓冲区半径为500 m时，主要土地利用类型为耕地；缓冲区半径为在1 000、1 500、2 000、2 500 m时，林地、草地、耕地3种土地利用类型面积之和占缓冲区总面积的90%以上。随着缓冲区半径的增大，耕地、水域、不透水面面积占比呈减小趋势，林地和草地面积占比呈增大趋势，其中水域在所有缓冲区中面积占比都较小。2015年、2018年、2020年不同尺度缓冲区土地利用类型发生变化，但变化不大。除半径100 m缓冲区外，不同尺度缓冲区耕地面积呈先减小后增大变化趋势；林地、水域、不透水面面积有不同程度增大，在半径为2 500 m缓冲区分别增加了8.4%、13.7%、18.2%;草地面积呈减小趋势，在半径为2 000 m缓冲区减少了9.6%;未利用地、灌木面积变化相对较小，未利用地基本保持不变，灌木面积呈减小趋势。

2.2 不同尺度缓冲区景观格局指数

对不同尺度缓冲区景观格局指数进行分析处理，结果见图1。面积-密度-边缘类指数主要描述景观中斑块的数量和面积，以及斑块的边缘长度，LPI最大值出现在半径为200 m的缓冲区，最小值出现在半径为2 500 m的缓冲区，表明半径2 500 m缓冲区人类活动对景观的干扰较大。随着缓冲区半径逐渐增大，PD逐渐减小，表明随着缓冲区半径的增大，研究区景观破碎程度减小。NP随着缓冲区半径的增大而增大。

图1 密云水库上游不同尺度缓冲区景观格局指数

蔓延度与离散度指数反映不同类型斑块之间的空间分布情况以及混合状况。LSI、DIVISION和COHESION变化趋势大致相同，3种指数均随着缓冲区半径的增大而增大，说明研究区斑块分布复杂程度、景观类型连通性和不同类型斑块之间的分离程度随着缓冲区半径的增大而增大。随着缓冲区半径的增大，AI、CONTAG整体变化较小，表明不同半径缓冲区优势景观连通性差别不大。

多样性指数SHDI、均匀度指数SHEI描述景观组成的丰富度和均匀度，不受斑块空间格局的影响。不同半径缓冲区SHDI、SHEI相差不大，即不同半径缓冲区景观组成丰富度和均匀度接近。半径为2 000 m的缓冲区SHDI均值最大，说明半径为2 000 m的缓冲区景观类型更多样。半径为500 m的缓冲区SHEI均值最大，说明半径500 m缓冲区景观类型更均衡。

2.3 水质分析

1)密云水库上游水质分析。

根据密云水库上游17个水质监测点监测数据，统计得到2015年、2018年、2020年研究区水质特征，见表2。密云水库上游，NH

-N和TP平均值满足地表水Ⅱ类水质标准，且呈减小趋势；DO平均值满足地表水Ⅱ类水质标准，且呈增大趋势；CODMn平均值先减小后增大，CODMn平均值满足地表水Ⅱ类水质标准；TN污染较为严重，TN平均值逐渐增大。从变异程度看，密云水库上游TP、NH

-N、DO、TN、CODMn属于中等变异(10%≤变异系数≤100%),说明5种水质指标具有很强的空间分异特征。

表2 密云水库上游水质特征

2)密云水库上游水质时空变化。

根据《地表水环境质量标准》(GB 3838—2002),对2015年、2018年、2020年5种水质指标进行时空变化分析，结果见图2。除潮河外，密云水库上游DO波动较小，呈增大趋势，2020年DO满足地表水Ⅰ类水质标准的监测点占94%。TP呈减小趋势，2020年满足地表水Ⅰ类水质标准的监测点占53%,所有监测点均满足地表水Ⅱ类水质标准。NH

-N呈减小趋势，满足地表水Ⅰ类水质标准的监测点占比从2015年的53%增大到2020年的100%,2015年、2018年、2020年大部分监测点满足地表水Ⅰ类水质标准，所有监测点均满足地表水Ⅱ类水质标准，其中潮河上游和白河上游NH

-N污染较严重。CODMn呈先减小后增大趋势，2015年、2020年大部分监测点满足地表水Ⅱ类水质标准，其中潮河和白河相较于其他支流CODMn污染严重。对比其他水质指标，TN污染严重，82%以上监测点TN超过地表水Ⅴ类水质标准，且TN呈逐年增大趋势。

2.4 水质与景观格局指数关系分析

1)不同空间尺度景观格局指数与水质相关性分析。

采用Pearson相关性分析方法对密云水库上游不同空间尺度景观格局指数与水质的关系进行分析，结果见图3,其中不同半径缓冲区中长椭圆形从左至右分别代表NP、SHDI、LPI、COHESION、DIVISION、SHEI、CONTAG、AI、PD、LSI;红色代表正相关，蓝色代表负相关，颜色越深相关性越强；*表示在0.05水平(双侧)显著相关。

图2 密云水库上游水质空间分布

图3 不同尺度缓冲区景观格局指数与水质相关性

NP与PD分别反映景观中斑块数量和破碎度，二者与TN在半径为1 000、1 500、2 000、2 500 m的缓冲区显著正相关，说明景观组成的分裂程度和破碎度影响密云水库上游水体氮的浓度。LPI反映最大斑块面积占比，在不同半径缓冲区，LPI除与NH+4-N正相关外，与其他指标大多为负相关，这与密云水库上游的优势景观类型是林地和草地有关，当优势景观类型为“汇”(能够拦截或消纳外来污染物的景观类型)时，水质会得到较大改善。

LSI反映景观中斑块形状特征，AI反映景观中斑块聚集程度。在半径为2 000、2 500 m的缓冲区，LSI与TN显著正相关，AI与TN显著负相关，斑块边界越不规则、景观破碎化越严重，水体越容易被污染，这一结果与乔郭亮等[6]在苏南地区的研究结论一致。CONTAG代表不同类型斑块的蔓延程度，反映景观之间的离散程度，其与TN、CODMn负相关，主要原因是优势景观林地和草地可以减轻水体污染情况。COHESION代表斑块之间的连通性，在半径为1 500、2 000、2 500 m的缓冲区与TN显著负相关，表明随着COHESION增大，景观朝着高连通方向发展，进而加剧TN的污染。

SHDI和SHEI反映各斑块类型的异质性与景观类型的丰富度、均匀程度。景观类型越多，河流流经“源”景观(能够促进污染物释放或引入外来污染物的景观类型)的可能性越大，水质下降的可能性越大。SHDI越大意味着人类活动对当地生态环境干扰强度越大，而河流水质变差与人类活动干扰有很大关系。SHDI、SHEI与TN、正相关，这一结果与刘可暄等[4]认为景观类型的复杂程度对水质影响较显著的结论一致。

2)不同空间尺度景观格局指数与水质的冗余关系分析。

运用Canoco5软件分析不同尺度景观格局指数与水质指标的冗余关系，2020年水质除趋势对应分析(DCA)结果见表3。所有排序轴中，最长梯度值为0.8(小于3.0),因此采用冗余分析(RDA)研究水质与景观格局指数的关系是适宜的。半径为2 000 m的缓冲区景观格局指数对水质的解释度最高，为47.0%,该尺度缓冲区第一轴特征值、相关系数分别为0.30、0.924,其排序图能够很好地反映景观格局指数与水质的关系，半径为2 000 m的缓冲区是改善水质的关键区域，这一结果与陈优良等[18]对东江源的研究结论一致。

表3 2020年水质除趋势对应分析结果

密云水库上游半径为2 000 m的缓冲区水质与景观格局指数RDA结果见图4,图中两类箭头之间夹角代表变量之间的相关性，角度越大相关程度越高，夹角的余弦代表相关系数，夹角大于90°为负相关，夹角小于90°为正相关，横纵轴均为相关系数，其中值越大相关性越高。DO、TN与LSI、PD、NP、DIVISION、SHDI、SHEI正相关，与其他景观格局指数负相关；TP、

-N与LPI、COHESION、CONTAG、AI正相关，与其他景观格局指数负相关；CODMn与AI、DIVISION、SHDI、SHEI正相关，与其他景观格局指数负相关；景观格局指数对CODMn与TN的影响较大。

图4 水质与景观格局指数RDA结果

3、讨论

1)景观格局指数与水质的相关性。

最大斑块指数LPI反映最大斑块面积占整个景观面积的比例。康文华等[19]研究以林地为主要景观的地区发现，LPI与污染物含量负相关；郭玉静等[20]研究以耕地为主要优势景观的普者黑岩溶湖泊湿地发现，LPI与水体污染程度正相关。表明当优势景观为“源”景观时，LPI与水质恶化有显著关系；当优势景观为“汇”景观时，则有利于水质改善。对于CONTAG、SHDI,也不能简单地从数值大小判断其对水质的影响，需要结合研究区优势景观类型一同分析[21]。王志芳等[22]研究表明，景观多样性丰富，若“源”景观(建筑用地、耕地)为优势景观类型，则会导致水质恶化；若“汇”景观(林地、草地)为优势景观类型，则有利于水质改善。

2)景观格局指数与水质相关性研究的最佳缓冲区尺度。

不同尺度缓冲区在景观类型组成与结构上具有异质性，导致与水质的相关性存在差异。密云水库上游景观格局指数与水质的相关性随着缓冲区尺度的不同而不同，大尺度缓冲区景观格局指数与水质相关性较强，最佳尺度为2 000 m, 这与以往研究结果一致，比如东江源半径为2 000 m的缓冲区[23]景观格局指数能够解释与水质的相关关系。但部分学者利用河岸带缓冲区分析景观格局指数与水质相关性的最佳尺度与本研究有所不同，河岸带缓冲区最佳尺度集中在100、200 m等[4,24]。出现差异的原因可能是，某些土地利用数据分辨率过低，圆形缓冲区不能充分诠释采样点周围景观格局。小尺度河岸带缓冲区能够有效诠释整条河流沿线景观格局，但河流过长就无法聚焦采样点周围景观格局指数与水质之间的关系，在后续研究中可以将圆形缓冲区和河岸带缓冲区结合起来进行景观格局指数与水质相关性研究。

4、结论

通过相关性分析与冗余分析对密云水库上游不同尺度缓冲区景观格局指数与水质的相关性进行了研究，结论如下：密云水库上游水质采样点不同尺度缓冲区主要土地利用类型为耕地、林地、草地，随着缓冲区半径的增大，耕地面积占比逐渐减小，草地和林地面积占比逐渐增大；随着缓冲区半径增大，缓冲区内景观破碎程度降低、斑块数量逐渐增多、不同类型斑块之间的分离程度逐渐提高，不同尺度缓冲区优势景观连通性差别不大，景观组成丰富度和均匀度接近；密云水库上游水质整体呈好转趋势，但TN污染依旧严重，2015年、2018年和2020年大部分监测点TN浓度超过了Ⅴ类水质标准，其余水质指标基本达到Ⅱ类水质标准；通过冗余分析得出半径为2 000 m的缓冲区景观格局指数对水质的解释度最高，作为改善水质的关键区域应重点关注。

参考文献:

[1]朱爱萍,原升艳,温姗姗,等.不同空间尺度的景观格局对流溪河水质的影响[J].生态学报,2023,43(4):1485-1495.

[3]王飞,陶宇,欧维新.景观格局变化的水质净化服务响应关系研究进展[J].地球科学进展,2021,36(1):17-28.

[4]刘可暄,王冬梅,常国梁,等.多空间尺度景观格局与地表水质响应关系研究[J].环境科学学报,2022,42(2):23-31.

[5]范雅双,于婉晴,张婧,等.太湖上游水源区河流水质对景观格局变化的响应关系:以东苕溪上游为例[J].湖泊科学,2021,33(5):1478-1489.

[6]乔郭亮,周寅康,顾铮鸣,等.苏南地区景观格局特征与坑塘水质关联关系[J].农业工程学报,2021,37(10):224-234.

[7]黄金良,李青生,洪华生,等.九龙江流域土地利用/景观格局-水质的初步关联分析[J].环境科学,2011,32(1):64-72.

[8]张大伟,李杨帆,孙翔,等.入太湖河流武进港的区域景观格局与河流水质相关性分析[J].环境科学,2010,31(8):1775-1783.

[10]夏叡,李云梅,王桥,等.京杭大运河无锡段水质和土地利用的响应关系[J].自然资源学报,2011,26(3):364-372.

[12]杨洁,许有鹏,高斌,等.城镇化下河流水质变化及其与景观格局关系分析:以太湖流域苏州市为例[J].湖泊科学,2017,29(4):827-835.

[13]方娜,刘玲玲,游清徽,等.不同尺度土地利用方式对鄱阳湖湿地水质的影响[J].环境科学,2019,40(12):5348-5357.

[14]曹灿,张飞,阿依尼格尔·亚力坤,等.艾比湖区域景观格局与河流水质关系探讨[J].环境科学,2018,39(4):1568-1577.

[15]王丽娟,杨扬,郑娟娟,等.施肥对北京山区农田地表氮磷流失的影响:以密云水库流域为例[J].中国水土保持科学,2017,15(1):97-104.

[16]林青战,王汉男,韩娜娜,等.密云水库流域土壤动物群落组成及多样性[J].江苏农业科学,2016,44(12):445-448.

[18]陈优良,邹文敏,刘星根,等.东江源流域不同空间尺度景观格局对水质影响分析[J].环境科学,2022,43(11):5053-5063.

[19]康文华,蔡宏,林国敏,等.不同地貌条件下景观对河流水质的影响差异[J].生态学报,2020,40(3):1031-1043.

[20]郭玉静,王妍,刘云根,等.普者黑岩溶湖泊湿地湖滨带景观格局演变对水质的影响[J].生态学报,2018,38(5):1711-1721.

[21]杨强强,徐光来,杨先成,等.青弋江流域土地利用/景观格局对水质的影响[J].生态学报,2020,40(24):9048-9058.

[22]王志芳,程可欣.北运河流域雨洪“源-汇”景观时空演变[J].生态学报,2019,39(16):5922-5931.

[23]张敏,李令军,赵文慧,等.密云水库上游河流水质空间异质性及其成因分析[J].环境科学学报,2019,39(6):1852-1859.

[24]张微微,李晓娜,王超,等.密云水库上游白河地表水质对不同空间尺度景观格局特征的响应[J].环境科学,2020,41(11):4895-4904.

基金资助:北京市科技计划课题(Z2211000000112006,Z221100005222013);

文章来源:范霄寒,薛万来,王赛超,等.密云水库上游不同空间尺度景观格局指数与水质相关性研究[J].人民黄河,2024,46(11):114-120.