郴州市郴江河流域水质与土地利用关联的空间异质性分析
Spatial heterogeneity of the correlation between water quality and land use in the Chenjiang River basin, Chenzhou City
通讯作者: 汤玉奇(1986-),女,副教授,主要从事遥感影像处理的研究与教学工作。Email:yqtang@csu.edu.cn。
责任编辑: 李瑜
收稿日期: 2022-05-9 修回日期: 2021-10-7
基金资助: |
|
Received: 2022-05-9 Revised: 2021-10-7
作者简介 About authors
周石松(1997-),男,硕士研究生,主要从事自然资源遥感方面研究。Email:
为了研究郴州市郴江河流域水质与土地利用关联的空间异质性特征,采集了2021年9月流域水质与高分遥感数据; 通过不同类型、不同尺度的流域缓冲区划分,计算土地利用组分及各地类斑块形状、尺寸及分布等景观配置; 利用单因子、多因子相关性分析,进行区域土地利用对水质的空间异质性影响分析。结果表明: ①在水质与土地利用单因子空间异质性关联分析中,圆形缓冲区的普适性优于河岸缓冲区,而河岸缓冲区在分析建筑用地占比对各水质参数的影响方面起补充作用。在多因子空间异质性关联分析中,河岸缓冲区土地利用对水质参数的解释率普遍高于圆形缓冲区,圆形缓冲区土地利用对水质参数的总解释率最大值出现在半径700 m尺度。②不同土地利用组分对水质的影响不同,耕草地、林地、水域面积占比与水质呈正相关关系,建筑物、裸地呈负相关关系,其中耕草地与建筑用地对水质影响最大; 在进行土地利用规划时,要考虑到土地利用组分的影响,合理分配人工地物和自然地物面积占比。③不同土地利用景观配置对水质影响不同。其中反映斑块形状与尺寸的景观形状指数(landscape shape index,LSI)、最大斑块指数(largest patch index,LPI)对水质存在负面影响,反映斑块分布的散布与并列指数(interspersion and juxtaposition index,IJI)和水体氨氮(ammonia nitrogen,NH3-N)含量在小尺度内呈负相关。在进行土地利用规划时更要考虑到景观配置的合理性,例如要控制缓冲区内优势地类斑块大小。研究结论找出了最能影响水质的缓冲区大小,揭示了最能解释郴江河水质的土地利用组分及其景观配置,为郴江河水质恶化的科学治理措施提供选择依据,对水环境保护有一定意义。
关键词:
This study investigated the spatial heterogeneity of the correlation between water quality and land use in the Chenjiang River basin of Chenzhou City using the water quality and high-resolution remote sensing data of the river basin in September 2021. Based on the division of river basin buffer zones by different types and scales, this study calculated land use components and their landscape configurations, such as the shape, size, and distribution of patches. Then, this study analyzed the spatial heterogeneity effect of regional land use on water quality through single-factor and multi-factor correlation analyses. The results are as follows: ① In the single-factor correlation analysis of spatial heterogeneity between water quality and land use, circular buffer zones showed higher universality than riparian buffer zones, which played a complementary role in analyzing the effect of built-up land occupation on water quality parameters. In the multi-factor correlation analysis of spatial heterogeneity, the land use presented generally higher interpretation rates for water quality parameters in riparian buffer zones than in circular buffer zones, with the maximum total interpretation rate of land use for water quality parameters in circular buffer zones occurring on a scale of a radius of 700 m. ② Different land use components had different effects on water quality. The proportions of cultivated land, grassland, forest land, and water areas correlated positively with water quality, while those of building land and bare land correlated negatively with water quality. Among them, cultivated land, grassland, and building land had the most significant effects on water quality. Land use planning should consider the influence of land use components and reasonably allocate the proportions of artificial and natural features. ③ Different landscape configurations had different effects on water quality. The landscape shape index (LSI) and the largest patch index (LPI), which reflect the shape and size of patches, had negative effects on water quality. Furthermore, the interspersion and juxtaposition index (IJI), which reflects the distribution of patches, correlated negatively with the ammonia nitrogen (NH3-N) content in water bodies on a small scale. Land use planning should consider the rationality of landscape configurations. For example, it is necessary to control the sizes of dominant land patches in the buffer zones. This study identified the buffer zone size that influences water quality the most and the land use components and landscape configurations that can explain the water quality of the Chenjiang River the best. Therefore, this study provides a basis for selecting scientific control measures for the deterioration of water quality in the Chenjiang River, showing certain significance for water environment protection.
Keywords:
本文引用格式
周石松, 汤玉奇, 程宇翔, 邹滨, 冯徽徽.
ZHOU Shisong, TANG Yuqi, CHENG Yuxiang, ZOU Bin, FENG Huihui.
0 引言
水资源是人类社会生存发展的重要资源。2020年《全球风险报告》指出,水资源相关环境问题是当前世界面临最严重的危机之一[1]。随着中国工业化、城镇化的快速发展,水质恶化和水资源短缺的现象日益严重[2]。为了实现可持续发展目标,我国有关部门做出巨大努力,对重点河流断面进行定期水质监测并发布水质数据,点源污染基本得到有效治理,非点源污染逐渐成为影响水质的主要因素[3]。流经人类活动区域的河流水质,受到与人类活动密切相关的土地利用模式的影响更为严重[4]。土地利用对水质指标的影响机制一直是国内外研究的热点问题[5]。已有研究表明,土地利用组成及其景观配置在污染物的形成、释放、吸收、拦截中起着巨大作用,进而影响水质[6]。在研究方法上,根据不同尺度的具体情况,多使用典型关联分析(canonical correlation analysis,CCA)、冗余分析(redundancy analysis,RDA)和多元逐步回归模型等方法进行研究。张以晖[7]在黄浦江河岸带利用回归分析法、RDA对河岸带景观格局与河流水质进行关联分析得出,景观配置相较于土地利用组成更能解释黄浦江水质的变化; 赵鹏等[8]以广东省淡水河流域为研究对象,利用多元线性逐步回归模型、Spearman秩相关分析、CCA对水质指标与景观格局指数之间的响应关系开展研究,得出流域缓冲区内的景观格局指数对水质的解释率较高; 张志敏等[9]则在典型网状河网区域利用偏最小二乘回归得出雨季大多数水质指标受到园地、林草地和香农均匀度指数(shannon’s evenness index,SHEI)的重要影响。而现有的土地利用对水质的影响机制分析方法,大多基于较大缓冲区尺度内某些特定水质参数与土地利用之间的关联分析,其有效性已受到广泛认可; 但不同区域之间由于地形地质、数据精度、气候条件差异较大,因此有必要针对水体保护典型区域展开进一步分析研究。
不同类型、不同尺度缓冲区内的土地利用是解释水质变化的主要因素[10]。李艳利等[11]设置以采样点为圆心的缓冲区,研究水质与土地利用之间的关系; 陆君等[12]、焦胜等[13]和Li等[14]将研究区流域划分为子流域来研究土地利用对水质的影响; 汪昱昆等[15]研究了上海市市内河岸带土地利用对水质的影响,发现土地利用在半径为500 m和1 000 m河岸缓冲区内能较好地解释水质; 蔡宏等[16]利用地理信息系统技术,得出在子流域缓冲区和河岸缓冲区内,水质指标对各土地利用类型响应规律一致的结论; 黄金良等[17]的研究表明,九江流域大部分水质参数与土地利用在河岸缓冲区尺度上相关性更强。本文研究区为复杂小流域,子流域划分不明显,采用圆形缓冲区与河岸缓冲区更能反映采样点周边土地利用状况。由此看出,从缓冲区类型、尺度上综合分析研究区水质与土地利用的空间异质性关联,对当地水资源保护和国土资源优化管理意义重大; 而现有研究的缓冲区尺度间隔较大,导致对复杂小流域内土地利用对水质的影响分析不足。
综上所述,为了量化评估郴州市郴江河流域土地利用组分及其景观配置对水质的影响,开展顾及空间异质性特征的流域水质与土地利用相关性分析,本文分别采用多尺度圆形(半径100~1 500 m)、河岸(离河岸的垂直距离100 m和200 m)缓冲区进行流域水质与土地利用空间异质性关联分析; 并综合单因子相关分析、多因子相关分析等方法,对该流域水质与土地利用进行关联分析,实现对复杂小流域范围内水质与土地利用关联的深入分析,为国家可持续发展议程创新示范区实现其生态文明建设等可持续发展战略目标提供技术支持。
1 研究区概况与数据源
1.1 研究区概况
1.2 数据获取
1.2.1 土地利用分类
为了定量评估郴州市郴江河复杂小流域内的土地利用组分及其景观配置对水质的影响,并考虑到现有土地利用产品不足以反映复杂小流域范围内的地物细节,本文在参考30 m空间分辨率土地利用产品的基础上,利用支持向量机(support vector machines, SVM)分类器对2021年9月获取的研究区哨兵2号卫星遥感数据(10 m分辨率)进行土地利用分类。参考《土地利用现状分类》标准[23]并结合郴江河流域低矮作物与草地混合分布等实际情况,本文将研究区土地利用类型分为建筑用地、裸地、林地、水域和耕草地5类,总体分类精度约95 %,Kappa系数为0.94,满足该区域分析需求。同时,利用30 m分辨率数字高程模型(digital elevation model,DEM)数据划分研究区流域范围,并辅助水系提取。
1.2.2 水样采集
考虑到郴江河流域降水集中于每年6—9月,水体非点源污染也主要发生在此时段,本文研究团队于2021年9月15日赴郴江河流域研究区内进行了水样采集。综合分层抽样、空间分布和交通可达性等因素,共选取14个点位进行采样,采样点分布如图1所示。
图1
图1
郴州市郴江河流域水体采样点分布
Fig.1
Distribution of water sampling points in Chenjiang River basin, Chenzhou City
水体采样点的选取充分考虑了周围土地利用类型的多样性,并避开了生活污水及农业废水排放口。每个样点采集3份500 ml水样,为避免水样微生物繁殖及化学反应,水样采集后用冰袋冷藏保存并于当日送检。由于郴江河流域水体污染以非点源污染为主,生活污水和农业废水是主要污染源,本文选取总磷(total phosphorous,TP)、总氮(total nitrogen,TN)、氨氮(ammonia nitrogen,NH3-N)和化学需氧量(chemical oxygen demand,COD)等4种水质指标进行水样污染评估。其中: 水体TP大多来源于农业废弃物、化肥和生活含磷洗涤剂; TN和NH3-N是反映水体富营养化的主要指标; COD是反映水体有机污染的重要指标之一。当上述水质污染指标较高时,会对流域水生生物造成影响,对流域内居民健康造成极大危害。
2 数据处理分析流程
为了量化评估郴州市郴江河流域土地利用组分及其景观配置对水质的影响,开展顾及空间异质性特征的流域水质与土地利用相关性分析,本文在分别进行流域水质指标空间分布分析和土地利用空间异质性分析的基础上,最终实现对复杂小流域范围内的水质与土地利用空间异质性关联分析(图2)。
图2
图2
流域水质与土地利用关联的空间异质性分析流程
Fig.2
Analysis process of spatial heterogeneity of water quality and land use in watershed
3 结果与分析
3.1 流域水质空间分布分析
表1 水质项目含量分类标准
Tab.1
项目/ 分类 | I类 (源头水) | II类(生 活饮用水) | III类(渔 业用水) | Ⅳ类(工 业用水) | Ⅴ类(农 业用水) |
---|---|---|---|---|---|
TN ≤ | 0.20 | 0.500 | 1.00 | 1.5 | 2.0 |
TP ≤ | 0.01 | 0.025 | 0.05 | 0.1 | 0.2 |
NH3-N ≤ | 0.15 | 0.500 | 1.00 | 1.5 | 2.0 |
COD ≤ | 15.00 | 15.000 | 20.00 | 30.0 | 40.0 |
图3为研究区各水质参数浓度空间分布情况。从中可以看出,TN,TP和COD等水质指标均表现出人工地表区聚集性, 且TN,TP和NH3-N等水质参数出现部分超标的情况,在郴江河流域各水质采样点中,人工地表区域的采样点水质多为Ⅳ和Ⅴ类水,而自然地表采样点水质多为Ⅲ类或更优水质。其中,7—9号采样点的TP,TN和NH3-N指标相比其他人工地表区域采样点稍低,可能由于附近有郴州市最大的公园——苏仙湖公园,污染排放较少且公园绿地具有一定水质净化作用。郴江河流域水体主要污染包括TP和TN,呈现多区域超标现象,浓度已超过可能发生水华的标准(TN>0.20 mg/L和TP>0.02 mg/L)[25],建议关注并预防水华等水体污染事件的发生,加强水体保护及治理措施,以保证流域内居民正常生产生活及生态资源保护利用。
图3
图3
流域采样点水质参数空间分布
Fig.3
Spatial distribution of water quality parameters at sampling points in the watershed
3.2 土地利用空间异质性分析
如图1所示,本文选取的水质采样点遍及郴江河主干流域,包括郴州市及其周边的自然地表及人工地表区域,河流宽度较窄,干流宽度多为30~50 m之间,周边土地利用组分及景观配置复杂多变。因此,在流域水质与土地利用关联分析中,合理的缓冲区类型及其尺度设置是分析的关键。
同时,为了进一步分析郴江河流域水质与土地利用景观配置的空间异质性关联,本文在计算各种类型、尺度缓冲区内土地利用组分的基础上,通过筛选并定义景观格局特征指数集,对流域土地利用景观配置进行多类型、多尺度区域分析。考虑到郴江河流域土地利用类型复杂多样,各地类占比变化及其破碎化程度较大,因此本文结合现有研究中对水质影响较大的景观格局指数[28],定义了郴江河流域土地利用景观格局特征指数集: 景观形状指数(landscape shape index,LSI)、散布与并列指数(interspersion and juxtaposition index,IJI)和最大斑块指数(largest patch index,LPI)。其中,LSI反映土地利用斑块形状复杂程度及分布离散程度; IJI测算不同土地利用斑块类型之间的总体离散情况; LPI为最大土地利用斑块类型在整体景观中的面积占比,一定程度上反映了缓冲区主体土地利用类型对水质的影响程度[29]。
各采样点在不同类型、不同尺度缓冲区内的土地利用类型占比如图4所示。从中可以看出: ①在圆形缓冲区内,建筑用地、耕草地和林地是优势地类。在较小尺度圆形缓冲区内水域占比仅次于优势地类; 在较大尺度圆形缓冲区内,水域占比最低。②在河岸缓冲区内,耕草地和建筑用地是主要地类,各地类面积占比表现为耕草地>建筑用地>林地>水域>裸地。综上可见,郴江河小流域地类较为复杂,以及本文选取较为密集的缓冲区的必要性。
图4
图4
郴江河流域各空间尺度土地利用类型占比
Fig.4
Proportion of land use types at different spatial scales in the Chenjiang River basin
在多数采样点中,①LSI表现为圆形缓冲区>河岸缓冲区,人工地表区>自然地表区; ②IJI表现为河岸缓冲区>圆形缓冲区,人工地表区>自然地表区; ③LPI表现为圆形缓冲区>河岸缓冲区,人工地表区>自然地表区。综上可见,3种景观格局指数在不同缓冲区、不同地表区表现不尽相同,在对研究区土地利用景观配置对流域水质污染的影响分析中,需结合多类型、多尺度缓冲区进行空间异质性关联分析。
3.3 流域水质与土地利用的空间异质性关联分析
3.3.1 单因子关联分析
表2 水质参数与土地利用单因子相关性分析结果①
Tab.2
水质参数 | 缓冲区 类型 | 缓冲区 尺度/m | 土地利用组分 | 土地利用景观配置 | ||||||
---|---|---|---|---|---|---|---|---|---|---|
林地 | 建筑用地 | 裸地 | 水域 | 耕草地 | LSI | IJI | LPI | |||
TN | 圆形缓冲区 | 100 | 0.64* | |||||||
200 | 0.56* | |||||||||
300 | ||||||||||
500 | ||||||||||
700 | 0.56* | -0.55* | ||||||||
1 000 | -0.52* | 0.61* | -0.60* | 0.55* | 0.58* | |||||
1 500 | -0.51* | 0.60* | -0.64* | 0.66* | 0.59* | |||||
河岸缓冲区 | 100 | 0.66* | ||||||||
200 | 0.65* | -0.55* | 0.64* | |||||||
TP | 圆形缓冲区 | 100 | -0.51* | 0.65* | ||||||
200 | -0.50* | 0.70** | ||||||||
300 | ||||||||||
500 | ||||||||||
700 | ||||||||||
1 000 | 0.58* | -0.55* | 0.53* | 0.51* | ||||||
1 500 | -0.56* | 0.69** | 0.55* | |||||||
河岸缓冲区 | 100 | 0.68** | ||||||||
200 | 0.61* | 0.61* | ||||||||
NH3-N | 圆形缓冲区 | 100 | -0.53* | -0.53* | ||||||
200 | -0.63* | |||||||||
300 | 0.61* | -0.54* | ||||||||
500 | ||||||||||
700 | ||||||||||
1 000 | 0.51* | |||||||||
1 500 | ||||||||||
河岸缓冲区 | 100 | -0.67** | ||||||||
200 | 0.63* | -0.56* | ||||||||
COD | 圆形缓冲区 | 100 | 0.82** | -0.73** | 0.71** | |||||
200 | 0.79** | -0.86** | 0.52* | |||||||
300 | 0.73** | -0.84** | 0.59* | |||||||
500 | 0.67** | -0.87** | ||||||||
700 | 0.67** | -0.87** | 0.58* | |||||||
1 000 | -0.53* | 0.72** | -0.84** | 0.67** | 0.71** | |||||
1 500 | -0.57* | 0.74** | -0.83** | 0.69** | 0.56* | |||||
河岸缓冲区 | 100 | 0.84** | -0.85** | 0.67** | ||||||
200 | 0.79** | -0.88** | 0.63* |
①表中只列出了相关性显著的解释变量; * 表示 P < 0. 05; **表示 P < 0. 01。
1)流域水质与土地利用组分的空间异质性关联分析。表2左侧所示为各水质指标在不同类型、不同尺度缓冲区内与土地利用组分的相关性(仅列出显著相关部分,即显著性水平P值小于0.01或0.05)。可以看出: ①采样点水体TN含量与耕草地、林地在较大圆形缓冲区尺度(大于700 m)呈显著负相关,而与建筑用地在较大圆形缓冲区尺度和河岸缓冲区尺度内呈显著正相关,其中与半径1 000 m圆形缓冲区内的建筑用地占比相关度最高(r=0.61,P<0.05); ②采样点水体TP含量与较大圆形缓冲区(半径大于1 000 m)内的耕草地、较小圆形缓冲区(半径小于200 m)内的水域呈显著负相关,与半径1 000 m圆形缓冲区及河岸缓冲区内建筑用地呈显著正相关,其中与100 m河岸缓冲区内建筑用地占比相关度最高(r=0 68,P<0.01); ③采样点水体NH3-N含量与小尺度圆形缓冲区(半径100 m)内的水域呈显著负相关,与较大尺度圆形缓冲区(半径1 000 m)内的裸地呈显著正相关。总体来看,合理规划郴江河流域内,特别是大范围流域内的植被及人工用地占比,将有利于对水体TN,TP和COD含量的有效抑制。同时,对未利用地的合理规划也是降低水体NH3-N污染的关键。
2)流域水质与土地利用景观配置的空间异质性关联分析。表2右侧所示为各水质指标在不同类型、不同尺度缓冲区内与土地利用景观配置的相关性(同样仅列出显著性水平P值小于0.01或0.05的显著相关部分)。可以看出: ①TN,TP与较小圆形缓冲区(半径100~200 m)和200 m河岸缓冲区内的LSI、与较大圆形缓冲区(半径1 000~1 500 m)内的IJI和LPI呈显著正相关; ②NH3-N与半径300 m圆形缓冲区和200 m河岸缓冲区内的LSI呈显著正相关,与较小圆形缓冲区(半径100 m)和300 m河岸缓冲区内的IJI呈显著负相关,地类分布越分散,越能有效抑制水体NH3-N的污染程度; ③相较于其他水质参数,COD与景观格局指数的显著相关性更强。总体来看,在主体地类显著相关的情况下,流域内各地类斑块越散布,越易加剧水体TN,TP和COD污染[32]; 各地类分布越分散,越能有效抑制水体NH3-N的污染程度。因此,各土地利用类型空间分布的合理配置,对抑制郴江河流域水体污染至关重要。
3.3.2 多因子冗余关联分析
由于单因子分析无法反映多变量间的相互关系,为综合分析土地利用组分及其景观配置对郴江河流域水体水质的空间异质性影响,定量评估土地利用组分和景观配置总体及个体对水质污染的解释率,本文进一步进行了流域水质与土地利用多元线性回归分析。RDA作为经典的多因子相关性分析方法,既能得到单个环境因子对流域水质的方差贡献率,也能反映多个环境变量对所有水质指标的影响。因此,本文首先对流域水质指标进行趋势对应分析(detrended correspondence analysis, DCA),其梯度长度(gradients length)第一轴较小,故选择RDA进行多元线性回归关联分析[33]。
水质总体污染与各类型、各尺度缓冲区内的土地利用组分及其景观配置的冗余关联分析结果如表3所示。
表3 水质与土地利用多因子相关性分析结果①
Tab.3
缓冲区 类型 | 缓冲区 尺度/m | 土地利用组分 | 土地利用景观配置 | 总解释率 | ||||||
---|---|---|---|---|---|---|---|---|---|---|
林地 | 建筑用地 | 裸地 | 水域 | 耕草地 | LSI | IJI | LPI | |||
圆形缓冲区 | 100 | 0.24* | 0.39** | 0.62* | ||||||
200 | 0.25* | 0.21* | 0.27* | 0.71** | ||||||
300 | 0.23* | 0.22* | 0.24* | 0.63* | ||||||
500 | 0.24* | 0.28* | 0.69* | |||||||
700 | 0.31* | 0.36* | 0.75** | |||||||
1 000 | 0.37** | 0.39** | 0.31* | 0.63* | ||||||
1 500 | 0.36** | 0.43** | 0.25* | 0.35* | 0.28* | 0.68* | ||||
河岸缓冲区 | 100 | 0.46** | 0.31* | 0.27* | 0.76** | |||||
200 | 0.42** | 0.37* | 0.39** | 0.75** |
①表中只列出了相关性显著的解释变量; * 表示 P < 0. 05; **表示 P < 0. 01。
从表3可以看出,建筑用地、耕草地占比及斑块散布程度在多数尺度的圆形、河岸缓冲区内对水质总体污染状况均有影响,IJI及LPI在较大缓冲区内对水体污染造成影响。因此,合理调控主体地类占比及各地类、斑块分布将有利于抑制郴江河流域水体污染。
在空间异质性多因子冗余关联分析中,随着圆形缓冲区尺度的扩大,各类土地利用组分及其景观配置对水体污染的总解释率呈现出先增大后减小的趋势,并在半径700 m尺度上达到最大解释率75%。其中,建筑用地、耕草地占比对水质总体污染状况的影响随着圆形缓冲区尺度的扩大呈现增大趋势。在河岸缓冲区与圆形缓冲区内,其土地利用组分与其景观配置对水质总体污染的最优解释率相当,在多因子冗余关联分析中稳定性更高。
为了进一步可视化展示并分析研究区土地利用组分及其景观配置对水体污染的空间异质性影响,图5示出各类型、各尺度缓冲区内土地利用组分及其景观配置对水体总体污染的贡献程度及其差异性。其中,箭头长度反映贡献程度大小,箭头越长,说明该土地利用组分或景观格局特征对水体污染的贡献程度越强; 箭头间的夹角表示变量间的相关性,夹角余弦值越大,则两变量对水体污染的贡献越相关。从图5可以看出,在不同缓冲区类型及其空间尺度上,各土地利用组分及其景观格局特征与水体污染程度的相关性存在差异。在多数缓冲区类型及其空间尺度内,耕草地、建筑用地占比与水体污染相关性最强。其中,人工地物占比、主体地类占比和斑块散布程度对水质污染有加剧作用,而植被占比能在一定程度上抑制水体污染程度。进一步证明了人工地物与植被的合理配置对于郴州市郴江河流域水体污染治理的重要性。
图5
4 问题与讨论
4.1 郴江流域水质与土地利用空间异质性关联特征
在水质与土地利用组分的单因子空间异质性关联分析中,较大圆形缓冲区内各水质参数与土地利用组分显著相关项更多,但河岸缓冲区内水质与建筑用地占比相关性强于圆形缓冲区; 圆形缓冲区的普适性优于河岸缓冲区,河岸缓冲区在对建筑用地占比、斑块及地类散布程度对下游水体污染的影响分析方面具有补充作用。
在水质与土地利用多因子空间异质性关联分析中,河岸缓冲区和较小的圆形缓冲区内建筑用地、耕草地和LSI是主要的水质影响因子,较大的圆形缓冲区内建筑用地、耕草地和各景观格局指数都是主要的水质影响因子。河岸缓冲区土地利用对水质参数的解释率普遍高于圆形缓冲区,圆形缓冲区土地利用对水质参数的总解释率最大值出现在半径700 m尺度。土地利用组分与景观配置对水质的总解释率最大值出现在100 m河岸缓冲区内(解释率达到了76%),更能体现水质参数对土地利用的响应。该结论与周俊菊等[33]关于近距离河岸缓冲宽度内草地、建成区和耕地具有较强水质解释能力的结论相符。
4.2 土地利用组分对水质的影响
水质与土地利用组分空间异质性关联分析显示: ①在多数空间尺度上,建筑用地面积占比与COD,TP和TN浓度正相关性较大(图5),体现了人为作用对水质的负面影响[34]。人类生产生活导致林地、耕草地和水域等地类破碎化程度增大,生活垃圾堆积存放,地表污染物浓度增大,进而通过水土径流渗入河流中,使得各水质指标恶化。相较于圆形缓冲区,在河岸缓冲区尺度下水质与建筑用地的相关性更强,这是因为河岸两侧的生活污染更容易迁移进入河流,甚至直接排放进入河流[2]。②在多数空间尺度上,耕草地与各水质参数呈负相关关系,这一结果与耕草地被认为是非点源污染源释放化肥中的氮磷元素,从而恶化水质的结论不同[35]。其原因一方面可能与土地利用分类有关[36],在耕草地中,可能更多的地类是不易与耕地区分但施肥量更少的草地; 另一方面可能由于郴江河流域贯穿郴州市区,耕草地区域人类活动相对较弱,由此对水质带来的正面效应相比于耕草地所带来的污染更明显。③林地具有截留、吸收污染物从而改善水质作用[25],本文在除半径100 m圆形和河岸缓冲区可能因区域内林地占比过小、影响不明显之外,林地都显示出与水质参数的负相关关系; 水域面积在地表径流、降雨的作用下由于其对水质污染具有稀释作用[37],与水质参数呈现负相关; 而裸地由于占比较小,且多为废弃荒地,水土流失较严重,因此在多数空间类型、尺度内都与水质呈微弱负相关关系。
4.3 土地利用景观配置对流域水质的影响
5 结论
本文为了量化评估郴江河流域土地利用对水质的影响,划分2种缓冲区类型和多种缓冲区尺度,并在不同缓冲区内分别利用单因子相关分析、多因子相关分析等方法进行二者的空间异质性关联分析。实现对郴州市郴江河这一复杂小流域范围内水质与土地利用关联的深入分析,为郴州市建设国家可持续发展议程创新示范区、实现其生态文明建设等可持续发展战略目标提供技术支持,并得出如下结论:
1)在单因子空间异质性关联分析中,圆形缓冲区在水体污染含量与土地利用组分及其景观配置的单因子关联分析中的普适性优于河岸缓冲区,较大的圆形缓冲区水质各参数与土地利用更显著相关,河岸缓冲区起补充作用。
2)在多因子空间异质性关联分析中,不同缓冲区类型及其空间尺度上,各土地利用组分及其景观格局特征与水体污染程度的相关性、对水体污染的解释率存在较大差异,最能影响水质参数的土地利用是建筑用地和耕草地,林地、水域和耕草地对水质有正面作用,而建筑用地、裸地、LSI和LPI对水质有负面作用。
3)建议发挥林地在较大圆形缓冲区范围内对水质的积极作用,通过植树造林、退耕还林还草等措施来保护水质。在进行土地利用规划时也要避免河流附近土地利用类型过于单一,尤其是建筑用地。建议严格管控近河岸区的生活污水排行,并合理规划对水质参数解释率较大的河岸缓冲区内地块形状与尺寸。
综上所述,通过对郴江河流域水质与土地利用关联的空间异质性分析,理清了流域水质与土地利用关联分析的缓冲区类型及尺度选择依据,并研究了不同土地利用指数对水质的影响。在未来的研究中,将进一步使用深度学习方法进行精确土地利用分类,并细化区分耕地与草地对流域水质的影响。
参考文献
The global risks report 2020
[J/OL].
Landscape ecological approach to the relationships of land use patterns in watersheds to water quality characteristics
[J].
山地城市河流土地利用结构对水质的影响——以重庆市为例
[J].
Effects of land use pattern on water quality in mountainous city:A case study of Chongqing City
[J].DOI:10.18307/2016.0211 URL [本文引用: 1]
流域尺度上的景观格局与河流水质关系研究进展
[J].
Advances in the studying of the relationship between landscape pattern and river water quality at the watershed scale
[J].
丹江上游土地利用/景观指数与水质关系初探
[J].
A preliminary study on the relationship between land use/landscape index and water quality in the upper reaches of Danjiang River
[J].
Linking landscape connectivity and ecosystem service provision:Current knowledge and research gaps
[J].DOI:10.1007/s10021-013-9647-2 URL [本文引用: 1]
流域景观格局与河流水质的多变量相关分析
[J].
Multivariate correlation analysis between watershed landscape pattern and river water quality
[J].DOI:10.5846/stxb URL [本文引用: 1]
典型网状河网区域土地利用和景观格局特征对地表季节水质的影响——以江苏省溧阳市为例
[J].
Effects of land use and landscape pattern characteristics on surface seasonal water quality in typical networked river networks:A case study of Liyang City,Jiangsu Province
[J].
Spatially varying relationships between land use and water quality across an urbanization gradient explored by geographically weighted regression
[J].DOI:10.1016/j.apgeog.2010.08.001 URL [本文引用: 1]
浑太河流域多尺度土地利用/景观格局与水质响应关系初步分析
[J].
A preliminary analysis of the relationship between multi-scale land use/landscape pattern and water quality response in Huntai River basin
[J].
黄山市太平湖流域土地利用结构与河流水质相关性分析
[J].
Correlation analysis between land use structure and river water quality in Taiping Lake basin,Huangshan City
[J].
沩水流域土地景观格局对河流水质的影响
[J].
DOI:10.11821/dlyj201412005
[本文引用: 1]
以红壤丘陵区沩水流域为研究区域,运用GIS空间分析及统计分析的方法,探讨了流域土地景观格局对河流水质的影响。结果表明:① 林地面积百分比与电导率(COND)、综合污染指数(Pr)相关性指数分别达-1.000、-0.997,呈显著负相关;耕地面积百分比与COND、总磷(TP)及土地利用程度综合指数(La)与COND、Pr呈显著正相关。② 斑块数量(NP)、香农多样性指数(SHDI)与COND、Pr相关性指数均在0.997以上,呈显著正相关;平均斑块面积(AREA_MN)、蔓延度指数(CONTAG)与COND、Pr及最大斑块指数(LPI)与TP呈显著负相关。③ 从时间上看,所选用的大部分土地景观格局指标与pH值、溶解氧(DO)在丰水期与枯水期的相关性状态相反,且枯水期土地景观格局对河流水质的影响较丰水期更为显著。
Effects of land landscape pattern on river water quality in Shuishui basin
[J].
Effects of catchment and riparian landscape setting on water chemistry and seasonal evolution of water quality in the upper Han River basin,China
[J].DOI:10.1371/journal.pone.0053163 URL [本文引用: 1]
上海地区河网水质空间分异及对河岸带土地利用的响应
[J].
Spatial differentiation of water quality in river network in Shanghai and its response to land use in the riparian zone
[J].
基于遥感和GIS的赤水河水质对流域土地利用的响应研究
[J].
Response of Chishui River water quality to land use in the basin based on remote sensing and GIS
[J].
九龙江流域土地利用/景观格局-水质的初步关联分析
[J].
Preliminary correlation analysis of land use/landscape pattern-water quality in Jiulongjiang River basin
[J].
“红三角”地区之郴州市的可持续发展评价研究——生态足迹法
[J].
Research on sustainable development evaluation of Chenzhou City in the Red Triangle area:Ecological footprint method
[J].
郴州工业场地重金属和PAHs复合污染特征及风险评价
[J].
Characteristics and risk assessment of heavy metal and PAHs composite pollution in industrial site in Chenzhou
[J].
为了郴州新世纪的明天
[J].
For the tomorrow of Chenzhou’s new century
[J].
郴州市主要河流水质评价和保护对策
[J].
Water quality evaluation and protection countermeasures of major rivers in Chenzhou City
[J].
《土地利用现状分类》国家标准的解读
[J].
Interpretation of the national standard for classification of land use status
[J].
Influences of the land use pattern on water quality in low-order streams of the Dongjiang River basin,China:A multi-scale analysis
[J].
Water quality in the upper Han River basin,China:The impacts of land use/land cover in riparian buffer zone
[J].
The impact of urban patterns on aquatic ecosystems:An empirical analysis in Puget lowland sub-basins
[J].DOI:10.1016/j.landurbplan.2006.08.001 URL [本文引用: 1]
不同尺度土地利用方式对鄱阳湖湿地水质的影响
[J].
Effects of different scales of land use patterns on water quality in Poyang Lake wetlands
[J].
关于用SPSS中单样本K-S检验法进行正态分布等的一致性检验时适用条件的研究
[J].
Study on the applicable conditions for the consistency test of normal distribution and the like by single-sample K-S test in SPSS
[J].
祁连山东部冰沟河流域景观格局与河流水化学特征关系
[J].以2016年9月冰沟河流域的土地利用数据和河流水化学数据为基础,以河水采样点为中心分别建立三种尺度的缓冲区,运用Fragstats 4.2软件计算景观指数,再采用冗余分析方法分析不同尺度下的景观格局与水化学特征关系。结果表明:流域内不同土地利用/覆被类型对河流水化学特征的影响存在差异,耕地和建设用地的增加会在一定程度上促使主要水化学元素向河水中汇集,草地对水化学元素起到固持作用,水域和荒漠具有稀释作用,上游地区林地与水化学特征主要表现为正相关关系;该流域内,斑块密度、景观分离度和斑块连通性指数在200 m和600 m尺度下对水化学特征具有较好的指示意义;1000 m缓冲区尺度下的面积加权平均斑块分维数及600 m尺度下的景观形状指数与水化学特征关系更符合客观规律;香农多样性指数在中下游600 m缓冲区尺度下对水化学特征具有更高的解释度。本研究结果揭示了祁连山区生态水文过程,为土地的合理利用和景观格局的优化提供了科学依据。
Relationship between landscape pattern and river water chemistry characteristics of Binggou River basin in the eastern Qilian Mountains
[J].
不同空间尺度土地利用结构与景观格局对龙川江流域水质的影响
[J].
DOI:10.16258/j.cnki.1674-5906.2021.07.016
[本文引用: 1]
为了探究龙川江流域不同尺度土地利用结构与景观格局对水质的影响,于2017年5和9月对龙川江流域的21个样点进行取样,测定水温、pH、溶解氧饱和度、氮、磷等水质指标。结合2017年10 m分辨率的土地利用数据,提取3种空间尺度(河段、河岸带、子流域)的土地利用结构与景观格局指数,运用相关分析、冗余分析(RDA)、“bioenv”分析和逐步回归模型揭示不同尺度土地利用结构与景观格局对水质的影响。结果表明,龙川江水体主要污染物为TN和TP,水质基本指标有月际性差异,N、P浓度在5月高于9月,但差异不显著。不同空间尺度土地利用结构与水质指标的相关关系不同,子流域尺度土地利用组合与总体水质之间的相关性最强,且林地对水质有正面影响,耕地、建设用地对水质有负面影响。不同尺度的景观格局指数组合与水质的相关性存在差异,河岸带景观格局组合与总体水质的相关性最高,RDA分析表明景观格局指数对总体水质的总解释率为河岸带>子流域>河段。整体上,9月土地利用结构与景观格局对水质的影响要大于5月。多元回归模型发现建设用地和耕地对养分的预测能力及预测指标依赖于空间尺度和取样月份。研究以龙川江为例,证实了土地利用结构与景观格局对水质的影响具有尺度效应,水质响应的空间尺度效应应纳入土地利用规划与水资源管理的过程中,提出长江上游更多的个案研究尤为必要,对长江上游河源区水资源保护及土地利用管理具有重要的意义。
Effects of land use structure and landscape patterns at different spatial scales on water quality in Longchuan River basin
[J].
Testing a dynamic complex hypothesis in the analysis of land use impact on lake water quality
[J].DOI:10.1007/s11269-009-9498-y URL [本文引用: 1]
Influence of landscape structures on river water quality at multiple spatial scales:A case study of the Yuan River watershed,China
[J].DOI:10.1016/j.ecolind.2020.107226 URL [本文引用: 1]
流溪河流域景观空间特征与河流水质的关联分析
[J].
Correlation analysis between landscape spatial characteristics and river water quality in Liuxi River basin
[J].
Catchment versus riparian buffers:Which land use spatial scales have the greatest ability to explain water quality changes in a typical temperate watershed?
[J].
DOI:10.3390/w13131758
URL
[本文引用: 1]
Identifying the multi-scale spatial relationship between land use and water quality is critical for determining the priorities and key areas of river management. To more accurately identify the scale effect of land-use patterns on water quality and quantitatively distinguish the difference in the impacts of land-use composition and configuration on water quality, we used 94 sites to extract the upstream catchment and riparian buffer zone with different widths. The results showed that the ability of land use variables with different buffer widths to explain water quality differed slightly from the ability of these variables at the catchment scale, and the joint explanatory ability of land use composition and configuration was greater than that of each individually. The patch density and landscape shape index of cultivated land, shrubland, and built-up land in the buffer area close to the water bodies were the main factors for the increase in the concentration of total nitrogen, nitrate nitrogen, total phosphorus, and suspended solids. As the width of the buffer increased, the role of the percent of land use increased. Our research indicates that water quality management needs to adopt a multi-scale perspective and focus on key local areas while coordinating at a broader scale.
/
〈 |
|
〉 |
