基于生态系统服务的流域分区规划与管理——以于桥水库流域为例

张海洋，杨琰瑛，刘凯，周滨，苏德岳，马成，师荣光

（1.农业农村部环境保护科研监测所，天津 300191；2.中国农业科学院研究生院，北京 100081；3.农业农村部农产品质量安全环境因子控制重点实验室，天津 300191；4.天津市生态环境科学研究院，天津 300191；5.天津市于桥水库管理中心，天津 300074）

生态系统服务是指人类直接或间接从自然生态系统中获得的各种物质或精神上的惠益[1-2]。随着社会经济发展速度的提高和人类活动强度的增加，流域生态系统受到的影响越来越大，生态环境逐渐恶化[3]。对生态系统服务进行有效的评价与管理，将有助于生态系统的可持续发展[4-5]。量化生态系统服务的时空分布，能够明确流域生态系统服务的空间分布特点，优化空间布局，常用的评价方法有当量因子法[6]、能值转化法[7]和生态模型法[8]等。其中，生态模型法随着近些年3S 技术的发展而得到了广泛的应用。生态模型能够定量评估流域生态功能时空分布特征，分析多种类型因素变化对生态系统服务功能的影响机制，对生态系统服务演变进行分析与预测[9-10]，突破传统价值评估方法无法表征生态功能空间分布情况的局限性[11]。

生态系统服务的多样性和空间分布的差异性，加之人类活动的影响，很容易造成区域生态系统的不稳定、不平衡，因此有必要进行合理的空间规划，以提高资源利用率，支撑社会经济发展，促进生态文明建设[12]。1976 年，美国生态学家Bailey[13]首次提出了基于生态系统角度的区划，此后，关于生态系统服务区域划分的研究逐渐增多[14]。我国在20 世纪末才开始对生态系统服务功能分区研究给予重视[15]。生态系统中往往存在着多种服务之间的权衡或协同关系[16]，一种生态系统服务的增强会导致另一种服务的增强或减弱[17]。如粮食产量与水质净化间会存在较强的权衡关系，因为粮食生产过程中肥料使用量的增加会导致土壤中更多营养物质随地表径流进入水体[18]；林地面积的增加则会提升流域的固碳量，但同时也会减少径流的产生，即产水量下降[19]。生态系统服务簇能够识别生态系统中不同分区的主导服务，它是指一系列在时间和空间上不断重复出现的生态系统服务[20]，可以为管理不同类别的区域提供足够的权衡信息。目前多采用K-means 聚类分析法[21]、自组织特征映射网络法（SOFM）[22]、主成分分析法（PCA）[23]和空间自相关分析法[24]等识别生态系统服务簇，以便对研究区进行生态功能分区。学者们在大量的深入研究下，取得了丰富成果[25-27]，Gong 等[28]利用K-means 聚类方法将甘肃白龙江流域划分为5 个生态功能分区；闫晓露等[29]使用自组织神经网络方法将大连市划分为4 类服务簇；Lin 等[30]借助K-means 聚类分析对西南三江流域进行了生态系统服务簇的识别；李慧蕾等[31]使用自组织特征映射网络对内蒙古自治区进行了生态功能区划分；生态系统服务间的权衡协同关系与服务簇识别存在尺度效应[21]，然而由于小尺度上数据获取难度较大，当下所研究的多为行政区域和大流域尺度[32-35]，小流域尺度的划分研究较为有限。因此本研究尝试在小流域尺度上对不同生态系统服务的空间分布及权衡/协同特征进行研究，在此基础上对流域的生态功能区进行划分，进而加强流域的生态管理。

于桥水库是引滦入津工程中重要的调蓄水库，是国家级重要饮用水水源地，也是南水北调前天津市唯一的饮用水源地，其生态环境状况备受关注。有研究显示，流域内存在着较为严重的面源污染问题，对水质及生态环境产生了一定的影响，并且农业发展与水质保护关系极不平衡[36-37]。近些年，对该流域的研究多为水质[38-39]和景观格局[40-41]等角度的分析，尚缺少对整体流域生态系统服务功能角度的分区管理探究。

鉴于此，本研究以于桥水库流域为研究对象，利用InVEST模型对整个流域的粮食产量、氮输出、磷输出、固碳量、产水量和土壤保持6 项生态系统服务功能进行评估，分析各项生态系统服务的空间分布差异，以此为基础探究生态系统服务之间的权衡协同关系，并利用生态系统服务簇对研究区进行生态功能区划分，以期为流域规划和管理提供决策依据。

1 研究区概况

于桥水库流域（39°56'～40°23'N，117°26'～118°12'E）位于天津市北部与河北省交界处，流域总面积约2 060 km2，包括淋河、沙河和黎河3条主要河流，是京津冀的重要生态屏障（图1）。其流域面积位于河北省的面积约占78%，涉及唐山市遵化市、玉田县和迁西县，承德市兴隆县和天津市蓟州区。流域地形自西北向东南倾斜，西北和北部为山地地区，中部为平原区，南部和东南地区为低山丘陵区。土地利用类型以耕地、林地为主。农业是该流域主要的支柱产业，但流域整体经济发展水平不高，是典型的农业流域。气候类型为温带大陆季风半湿润气候，年平均气温10～13 ℃，年平均降雨655.5 mm。土壤类型主要为褐土，植被类型多种多样，其中灌木林分布最多。

图1 研究区地理位置及土地利用图Figure 1 Geographical location and land use of the study area

2 材料与方法

2.1 数据来源

本研究采用2020 年土地利用遥感监测数据，精度为30 m×30 m，来源于中国科学院资源环境科学与数据中心（https：//www.resdc.cn/）的中国土地利用遥感监测数据库；数字高程（DEM）数据来源于Quick Bird卫星，精度为37.5 m×37.5 m。土壤最大根系深度数据和土壤有机质含量数据来源于国家地球科学数据中心-土壤分中心的二普土种剖面数据[42]。土壤质地数据与气象数据均来源于中国科学院资源环境科学与数据中心（https：//www.resdc.cn/）。

2.2 研究方法

2.2.1 生态系统服务评估

本研究利用于桥水库流域数字高程（DEM）数据，借助GIS Hydrology 功能提取得到流域边界和79个子流域，并根据以往的研究，将河网流量阈值调整为10 000，为避免模型运行时间过长，通过编辑器功能，结合河网分布情况，最终将子流域合并为26 个。本研究利用InVEST模型对研究区进行生态系统服务评估，基于2020 年土地利用数据，划分出26 个子流域，并选取Water yield、Sediment delivery ratio（SDR）、Nutrient delivery ratio（NDR）、Crop production（percentile）和Carbon 5 个模块，分别对于桥水库流域内产水量、土壤保持、氮输出、磷输出、粮食产量和固碳量6种生态系统服务进行模拟。

（1）产水量

Water yield 模块是一个简化的水循环模型，基础公式为：

式中：Yxj为第j类土地利用类型栅格x的产水量，mm；AETxj为第j类土地利用类型栅格x的每年实际水分蒸发量，mm；Pxj为第j类土地利用类型栅格x的年降雨量，mm。

（2）水土保持

SDR 模块的原理是模拟陆地表面沉积物的产生以及沉积物向河流的运移，对传统通用土壤流失方程进行修正后，公式[43]为：

式中：Rx为降雨侵蚀力因子；Kx为土壤可侵蚀性因子；LSx为坡度坡长因子；Cx为覆被管理因子；Px为措施因子。Cx和Px参考以往研究[44-47]进行整理。

（3）水质净化

NDR 模块用于估算一定区域内氮、磷的输出量，原理为分别计算研究区地上部分和地下部分氮和磷的输出量，两部分之和即为氮、磷的输出总量。氮、磷的输出量越大，水质净化功能就越弱。

地上部分养分运移率计算公式为：

式中：NDR0,i是i栅格中未被下游栅格所截留的养分比例；ICi是栅格i的地形指数；IC0和k为标定参数。

地下部分养分运移率计算公式：

其中：effsubs是地下所能达到的最大养分拦截效率；lsubs是地下径流的截留长度，m；li是栅格i到径流的距离，m。

最终的养分运移总量为地上运移量和地下运移量之和。其中，栅格i处的养分运移量为：

式中：Loadsurf，i是地上部分养分负荷量；Loadsubs，i是地下部分养分负荷量；NDRsurf，i是地上部分养分运移率；NDRsubs，i是地下部分养分运移率。

（4）粮食产量

Crop Production（percentile）模块根据现有的土地利用数据，并结合数据库中已有的175种作物在2000年观测的产量，以气候作为驱动去估算作物产量，并计算集约化分别提升25%、50%、75%和95%情况下的作物产量。Crop Production（percentile）模块工作原理为在数据库中裁剪出研究区的气候地图，通过已裁剪出的气候地图进行重分类，计算研究区集约化提升水平下的产量，并将研究区非耕地部分的作物产量值赋为0。

前期调研发现，于桥水库流域以玉米为主要粮食作物，因此用玉米产量表征流域粮食产量。通过与统计年鉴中流域内实际玉米产量进行比较，选择集约化提升25%水平下的玉米产量作为模拟结果。

（5）固碳

Carbon 模块基于规则单元的栅格土地利用图运行，每个单元分配了一种土地利用类型。模块需要每个土地利用类型的四种基本碳库（地上部分碳储量、地下部分碳储量、土壤碳储量和死亡有机物碳储量），四个碳库之和即为总碳储量，公式[48]为：

式中：Ctotal为总碳储量，t·hm-2；Cabove为地上部分碳储量，t·hm-2；Cbelow为地下部分碳储量，t·hm-2；Csoil为土壤碳储量，t·hm-2；Cdead为死亡有机物碳储量，t·hm-2。

由于死亡有机物碳储量在碳库中占比较小，故本研究不予考虑，统一设置为0。其他三大碳库参考以往研究[49-52]进行设置。

2.2.2 生态系统服务权衡协同关系

本研究以子流域为分析单元，基于流域中26 个子流域的6 种生态系统服务单位面积服务量进行相关性分析，利用SPSS 软件中的Spearman 相关系数法来衡量生态系统服务间的权衡与协同关系[53]。当相关系数r＜0 时，生态系统服务间为权衡关系；反之则为协同关系。

2.2.3 生态系统服务功能分区

利用主成分分析和K-means 聚类分析法进行生态系统服务功能区的划分。本研究中的主成分分析以26 个子流域的6 种生态系统服务总量进行统计分析。主成分分析可以去除变量中的冗余信息[54]，判断各生态系统服务在生态系统服务簇中共存关系，并对生态系统服务间的重复供给情况进行排除，使得聚类的稳健性提高[55]。K-means 聚类分析法能够识别具有相似生态系统服务的组合。

3 结果与分析

3.1 生态系统服务空间分布格局

于桥水库流域生态系统服务的空间分布具有明显差异（图2）。产水量服务整体呈现北高南低的分布格局，总产水量为2.28×108m3，平均产水深度111.55 mm。高值区主要分布在河北省遵化市城区和流域北部的连片山区以及水库西北部局部地区。三条主要河流中沙河和黎河上游产水量较大，产水深度在300～400 mm之间，淋河上游的产水量较小，产水深度低于50 mm。1～3 号、8 号子流域产水量较高，7 号、20 号和25 号子流域产水量较低（表1）。土壤保持服务表现出西北高、东南低的分布趋势，土壤保持量较高的区域集中在西北部山区，即天津市蓟州区山区和河北省兴隆县山区，整个流域的土壤保持量为1.75×107t。高值区主要包括1 号、8 号、11 号子流域，低值区包括17号、20号、26号子流域。流域氮输出总量为1 978.16 t，输出量在流域中部和南部平原地区较高，并且在高值区内呈现由东北部平原地区向西南部水库区逐渐减少的趋势。流域四周的环中部平原山区、丘陵地带和水库库区氮输出较低。氮输出最高的是3号子流域，输出最低的为20号子流域。流域磷输出总量为30.03 t，其空间分布与氮输出极为相似，高值区亦为流域中部及南部平原区，说明中部和南部平原区水质净化服务较差，需重点关注。流域周边的山地丘陵区与水库库区为低值区，磷输出最高的是3 号子流域，最低的为20 号子流域。粮食产量服务的空间分布与氮磷输出服务相似，高值区分布在流域中部和南部平原区，两种服务均与耕地存在密切关系，而流域中部和南部主要土地利用类型即为耕地。于桥水库流域粮食（玉米）产量为2.38×105t。产量较高的为10 号、22 号子流域，较低的为20 号、26 号子流域。固碳服务整体呈现西北高、东南低的分布趋势，总固碳量达1.80×107t，西北部的固碳强度最高，其次为北部，流域中东部以及水库库区为低值区。固碳服务的高值区与土壤保持服务的高值区有较高重合，基本分布在流域中的山区，而流域中山区土地利用类型多为林地，因此固碳能力与土壤保持能力较强。固碳量较高的为8 号、11 号子流域，较低的为20 号、26 号子流域。

表1 于桥水库流域子流域生态系统服务量Table 1 Ecosystem services in sub-watersheds of Yuqiao Reservoir watershed

图2 于桥水库流域生态系统服务空间分布Figure 2 Spatial distribution of ecosystem services in Yuqiao Reservoir watershed

3.2 生态系统服务的权衡协同关系

相关性分析（表2）显示，氮输出与磷输出之间相关系数为0.869（P＜0.01），呈显著正相关，互为协同关系；氮输出与土壤保持相关系数为-0.497（P＜0.01），呈显著负相关，磷输出与土壤保持并未表现出显著相关性，但存在一定的负相关关系，因此，可以说明水质净化服务与土壤保持服务互为协同关系；粮食产量与氮输出和磷输出的相关系数分别为0.963（P＜0.01）和0.795（P＜0.01），皆呈显著正相关，表明粮食产量服务与水质净化服务间存在较强的权衡关系；土壤保持与粮食产量相关系数为-0.490（P＜0.05），表明土壤保持服务与粮食供给服务为权衡关系；土壤保持与固碳量相关系数为0.839（P＜0.01），呈显著正相关，说明土壤保持服务与固碳服务为协同关系。

表2 生态系统服务间权衡协同关系Table 2 Trade-offs and synergies among ecosystem services

3.3 生态系统服务功能分区

主成分分析筛选出两个主成分因子（表3），主成分1 解释了60.42%的方差变异，粮食产量、产水量、氮输出和磷输出具有较高的正荷载，由于氮、磷输出与粮食产量及产水量间均呈正相关关系（表2），因此主成分1 反映出水质净化服务与产水量及粮食产量服务间是权衡关系，但氮和磷是以输出量来表征服务的功能，所以指标间关系为协同关系，能够在同一生态系统服务簇中共存；主成分2 解释了26.25%的方差变异，土壤保持服务与固碳服务具有较高的载荷，所以主成分2 反映的是二者间的协同关系，且二者呈正相关关系，因此能够在同一生态系统服务簇中共存。

表3 生态系统服务因子载荷Table 3 Factors loading of ecosystem services

通过K-means 聚类分析法对于桥水库流域进行生态系统服务集划分，并根据R语言中的聚类方差解释度（TWSS）[56]以及流域内的生态特征，结合不同服务集内提供的生态系统服务的差异性，最终划分为水源保护区、森林生态区、粮食生产区和综合发展区4个分区，并用玫瑰图表示不同分区内生态系统服务间权衡与协同关系（图3）。

图3 生态系统服务功能分区及生态系统服务簇特征Figure 3 Ecosystem services function zones and their ecosystem services bundles

3.3.1 水源保护区

水源保护区是以于桥水库库区（26 号子流域）为核心区域，由周边的15 号、17 号、20 号、21 号、24 号和25号子流域组成（图3），总面积382.85 km2，占流域面积的18.70%。整个区域位于流域下游地区，行政区划上大部分位于天津市蓟州区内。区域内总体上以产水服务和粮食供给服务为主，其他服务水平均较低。氮输出高于磷输出，说明在该分区内氮元素对库区水质的影响大于磷元素。该分区内土地利用类型以水域和耕地为主，耕地提升了粮食产量服务，并且蓟州区部分城区和上游零散的农村居民点及不透水地面也促进区域内的产水量增加。

3.3.2 森林生态区

森林生态区位于于桥水库流域内淋河上游的西北部山区，由8～11 号子流域组成（图3），行政区划上主要位于河北省兴隆县境内，区域总面积600.04 km2，占流域面积的29.30%。该分区的土壤保持服务和固碳服务与其他分区相比均处于最高水平，氮和磷输出量低，说明水质净化服务较高。玫瑰图显示区域内的土壤保持服务及固碳服务与氮输出、磷输出、粮食产量服务存在强烈的权衡关系（图3），这可能与该分区大部分土地利用类型是林地有关，大面积的树木能够提高区域内的土壤保持和固碳能力，同时对水的需求量较大，因此使得产水量服务处于较低水平。此外，该分区内的氮、磷的输出源很少，并且林地本身就拥有较高的水质净化能力，因此森林生态区内水质净化服务极高。

3.3.3 粮食生产区

粮食生产区主要包括于桥水库流域的中部和南部平原地区，共10个子流域，总面积457.11 km2，占流域总面积的22.32%，行政区划上基本位于河北省遵化市境内。该分区内的各生态服务间关系与森林生态区内恰好相反，产水量、氮输出、粮食产量服务皆为4 个分区中的最高水平，土壤保持与固碳服务很低，粮食生产与固碳服务和土壤保持之间存在强烈的权衡关系。粮食生产区内耕地为主要土地利用类型，因此成为于桥水库流域的主要粮食供给区域。而许多农村居民点零散分布在区域内，形成了大面积不透水地面，增加了区域内的产水量，并且区域内的主要农产品作物为低耗水的玉米和板栗，使该区域的产水量服务非常高。大面积的耕地意味着该区域中会有大量的氮输出和磷输出，使该区域成为流域内主要的氮和磷输出区。

3.3.4 综合发展区

综合发展区主要包括流域北部和东部的1～5 号子流域，总面积为607.64 km2，占整个流域总面积的29.67%，为面积最大的分区，行政区划上属河北省遵化市和兴隆县。区域内的各生态系统服务间关系比较均衡，但服务水平均较低，其中，磷输出的强度大于氮输出，与水源保护区相反，可能是由于该区域内包含大量城镇，居民生活使用的洗涤剂含有大量磷。综合发展区中的土地利用类型较为全面，包括林地、草地、耕地和城镇用地等，遵化市城区位于综合发展区内，城区北部的燕山余脉山区以灌木林为主，对水量消耗较大，因此区域内整体产水量较低，同时山区的林地提供了一定的土壤保持和固碳服务。

4 讨论

本研究发现，产水量的高值区与城镇用地和农村居民用地基本重合，可能因为城镇与农村用地的人类活动程度高于其他土地利用类型，使得水分蒸散量较少。土壤保持量的高值区基本分布在流域西北部，此区域存在大面积林地，而林地对土壤流失有着拦截作用，因此土壤保持量较高；氮输出与磷输出的高值区与耕地的分布极为相似，这与农业生产中投入的大量化肥、农药有直接关系；粮食产量的高值区同样为耕地所在地，而流域中南部耕地粮食产量更高的原因可能是此区域降雨量较高；固碳量的高低直接受植被覆盖影响，而林地固碳作用的效果更好，因此高值区基本分布在林地。

在生态系统服务的权衡与协同关系中，粮食产量服务与水质净化服务间存在较强的权衡关系，粮食产量的提升往往依靠化肥、农药等农业投入品的增加，导致了氮、磷输出的增加，这与刘晓娜等[57]、Chen 等[58]的研究结果一致；水质净化服务与土壤保持服务间为协同关系，因为氮、磷的输出是以降雨冲刷和水土流失为基础，若土壤保持量增加，则氮、磷输出会随之减少；土壤保持服务和粮食供给服务为权衡关系，因为耕地中的翻耕、种植与收获过程都会导致土壤表层松动，这与汪仕美等[21]、闫晓露等[29]的研究结果一致。

基于生态系统服务簇的功能区划分，分析区域内各服务间的权衡与协同关系，明确其中主导生态系统服务，因地制宜地实施生态保护规划与管理措施，可为区域发展与生态建设提供决策支持。本研究针对于桥水库流域内的粮食产量、氮输出、磷输出、固碳量、产水量和土壤保持6 种生态系统服务，使用Kmeans 聚类分析法划分出4 类生态功能区，并提出发展管理建议。

水源保护区位于于桥水库库区，是水源保护关键的最后“一公里”。2000 年以来，为保护水库水质，政府对库区实施了一系列生态保护措施，如对水库周围进行环湖生态防护带建设，入库河段建设前置库湿地[59]。有研究发现，2010—2020年间流域内增加的林地基本都分布在于桥水库河流入库口处[60]，同时实施“南迁北管”措施，对水库南岸农村居民点进行搬迁，严控北岸开展对水质有害的经济生产活动，这些生态保护措施的实施使得水库水质得到了很大改善。于桥水库主要保护措施基本落实在水库入水口和水库南岸地区，针对水库北岸的生态保护措施略显不足，建议将今后的管理重心放在水库北岸的生态治理和区域内的农业面源污染治理，具体措施包括在17 号子流域内增设淋河入库林草缓冲带，在20号、24号和25 号子流域内进一步加强种植结构调整，发展有机农业、休闲观光农业；对北岸临近库区的农村居民点进行异地搬迁，并做好生态补偿工作。

森林生态区内的森林对整个流域的土壤保持、水源涵养、水质净化、气候调节发挥着关键作用，因此保护好流域内的森林资源十分重要。建议基于天津市和河北省的生态红线划分区[61]，加强流域内九龙山国家森林公园、上元古界国家自然保护区等自然公园的保护以及对龙门口水库和上关湖风景区的水质保护工作。严格控制淋河上游工矿企业，适当发展旅游、文化、康养产业，但需要注意避免清东陵及其周边地区的过度开发，建立生态保护和文化旅游增收平衡体系[62]。

粮食生产区是于桥水库流域中粮食生产、水量供给和氮、磷输出的主要区域，这一区域生态系统服务之间的主要矛盾是粮食生产和氮、磷输出二者间的强烈权衡。因此，本区域的治理应以保证粮食供给为基础，同时尽量减少氮、磷输出对下游水库水质的影响。具体措施可分为源头防控和过程防控。在源头防控方面，积极推广规模集约化农业[56]，提高水肥的施用效率，加强对当地农户科学用肥、用药的知识培训，尽量避免滥用化肥、农药，加大力度推广测土配方施肥、有机肥、缓释肥、生态沟渠、林草隔离带等面源综合控制措施；在过程防控方面，以治理水土流失为切入点，开展重点湖泊河流的岸边缓冲带建设，如在沙河支流汇流地水平口湿地开展植树造林，在沙河和黎河的干流附近设置林地缓冲带[63]。此外，加强该分区中部丘陵地区林草地的保护，防止其被转化为其他用途。

综合发展区是于桥水库流域的经济发展核心区，此区域内的各生态系统服务间关系较为均衡，因此，该区的管理方式应是在保证经济发展的同时，稳定各服务间的均衡关系。具体的措施包括控制城市化进程、避免城市化粗放式扩张。有研究表明，虽然城镇用地的扩张会增加产水量服务，但随之也会增加城市内涝的风险[64]。该区域内磷的输出量较大，主要是居民生活洗涤造成的，建议推广无磷洗剂，加强生活污水的收集和集中处理；此外，建议对沙河和黎河支流上的工矿企业进行管控，尽量避免工矿业污染物进入水体造成污染。

囿于数据的可获取性，本研究的粮食产量服务仅用玉米作物进行表征，通过土地利用变化直观反映粮食产量的变化，有一定的参考意义，但当下的模型无法将技术进步对粮食产量的影响纳入考虑范畴，今后的研究应考虑技术进步对粮食产量的影响，进一步提升模型模拟的准确性。此外，化肥减施增效、化肥零增长行动的实施，以及测土配方施肥等技术和管理手段提升所带来的氮磷输出的减少，也未在模型中进行精确的模拟，因此下一步研究中，开展中长期的实地监测，有利于更加精准模拟流域内的生态系统服务。

5 结论

本研究以于桥水库流域为例，采用InVEST 模型定量评估了流域2020 年6 项关键生态系统服务的空间分布格局，利用相关分析法揭示了各项服务之间的权衡协同关系，最后通过主成分分析和K-means聚类分析将于桥水库流域划分为4 个生态系统服务功能区，提出因地制宜的管理措施建议。主要结论如下：

（1）于桥水库流域生态系统服务的空间分布具有明显差异。产水量服务整体呈现北高南低的分布趋势，高值区主要分布在遵化市城区和流域北部的连片山区以及水库西北部局部地区；土壤保持服务表现出西北高、东南低的分布格局，高值区集中在西北部山区；水质净化服务中，氮输出与磷输出空间分布高度重合，高值区皆为流域中部及南部平原区；粮食产量服务与氮输出、磷输出的空间分布重叠，高值区分布在流域中部和南部平原区；固碳服务整体呈现西北高、东南低的分布格局，西北部是燕山余脉，森林分布广，固碳强度最高。

（2）各项生态系统服务之间的相关性分析表明：氮输出、磷输出与粮食生产之间互为协同关系；土壤保持服务和固碳服务互为协同关系；粮食产量服务与土壤保持服务间存在较强的权衡关系。生态系统服务之间的权衡协同关系和主成分分析显示，保护好流域水质的关键在于处理好粮食生产和水质净化的权衡关系，而加强森林的保护可以协同增加固碳和土壤保持服务。

（3）依据生态系统服务簇的特征，将于桥水库流域划分为水源保护区、森林生态区、粮食生产区和综合发展区4 个功能区。其中：水源保护区以产水服务和粮食产量为主，建议加强北岸的管控工作；森林生态区以土壤保持和固碳服务为主，建议加强森林的保护；粮食生产区以产水量、氮输出、磷输出及粮食产量服务为主，需要加强区域内农业面源的管控；综合发展区需要控制城镇化用地过度扩张，加强生活源和工业源污染物的管控。