城镇化流域“源-汇”景观格局对河流氮磷空间分异的影响——以天津于桥水库流域为例

张亚娟,李崇巍,*,胡蓓蓓,谢慧君,宋爱云

1 天津师范大学城市与环境科学学院, 天津 300387 2 滨州学院山东省黄河三角州生态环境重点实验室, 滨州 256603



城镇化流域“源-汇”景观格局对河流氮磷空间分异的影响
——以天津于桥水库流域为例

张亚娟1,李崇巍1,*,胡蓓蓓1,谢慧君1,宋爱云2

1 天津师范大学城市与环境科学学院, 天津 300387 2 滨州学院山东省黄河三角州生态环境重点实验室, 滨州 256603

景观格局；景观空间负荷对比指数；城镇化流域；非点源污染

随着社会经济的快速发展,流域城镇化的趋势加快[1- 3],人为活动加剧引发了诸多的河流生态环境问题[4- 5]。各种点源、非点源污染物以不同的方式随水体在各类型景观中迁移转化,导致流域景观组成、空间格局、不同空间尺度对河流水质均有影响[6- 7]。景观组成和空间结构指数缺乏一定的生态学意义,而基于“源-汇”过程的景观空间负荷对比指数(LWLI),融合了景观空间结构和景观性质,可以定量评价景观格局对污染过程的影响[8],该指数已在长江流域[9]、海河流域[10]和太湖流域[11]等多个流域进行了验证与应用。近几年来,该方法应用更加深入,结合不同区域特点,将栅格化LWLI应用于河口区域[12- 13],将基于景观坡度指数和水文响应单元的空间负荷对比指数应用于水土侵蚀区域[14]。目前,城镇化引发的污染形式多样,且随机性、复杂性、不确定性更强,而单纯的对流域农业面源污染引发的水体污染进行分析已不能满足解决城市化带来的水环境恶化问题,对城镇化流域景观格局与水体污染关系的研究已成为近年来关注的焦点[15-17]。

流域景观格局具有高度的空间变异性和空间关联性,而水体污染物的空间异质性与之密切相关。景观格局特征的空间自相关分析,可以揭示景观格局空间聚集特征,有助于理解污染物的空间分布特征[18]。目前,景观格局与水质关系的研究多采用模型研究和数理统计分析方法,基于生态水文过程的模型机理性强,模型结构较为复杂,研究需要大量基础空间数据和长时间序列的监测数据[19-20]。而时间尺度较短的监测研究多采用相关分析、回归分析等数理统计方法来探讨景观格局与水体污染的关系。由于流域污染物具有复杂性和空间异质性,景观格局与流域水质存在一定的不确定性。近年来,结合多元线性回归分析、典范对应分析和冗余分析等研究多变量关系的方法,对影响水质的多因素关系进行了较为深入的研究[21-22]。

于桥水库是引滦入津工程中的重要调蓄水库,其水环境状况直接影响着天津市工农业用水和生活饮用水。近年来流域城镇化趋势明显,居民建设用地从1984年的11.7%增加到2013年的18.1%,城镇居民点空间异质性显著增加,由平原向山区扩展明显[23]；流域的污染形式更为复杂,氮、磷等营养盐严重富集,蓝藻水华频繁暴发,这对天津市饮用水源地的安全构成了威胁[24]。对该流域水环境的相关研究从20世纪80年代一直延续至今,过去对流域中小尺度的研究较多,而对整个流域尺度的研究相对较少。陈利顶等[25]在本流域研究中应用“源-汇”景观格局理论,对流域农业面源污染进行了系统的研究,而对流域城镇化背景下,“源-汇”格局特征与流域水质空间耦合关系的研究相对较少[26]。

总之,流域城镇化过程中,河流水体污染的复杂不确定性更强。在多空间尺度上,揭示“源-汇”景观格局与河流水质空间异质关系,对流域景观格局优化具有重要意义。本研究选取于桥水库流域为研究区,基于3年来流域实地监测数据,采用空间分析及冗余分析的方法,探讨城镇化流域平水期与丰水期景观格局特征对水质的影响,为于桥水库流域的水环境保护与管理提供科学依据。

1 研究区域

于桥水库流域(39°56′—40°23′N,117°26′—118°12′E)总面积约2060km2,位于天津市北部蓟县城东4km处,其中约80%位于河北省遵化市地区。地形地势从东北向西南呈阶梯下降,西北部是山区,中部为平原,南部为丘陵和低山区。流域属温带大陆季风型半湿润气候区,年平均降水量在700 mm以上[24]。汇入该流域的主要河流有沙河、淋河和黎河。流域附近广泛分布着村镇,产生的各种污染物直接或间接影响着流域的水质。为了能全面的反映流域内水环境特征,本研究在全流域内划分出33个子流域开展水质监测与采样分析(图1)。

图1 于桥水库流域地理位置及采样点分布Fig.1 Geographical position and sampling points in the Yuqiao Reservoir watershed

2 研究方法

2.1 土地利用类型的划分

利用2013年8月美国陆地卫星Landsat8遥感影像解译获得于桥水库流域土地利用数据,在ENVI 5.1支持下对其进行大气辐射校正、几何校正等预处理,对解译结果进行实地调查校正,解译精度为87%,将流域土地利用类型划分为水体、林地、城镇居民用地、园地、裸地和耕地6种类型[26]。

2.2 水样点设置与测定

表1 水质数据描述性统计

2.3 景观格局指数的获取

基于研究区的实际情况及景观指数的生态学意义,本研究选取景观多样性指数(SHDI)、斑块密度(PD)、面积加权分形维数(area-weighted mean patch fractal dimension, AWMPFD)和最大斑块指数(LPI)与水质相关性较高的景观指数,用以表征流域景观格局特征[31-32]。其中AWMPFD反映人类活动对景观格局的影响程度,LPI是斑块聚集程度及优势景观类型的指数,而PD及SHDI都是景观破碎化和丰富程度的指数。采用Fragstats 4.2软件计算各种景观格局特征指数。

陈利顶、傅伯杰等[33- 34]在“源-汇”景观理论的基础上,借用洛伦兹曲线提出了源汇景观空间负荷对比指数,来定量描述景观空间格局与非点源污染之间的过程。根据前人的研究[10,29],结合土壤侵蚀通用方程中的植被覆盖与管理因子(C值),对其“源”和“汇”的贡献权重进行了赋值。各土地利用类型赋值依次为林地(0.8)、耕地(0.4)、城镇居民用地(1.0)、园地(0.4)、裸地(0.5)。基于洛伦兹曲线的景观空间负荷对比指数模型如下：

(1)

LWLI=LWLI′distance×LWLI′elevation÷LWLIslope

(2)

式中,LWLI代表综合的景观空间负荷对比指数,Asourcei和Asourcej分别指第i种“源”景观和第j种“汇”景观面积累积曲线组成的不规则三边形面积,Wi和Wj指“源”景观和“汇”景观类型的权重[10],APi和APj指第i种源景观和第j种汇景观类型的面积百分比,m和n分别为“源”景观和“汇”景观类型的总数目。

2.4 统计分析

水质数据为2013、2014和2015年每年5—10月份于桥水库流域的监测数据,分别选取每年丰水期和平水期各监测点数据进行统计分析。利用ArcGIS 10.2空间统计工具模块,对研究区源汇景观空间负荷对比指数(LWLI)进行空间自相关(Moran′sI)和聚类空间关联性(Anselin Local Moran′sI))分析[18]；子流域“源-汇”景观特征指数与子流域水质数据采用Spearman秩相关性分析在SPSS17.0中完成；逐步回归分析用来确定最能解释单个河流水质变量的环境或景观变量因子；CANOCO 4.5软件进行冗余分析(Redundancy analysia, RDA),用以揭示单个景观变量对河流水质变化的贡献。

3 研究结果

3.1 于桥水库子流域景观格局分析

33个子流域内景观类型的构成(图2)按照“源”型景观类型进行面积比例排列。流域内耕地和城镇居民用地面积比分别为7.9%和17%,“汇”型景观中林地面积占流域总面积的33.5%。城镇居民用地在各个子流域中均有分布,城镇居民用地所占子流域面积比例在0.1%—56.4%,中下游地区城镇居民用地高于上游地区,城镇居民用地的面积百分比在中下游子流域平均值达34.6%,个别子流域超过56%。受人为因素影响,耕地主要分布在中下游地区,占子流域面积比例平均值为4.6%。城镇居民用地和耕地所占面积比例从上游到下游逐渐增大,整个流域从上游到下游呈现“汇”景观减小,“源”景观增大的趋势。

城镇居民用地分布广泛且空间聚集特征明显,上游地区以林地和园地为主要土地利用类型,虽然城镇居民用地面积比例低,但分布范围很广泛,如地势较为陡峭,海拔高度为700多米的四拨子仍然分布有居民地。城镇居民用地分布主要集中于中下游遵化市周边。15号流域包含遵化市区,其城镇化面积比最大,为56.4%,其次为16、12、17、8和5号流域,主要分布在遵化市周边地区。林地面积比例最大的流域为27号,其次为33、20、31、26、32和21号流域,分布在上游流域主要集中于清东陵等山区。

图2 各子流域景观类型构成Fig.2 Composition of landscape types in different watersheds

3.2 流域景观空间负荷对比指数空间关联性及其与流域氮磷空间分布的相关性

空间自相关分析是一系列空间数据分析方法和技术的集合[18],对LWLI进行空间自相关性分析,可揭示LWLI的空间集聚特征与流域“源-汇”景观格局及流域城镇化的空间分布的内在联系,而流域城镇居民用地和汇景观的分布与河流水质关系密切。首先使用全局 Moran′sI指数用于验证整个研究区域LWLI的空间相关关系,LWLI全局Moran′sI的值为0.637,P<0.01,表明研究区域LWLI存在较强的正相关情况,其空间分布存在一定的内在联系,即LWLI在空间上存在趋于集群的现象。

全局空间自相关指标可以检验整个区域某一要素的空间分布模式,但全局Moran′sI不能用来测度相邻区域之间要素或属性的空间关联模式,因此需要讨论局部子流域LWLI的相关程度。为此,本文进一步对子流域LWLI进行局部空间关联格局分析,得出 LWLI集聚图(图3),高-高和低-低聚集模式的区域代表具有显著的空间自相关性。位于高-高聚集模式的区域主要集中在中部遵化市周边的15西留村、16张七各庄村、12水平口、8鄂庄子村、5东双城村和11堡子店,这些子流域的城镇居民用地面积比分别为56.4%、50.2%、43.7%、38.5%、31.2%、11.9%和6.8%,均属于“源”型景观起主导作用城镇化较为突出的区域,LWLI高-高区域与城镇化集中区域具有一致性,表明LWLI能较好地反映流域城镇化趋势。低-低聚集模式主要分布在西北部地区淋河流域的28大鹿圈、20金山子村、31后龙台、26小溪、21四拨子和32二拨子流域,主要以灌草地和林地为主,居民建设用地仅有0.08%、0.21%、0.25%、0.38%、0.44%和0.83%,林地占79.5%—92.5%,人为干扰较少,“汇”型景观在流域内起主导作用。LWLI的空间分布与流域“源-汇”景观格局形成对应关系,在空间分布上与城镇化集中区域具有一致性。

图3 源汇景观空间负荷对比指数局部空间自相关图Fig.3 The LWLI of local spatial autocorrelation in the Yuqiao Reservoir watershed

表2 LWLI与水质指标之间的Spearman相关性分析

**表示P<0.01；LWLI：景观空间负荷对比指数

3.3 水质变化的主导景观变量辨析

为了更直观反应景观特征指数与水质指标之间的关系,选取RDA线性模型,以33个检测点的水质数据和10个景观特征指数为样本对其进行冗余分析(图4)。当景观特征指数与水质指标之间的箭头夹角小于90度,两者之间为正相关；夹角等于90度,两者不相关；夹角大于90度,两者之间为负相关。由RDA排序图可看出,丰水期与平水期水质指标与城镇居民用地比例(Red%)、裸地比例(Bar%)、园地比例(Ord%)、耕地比例(Cul%)、斑块密度(PD)、景观多样性指数(SHDI)和面积加权分形维数(AWMPFD)呈正相关,与林地比例(For%)、最大斑块指数(LPI)呈负相关。另外,箭头的长度反应该景观特征指数对水质指标的影响程度,箭头越长影响程度越大,可看出土地利用类型中Red%和For%对水质指标的影响最为显著,与Bar%和Cul%的相关性并不显著。

For%、Ord%、LPI、AWMPFD与上游地区的子流域主要分布在第二排序轴的左侧,结合子流域景观类型构成可知,这些子流域林地面积百分比相对较高,受人为干扰较小。分布在排序图右侧的子流域主要受SHDI、PD、Red%、Bar%、Cul%的影响。结合子流域景观类型构成分析,这些子流域城镇居民用地面积比例较高,且景观破碎化严重,人为干扰强烈,其中15、8、16、12、5和11号子流域对氮、磷影响较大,这与城镇居民用地面积百分比大小具有一致性,表明城镇居民用地是流域水质污染的重要贡献源。

图4 不同水期的RDA排序Fig.4 RDA ordination in different water periods

4 讨论

4.1 景观特征指数与水质之间的关系

许多研究表明城镇居民用地是水质退化最主要的贡献源,森林能够有效控制水体中污染物的浓度[2,27-28],本研究中RDA排序图直观展示了平水期与丰水期城镇居民用地与水质指标之间的相关性均为最高,由逐步回归分析和Spearman相关分析可知,林地与水质指标呈显著负相关,说明“汇”景观对水体的净化功能较强,即水质较好的流域具有较高的森林覆盖率和较低的城镇化水平。由于城镇居民用地遍及所有的子流域,居民生活区的污染物经不透水面排入河道更为便捷,加剧了水质的恶化。位于流域中下游地区的遵化市周边城镇居民用地面积比例较大,河道水质相对较差,城镇化较高,人为干扰严重,除企业排污及矿业废水等点源污染外,居民生活污水及家禽养殖对流域面源污染贡献较大。

表3 景观特征指数与水质指标之间的逐步回归分析

For%：林地比例forest percent；Red%：城镇居民用地比例residential percent；Ord%：园地比例Orchard percent；SHDI：景观多样性指数shannon′s diversity index

由RDA排序结果可知,斑块密度(PD)、景观多样性指数(SHDI)和面积加权分形维数(AWMPFD)与水质指标呈正相关,PD、AWMPFD和SHID的增高意味着景观的破碎化程度越高,受人类干扰强度越大,越不利于控制污染物的迁移和转化；最大斑块指数(LPI)与水质指标呈负相关,高的LPI意味着受人类干扰强度越小,水质往往较好；城镇化引发的污染随机性、复杂性、不确定性更强,水体中的污染物输出只能被少数景观指数解释,甚至有些水质指标与景观格局指数之间的关系无法被合理解释．多数景观指数是对景观格局几何特征的描述和分析,忽略了对景观格局特征的内涵理解[36]。在“源-汇”理论基础上建立的LWLI,能够较好的融合景观类型、面积、空间位置和地形特征[8],在本研究中和氮磷指标有显著的相关性,有效地刻画水体污染的空间异质性。因此,适宜的景观格局指数能有效说明流域水质的空间分布特征。

4.2 LWLI与水质之间的关系

在逐步回归分析中,LWLI均为解释变量,表明LWLI能有效地反映氮、磷输出特征,结合Spearman相关性分析结果(表1),表明LWLI与流域氮、磷空间分布存在极显著的相关性,这与已有研究结果是相似的[9- 10]。本研究中(LWLI)与流域氮、磷的相关性要高于其他景观特征指数,这与孙然好等[10]和李晶等[14]的观点相同。LWLI指数是基于高程、坡度、距离和景观类型的综合指标[10],较好地融合了景观组成结构、空间结构和景观性质,可以定量评价水体的污染程度。从平水期到丰水期LWLI与磷污染相关性增加,而与氮污染相关性较低,可能由于LWLI中体现了地形因素的影响,河流中水质的空间分布与地形因素密切相关,研究表明强降雨作用下,地势较高将引发较高的侵蚀率,使得土壤中颗粒态的磷大量的流入河流中,同时流域内园地施加的磷肥在丰水期也会被雨水冲出进入河道,导致河流中的磷浓度增加。LWLI涵盖地形因子和景观类型因子,在汛期LWLI能更好的解释磷污染的状况。氮含量与人为因素密切相关,降水较少的情况下,居民生活区是氮输出的主要形式[36],本流域内城镇居民用地遍布所有子流域,同时城镇居民用地与水质指标之间的相关性最显著。由于LWLI考虑了景观类型和汇水口距离的影响,因此,在平水期相对传统景观指数,LWLI能更好地解释氮的输出。

5 结论

(1)整个流域从上游到下游呈现“汇”景观减小,“源”景观增大的趋势。城镇居民用地主要分布在中下游地区遵化市周边,城镇居民用地的面积百分比在中下游子流域平均值达34.6%,个别子流域超过56%。城镇用地和耕地所占面积比从上游到下游逐渐增大,林地所占面积比例呈现相反的趋势。

(2)LWLI全局Moran′sI的值为0.637,P<0.01,在空间上存在趋于集群的现象。高-高聚集模式的区域主要集中在中部遵化市区周边,这与城镇居民点主要集中分布于遵化市周边一致。低-低聚集模式主要分布在流域上游西北部地区,而林地等“汇”型景观也主要集中在清东陵等山区。

在流域城镇化发展中,城镇居民用地分布广泛,源汇景观破碎化程度较高。建议优化流域“源-汇”景观格局。控制“源”景观的排放,提高村镇景观的连通性,有利于污水和垃圾的收集和处理；优化农林格局,推广生态农业,减少化肥农药的使用。同时,提高“汇”景观的调节功能,退耕还林,降低高山区林地的景观破碎化；流域下游应增加林地、草地的斑块,可以改善水质优化环境。

[1] Xian G, Crane M, Su J S. An analysis of urban development and its environmental impact on the Tampa Bay watershed. Journal of Environmental Management, 2007, 85(4): 965- 976.

[2] Evans D M, Schoenholtz S H, Wigington Jr P J, Griffith S M, Floyd W C. Spatial and temporal patterns of dissolved nitrogen and phosphorus in surface waters of a multil and use basin. Environmental Monitoring and Assessment, 2014, 186(2): 873- 887.

[3] 李静芝, 朱翔, 李景保, 徐美. 洞庭湖区城镇化进程与水资源利用的关系. 应用生态学报, 2013, 24(6): 1677- 1685.

[4] Shen Z Y, Houa X S, Li W, Aini G, Chen L, Gong Y W. Impact of landscape pattern at multiple spatial scales on water quality: A case study in a typical urbanised watershed in China. Ecological Indicators, 2015, 48: 417- 427.

[5] Sun R H, Chen L D, Chen W L, Ji Y H. Effect of Land-Use Patterns on Total Nitrogen Concentration in the Upstream Regions of the Haihe River Basin, China. Environmental Management, 2013, 51(1): 45- 58.

[6] Bu H M, Meng W, Zhang Y, Wan J. Relationships between land use patterns and water quality in the Taizi River basin, China. Ecological Indicators, 2014, 41: 187- 197.

[7] 李艳利, 李艳粉, 徐宗学, 关中美. 浑太河上游流域河岸缓冲区景观格局对水质的影响. 生态与农村环境学报, 2015, 31(1): 59- 68.

[8] 陈利顶, 李秀珍, 傅伯杰, 肖笃宁, 赵文武. 中国景观生态学发展历程与未来研究重点. 生态学报, 2014, 34(12): 3129- 3141.

[9] 刘芳, 沈珍瑶, 刘瑞民. 基于“源-汇”生态过程的长江上游农业非点源污染. 生态学报, 2009, 29(6): 3271- 3277.

[10] 孙然好, 陈利顶, 王伟, 王赵明. 基于“源”“汇”景观格局指数的海河流域总氮流失评价. 环境科学, 2012, 33(6): 1784- 1788.

[11] 王瑛, 张建锋, 陈光才, 单奇华, 李泽波, 徐永辉, 陈云. 太湖流域典型入湖港口景观格局对河流水质的影响. 生态学报, 2012, 32(20): 6422- 6430.

[12] Jiang M Z, Chen H Y, Chen Q H. A method to analyze “source-sink” structure of non-point source pollution based on remote sensing technology. Environmental Pollution, 2013, 182: 135- 140.

[13] Jiang M Z, Chen H Y, Chen Q H, Wu H Y. Study of landscape of patterns variation and optimization based on non-point source pollution control in an estuary. Marine Pollution Bulletin, 2014, 87(1/2): 88- 97.

[14] Li J, Zhou Z X. Coupled analysis on landscape pattern and hydrological processes in Yanhe watershed of China. Science of the Total Environment, 2015, 505: 927-938.

[15] 冯源嵩, 林陶, 杨庆媛. 百花湖周边城市近郊小流域氮、磷输出时空特征. 环境科学, 2014, 35(12): 4537- 4543.

[16] 李昆, 王玲, 李兆华, 王祥荣, 陈红兵, 吴忠, 朱鹏. 丰水期洪湖水质空间变异特征及驱动力分析. 环境科学, 2015, 36(4): 1285- 1292.

[17] Shen Z Y, Hou X S, Li W, Aini G. Relating landscape characteristics to non-point source pollution in atypical urbanized watershed in the municipality of Beijing. Landscape and Urban Planning, 2014, 123: 96-107.

[18] 孙才志, 闫晓露, 钟敬秋. 下辽河平原景观格局脆弱性及空间关联格局. 生态学报, 2014, 34(2): 247- 257.

[19] Chiang L C, Yuan Y P, Mehaffey M, Jackson M, Chaubey I. Assessing SWAT′s performance in the Kaskaskia River watershed as influenced by the number of calibration stations used. Hydrological Processes, 2014, 28 (3): 676-687.

[20] Ouyang W, Song K Y, Wang X L, Hao F H. Non-point source pollution dynamics under long-term agricultural development and relationship with landscape dynamics. Ecological Indicators, 2014, 45: 579- 589.

[21] 赵鹏, 夏北成, 秦建桥, 赵华荣. 流域景观格局与河流水质的多变量相关分析. 生态学报, 2012, 32(8): 2331- 2341.

[22] 刘丽娟, 李小玉, 何兴元. 流域尺度上的景观格局与河流水质关系研究进展. 生态学报, 2011, 31(19): 5460- 5465.

[23] 陈萍, 李崇巍, 王中良, 孟伟庆, 包姗姗, 谢慧君. 天津于桥水库流域典型景观类型时空演变分析. 生态学杂志, 2015, 34(1): 227- 236.

[24] 徐媛, 谢汝芹, 卢蔚, 王玉秋. 于桥水库富营养化评价及空间分布特征研究. 水资源与水工程学报, 2014, 25(1): 1- 6.

[25] 陈利顶, 傅伯杰, 张淑荣, 丘君, 郭旭东, 杨福林. 异质景观中非点源污染动态变化比较研究. 生态学报, 2002, 22(6): 808- 816.

[26] 包姗姗, 李崇巍, 王祖伟, 陈萍, 郭继发. 天津于桥水库流域水体氮磷空间分异与景观格局的关系. 农业环境科学学报, 2014, 33(8): 1609- 1616.

[27] 王琼, 姜德娟, 于靖, 张华. 小清河流域氮磷时空特征及影响因素的空间与多元统计分析. 生态与农村环境学报, 2015, 31(2): 137- 145.

[28] 王欢, 袁旭音, 陈海龙, 许海燕, 李正阳. 太湖流域上游西苕溪支流的营养状态特征及成因分析. 湖泊科学, 2015, 27(2): 208- 215.

[29] Chen L D, Tian H Y, Fu B J, Zhao X F. Development of a New Index for Integrating Landscape Patterns with Ecological Processes at Watershed Scale. Chinese Geographical Science, 2009, 19(1): 37- 45.

[30] 国家环境保护总局. 水和废水监测分析方法(第四版). 北京: 中国环境科学出版社, 2002.

[31] 焦胜, 杨娜, 彭楷, 郭谌达, 李振民, 周怀宇. 沩水流域土地景观格局对河流水质的影响. 地理研究, 2014, 33(12): 2263- 2274.

[32] 欧洋, 王晓燕, 耿润哲. 密云水库上游流域不同尺度景观特征对水质的影响. 环境科学学报, 2012, 32(5): 1219- 1226.

[33] 陈利顶, 傅伯杰, 徐建英, 巩杰. 基于“源-汇”生态过程的景观格局识别方法—景观空间负荷对比指数. 生态学报, 2003, 23(11): 2406- 2413．

[34] 陈利顶, 傅伯杰, 赵文武. “源”“汇”景观理论及其生态学意义. 生态学报, 2006, 26(5): 1444- 1449.

[35] 吉冬青, 文雅, 魏建兵, 吴志峰, 刘庆, 程炯. 流溪河流域景观空间特征与河流水质的关联分析. 生态学报, 2015, 35(2): 246- 253.

[36] 韩黎阳, 黄志霖, 肖文发, 田耀武, 曾立雄, 吴东. 三峡库区兰陵溪小流域土地利用及景观格局对氮磷输出的影响. 环境科学, 2014, 35(3): 1091- 1097.

Impact of a “source-sink” landscape pattern in an urbanized watershed on nitrogen and phosphorus spatial variations in rivers: a case study of Yuqiao Reservoir watershed, Tianjin, China

ZHANG Yajuan1, LI Chongwei1,*, HU Beibei1, XIE Huijun1, SONG Aiyun2

1CollegeofUrbanandEnvironmentalScience,TianjinNormalUniversity,Tianjin300387,China2ShandongProvinceKeyLaboratoryofEco-environmentalScienceforYellowRiverDelta,BinzhouUniversity,Binzhou256603,China

landscape pattern; location-weighted landscape contrast index; urbanized watershed; non-point source pollution

10.5846/stxb201512242566

国家自然科学基金项目(31270510, 41301221)；滨州学院博士基金项目(2014Y02,2014Y14)

2015- 12- 24; 网络出版日期:2016- 08- 30

张亚娟,李崇巍,胡蓓蓓,谢慧君,宋爱云.城镇化流域“源-汇”景观格局对河流氮磷空间分异的影响——以天津于桥水库流域为例.生态学报,2017,37(7):2437- 2446.

Zhang Y J, Li C W, Hu B B, Xie H J, Song A Y.Impact of a “source-sink” landscape pattern in an urbanized watershed on nitrogen and phosphorus spatial variations in rivers: a case study of Yuqiao Reservoir watershed, Tianjin, China.Acta Ecologica Sinica,2017,37(7):2437- 2446.

*通讯作者Corresponding author.E-mail: yllcw13@126.com