生态系统服务是社会经济发展的重要物质基础和支持^[1]。由于城市扩张、自然资源开采和不合理区域发展导致生态系统服务普遍下降, 保护和提升生态系统服务已经成为自然空间管制的重要目标^[2]。生态系统服务制图是根据决策需求将生态系统服务空间分布及相关要素进行可视化展示的过程。其制图方法和内容丰富, 如生态系统服务供给、需求及服务流制图等^[3-4], 但制图最终目的是将生态系统服务纳入生态安全规划和自然资产管理^[5]。

生态系统服务空间规划和管理需兼顾自然-社会经济系统和多方面利益主体进行统筹规划和决策。这需要制图反映生态系统服务依存的自然条件、人为干扰因素、及其受益区的社会经济状况等详细的背景信息。绘制生态系统服务分布图直接作为辅助决策的工具, 是综合制图学、地理信息学和生态学理论方法系统性的工作, 是对传统制图方法和技巧的挑战。目前生态系统服务制图普遍存在缺乏与生态系统服务相关的高精度的地理空间信息^[6], 及忽视制图基本规范而达不到制图较高的质量要求^[7]等问题。因此, 生态系统服务分布图决策应用性受到较大的局限。大比例尺生态系统服务分布图绘制方法还缺少统一规范和要求, 极大限制了生态系统服务评价结果直接服务于地方生态环境保护规划。

因此, 为了提高生态系统服务分布图的质量、规范性和决策应用性, 本研究系统的结合制图学理论方法及国家、行业规范, 提出生态系统服务分布图的编制方法。以北京市延庆区土壤保持服务为研究案例, 评价了生态系统土壤保持服务分布图编制方法的可操作性和应用价值。探寻生态系统服务分布图的编制方法以求从图面识别生态系统服务重要区域是本研究拟解决的科学问题。从而使编制的生态系统服务分布图可直接作为生态保护规划的基础图件。

1 生态系统服务制图方法

为提高制图的规范性、科学性和通读性, 本研究参考制图学理论和国家行业相关制图规范, 提出了生态系统服务制图的目的、原则、流程和基本方法。

1.1 制图目的和原则 1.1.1 制图目的

(1) 反映生态系统服务的空间分布规律及与相关自然、人文基础要素之间的关系。

(2) 以明确的应用性目的为导向, 方便使用者直接查询生态系统服务高低值分布区域和关注地区的生态系统服务分布状况。

(3) 生态系统服务图件直接作为生态保护规划的基础图件, 为保护优先区划定等生态保护规划决策服务。

1.1.2 制图原则

(1) 制图科学性。制图遵循数据资料的现势性, 基本符号与国家、行业规范的统一性, 空间位置精确性, 表达艺术性和制图信息量合理性等要求。

(2) 尺度匹配性。由于生态系统服务与一定的自然和社会过程时空尺度相匹配^[8], 生态系统服务制图尺度要与专题要素、地理基础要素以及决策的空间尺度相一致。

(3) 地理要素相关性。在有限的地图空间内, 展示生态系统服务与地理要素的关系, 需选择与生态系统服务、决策规划目标直接或间接关系密切的自然、人文要素。

(4) 重点突出性。突出反映生态系统服务空间差异特征及重要过程因素, 服务于生态环境保护决策。

1.2 制图流程

本研究采用专题制图的方法编制土壤保持服务分布图, 主要的制图过程为(图 1)：(1)确定制图的数学基础, 以精确定位生态系统服务及其相关要素的地理位置。(2)选择自然、人文基础要素。并用地图学中的综合方法, 基于比例尺、制图目标和图面信息量控制等, 对要素进行分类、简化。采用适当的符号和注记表示, 展现各要素大小、密度、级别等主要信息。(3)通过生态系统服务数据分级和色彩组合绘制生态系统服务专题图层, 展现生态系统服务的数量空间分布特征。(4)综合基础底图和土壤保持服务专题图层, 进行图面配置。

1.3 基础底图设计 1.3.1 数学基础设置

在县域尺度的生态保护规划和落实需要协调村落管制单元, 在合适的图幅范围内展现村落居民点及名称信息, 制图适宜采取1 ∶ 250000比例尺, 其他数学基础设置参考1 ∶ 250000比例尺国家、行业相关规范、图式^[9-10]。本图采用Gauss Kruger 6分带投影, 1985年国家高程基准、2000年国家大地坐标系。设置4km间距公里网, 及大格间距10′, 小格间距1′的经纬网。

1.3.2 地理要素选择与综合

(1) 地理基础要素选择。

自然基础要素选择首先考虑县域尺度影响生态系统服务分布格局的生态因子。植被是生态系统服务产生和传递的基础^[11], 中小尺度影响植被生长格局的主要生态因子有地形、降水和气温^[12], 进而影响着生态系统服务。地形是其他因子的控制因素, 因此, 选为本图重要的基础自然要素。由于气温和降水变化较大, 考虑降水和地表径流有密切关系^[13], 选择较稳定的水系为另一自然要素。建设用地的扩张是生态系统服务变化的主要人为因素之一^[14], 居民地(包括居民点和居民面)和交通要素为基础建设用地, 作为制图的基础人文要素。进一步结合生态系统服务相关性、规划决策的可参考性及制图可行性三个方面阐述所选要素的代表性, 如表 1。此外, 地名是实现制图空间定位、查询功能必不可少的要素, 也可显示空间管制单元信息。因此, 也选取详细到村落的居民点名和自然地名要素为基础要素。

(2) 地理要素分类和筛选。由于地理基础要素信息量大, 需对其分类和筛选, 特别针对公路、居民点、自然地名和水系线。主要基础要素的分类依据和分类结果如表 2所示。

(3) 地理要素简化。地形要素主要采用适于比例尺精度要求的100m高程数据^[22]通过高程点、等高线、山体阴影提取三种方式, 增加地形要素的信息量和表现力。行政注记是体现管制单元的重要注记, 首先显示镇、乡等重要级别的注记。化简除去乡、镇政府驻地周边较密集的村落注记, 避免注记压盖。

(4) 符号注记表示。符号注记基本样式参考相应比例尺的国家、行业规范图式^{[9-10, 23]}, 以保证符号注记的通读性。对以上规范中没有涉及或与ArcMap平台不匹配的符号注记参数, 通过制图实践进行调整补充。变量设置遵循符号注记设置原则^[24], 由于注记较多, 通过大小、位置、权重等多变量设置避免注记符号之间的压盖。

1.4 生态系统专题图层

生态系统服务分布专题图层主要通过等值线方法表现, 包括分级方法选择和分级色彩模式设计。

1.4.1 生态系统服务数据分级方法

分级方法包括分级界限和分级数的选择。郭庆胜等研究提出为符合原始数据空间自相关分布的特征, 分级数应设置在5—8级为宜^[25]。增加级数会减少专题数据信息量的损失, 但过大也会降低读图效率^[26], 增加颜色模拟的难度。级数选择也与制图尺度有关^[27]。在知网以“(生态功能OR生态系统服务)AND县域”关键词进行检索, 统计县、市、区尺度的生态功能或服务评价使用的分级数, 共搜到的19篇文献中使用的分级数为3—5级^[28-46]。因此, 综合空间自相关合理性, 县市尺度生态系统服务评价常用分级数的范围, 及降低制图的难度和提高读图效率, 将分级数设为5级。表示为低级、中低级、中级、中高级和高级。分级界限设置方法参考《生态红线划定指南》^[47]生态功能重要性评估使用的累积比例法。

1.4.2 生态系统服务分级色彩模式

相关生态系统服务制图研究中少有学者对分级色彩模式进行探讨。而本研究方法编制生态系统服务分布图, 将绘制出大比例尺生态系统服务规律, 又能展示级别间破碎化小斑块的差异。同时色彩模式考虑与基础符号色彩的协调性, 兼顾分布图表现的艺术性。由于多方面综合性要求, 专题色彩模式的选择是本研究制图的难点。

在制图过程中, 专题图颜色组合设置遵循用色逻辑^[48-50], 选择亮度变化由浅到深的连续型逻辑颜色组合^[51]表示生态系统服务低值到高值的次序级别。综合考虑颜色的寓意和心理、“视觉对比”、“视觉平衡”等视觉效应。使用自然色调、与基础符号暖色产生对比的中冷色相。由于色相逐渐过渡的颜色组合各级别间辨识度差, 设置中间色相为灰色^[50], 以增加级别间的辨识度。

2 延庆区土壤保持服务制图 2.1 延庆区介绍

延庆区隶属北京市, 东经115°44′—116°34′, 北纬40°16′—40°47′, 属于温带与中温带, 半干旱和半湿润的过渡带, 山区面积占72.8%, 平原占26.2%^[52]。延庆是北京乃至整个京津冀地区的重要生态屏障和生态涵养区, 是保障生态安全的关键区域。在延庆各种生态系统服务经济价值核算中, 调节服务占总价值的91.4%。其中土壤保持是延庆生态系统十分重要的调节服务^[53]。由于延庆地貌复杂破碎, 山区植被郁闭度低, 雨季集中于夏季, 容易造成水土流失。近年来延庆加快设施建设, 大力发展旅游业。城镇、村庄、交通用地增加迅速, 天然草地和林地破坏, 水土流失导致地质灾害, 水库泥沙淤积问题增加^[54]。延庆以建设国际一流生态文明示范区为目标, 在补齐设施短板的同时, 保护山区的生态系统服务功能面临较大的挑战。通过本研究对延庆土壤保持服务分布图的绘制, 可为延庆的生态保护规划管理和设施建设提供参考和借鉴。同时可为其他县区生态保护重要区域识别提供方法和技术参考。

2.2 数据来源

专题数据使用优化的土壤流失方程测算的2015年90m精度的土壤保持服务强度数据(单位：t hm^-2a^-1)代表土壤保持服务的功能量^[14], 数据来源于中国科学院生态环境研究中心。底图数据使用了地理空间数据云网站(http://www.gscloud.cn/)ASTER GDEM V2 2009 30m精度的高程数据及地理专业知识服务系统网站(http://kmap.ckcest.cn/)2015年公众版1 ∶ 25万矢量数据。土壤侵蚀强度数据来源于考虑提拔中国科学院生态环境研究中心(http://www.ecosystem.csdb.cn), 强度分级参考《土壤侵蚀分类分级标准(SL 190-2007)》。

2.3 制图方法

基于ArcMap平台绘制延庆土壤保持服务分布图, 主要的制图方法、工具和参数如表 3所示。

2.4 土壤保持服务分布图

从土壤保持服务分布图看(图 2), 各地理基础要素分布信息清晰易读。其中, 等高线、山体阴影和高程点综合展现了宏观立体的地形特征和微观的地貌细节。从图面(图 2)容易读取河流、湖库等水系分布及名称信息。居民面和居民点要素多分布于800m高程以下的盆谷区。居民面多聚集分布于延庆盆地, 也有小面积分布于东北部山区的谷地如千家店镇、四海镇、珍珠泉乡等, 居民点分布较为分散。级别较高的铁路、国道、省道多分布于人口集聚区或较开阔的山谷区。县道数量较多, 在人口聚集区及地势较陡的山区都有分布。街道和乡村道路多分布于城区及山区谷地。图中的行政注记层次分明, 村落等居民点、自然地点名及道路注记都可清晰辨识。

图面(图 2)显示的专题要素看, 土壤保持服务各级别颜色组合能够清晰详尽的显示出土壤保持服务的分布特征及其与底图要素的分布关系。总体上, 土壤保持服务整体呈现低值区面积大, 聚集分布在延庆区西南部及中部盆地, 高值区面积小且破碎化分布在盆地周边山区的特征。最高级别土壤保持服务相对累积量为24.9%(表 4), 所占面积为82.5km², 仅占全区面积6.9%, 远小于其他级别。最低级别相对累积量为10.0%, 总面积为618.2 km², 约占全区面积51.8%, 占比远高于其他级别。

从图 2反映的土壤保持服务分布与地形、水系、公路和居民地等要素关系看, 土壤保持服务低值区主要分布在等高线较为稀疏的盆谷地, 及狭窄的沟谷区, 土壤保持服务较高值区主要分布在等高线较密集陡峭的山区, 但山区等高线多蜿蜒曲折, 地形复杂, 土壤保持服务高值区呈破碎化分布。居民面和居民点周边土壤保持服务多为低值, 尤其在居民地密集区土壤保持服务呈现低值区大面积聚集分布。河流沿线土壤保持服务多为低值区。省道、高速、国道和县道等公路周边土壤保持服务多为低值区, 山区村道等公路两侧土壤保持服务低值区较窄或与周边土壤保持服务没有明显差别。

综合分析各要素与土壤保持服务的关系, 总体上呈现出土壤保持服务中高值区分布与交通要素、居民点和居民面等基础要素空间分异的特征。

3 讨论 3.1 地图要素综合和呈现

地图要素制图采用500m PEAK线平滑法绘制等高线, 与其他等高线制图方法相比较, 地貌宏观特征保留程度较好且较为平滑, 高程标准差^[57]为17.9m, 符合高程精度的质量要求。从土壤保持服务分布图分析(图 2), 等高线较符合实际的展现了土壤保持服务发生发展的地形因素。同时, 地形要素结合水系要素可直观展示县区流域内由于水土流失导致的泥沙运移、堆积的可能路径, 有助于识别土壤保持服务下降潜在的风险。目前也有学者绘制精细到村落注记的生态系统服务分布图^[45], 绘制高程阴影展现地貌特征^{[3, 34]}。本图较其他学者更精细绘制公路、等高线等诸多要素符号注记, 提供了更精准的定位和土壤保持服务发展的环境信息。

本图专题图层设计采用的累积比例分级方法相较于聚类法制图^[58]相比, 可识别出破碎化的热点和被干扰的斑块, 直观展现小尺度上土壤保持服务的干扰因素。颜色组合是呈现土壤保持服务分布特征的关键符号, 本图的颜色组合辨识度高、次序感强、和谐美观, 且避免了与底图符号色彩的冲突, 有助于提高土壤保持服务分布信息读取的效率和准确性。

制图通过以上易于实现的多种要素的综合方法和制图技巧, 有效的提高了图面信息量和精细程度, 展现了土壤保持服务与基础自然和人文要素的关系。坡度和坡形对土壤保持的分布格局产生着关键的作用。山区复杂的地貌如沟谷是土壤保持服务重要的自然干扰因素, 加上山区水系丰富, 河流侵蚀加剧了对土壤保持服务的干扰。缓坡区土壤保持服务普遍较低, 有村落、景区散落分布, 村道等低级别公路交错分布, 虽然规模相对较小, 这些基础建设及周边其他人为因素对土壤保持服务产生了较大的干扰。盆谷地区基础建设规模较大, 及其他未绘制的如农业要素是土壤保持服务普遍较低的重要人为因素。

3.2 土壤保持服务重要区域的识别

土壤保持服务分布图(图 2)显示出土壤保持服务自然、人为干扰要素分布情况, 及其地形、水系、交通和居民地要素综合展现了供给区-流-受益区的关联信息。进一步通过这些信息, 结合土壤侵蚀强度分布图(图 3), 延庆地质灾害研究等资料, 可通过图面识别土壤保持服务重要的分布区域。

从图面(图 2)分析, 张山营镇周边山区是其水系下游重要经济建设区与水利设施的土壤保持服务重要供给区, 如延庆城区、龙庆峡水库、香村营水库及北京市的饮用水源地之一官厅水库等。该山区水系丰富, 多处为陡坡区, 土壤受到的侵蚀力强。尤其是山区与盆地交界处的陡坡区, 多为土壤侵蚀强度、中强度地区, 如白石头坑附近为土壤侵蚀剧烈区, 应梦寺山附近为极强区。相关研究也表明泥石流灾害频发区^[59]临近其山下G110国道西南段。同时, 龙庆峡4A级旅游景区、玉渡山风景区(大北梁西部)等公路设施对山区的土壤保持服务产生了较大的干扰^[59], 周边多为土壤保持服务低值区。综合来看, 张山营镇北部山区为土壤保持服务重要区域, 需重点协调旅游开发建设和山区土壤保持服务的保护。

此外, 四海镇南部山区沟壑险峻, 水系呈放射状汇集于菜食河, 里程较长, 沿线分布着较多的居民建设区。研究也表明该山区是泥石流灾害的频发区^[59], 如省道213四海镇南段海字口、省道213永安堡村附近。山区中凤凰驼、长城等自然景区受游人登山远足等活动干扰频繁, 土壤保持服务多处为中低值区。从图面综合分析, 四海镇南部山区也为土壤保持服务重要区域。

八达岭长城风景区为延庆长城文化带的核心区域, 相较于其他山区是受旅游开发建设活动影响最剧烈的地区^[60], 交通要道和观光线路纵横交错, 土壤保持服务被分割为破碎化分布, 研究也表明石峡-八达岭一带是地质灾害高发区^[59], 多处发生崩塌, 保护该地区的土壤保持服务对防止滑坡等地质灾害, 保证游客安全有重要意义。

此外, 其他的交通要道如省级公路S309千家店西部山区后坑、干沟等, 及千家店镇东端上奶山村等, 省道S212大庄科乡北段和南段这些交通要道临地势险峻的山区分布, 研究也表明这些地区多处易发生崩塌、泥石流等地质灾害^[59]损害公路设施, 这些要道周边山区也为土壤保持服务重要区域。

3.3 决策服务功能评价

本研究编制的土壤保持服务分布图决策应用功能主要体现在以下三个方面：

(1) 通过详实的专题要素和基础要素信息及联系, 可实现从图面识别生态系统服务重要区域。这很大程度上提高了分布图辅助生态保护规划决策的效率。制图显示出延庆山区土壤保持服务功能的重要性和脆弱性, 延庆土壤保持服务的保护应该注重山区保护, 应对山区缓坡区村道等公路、景区建设及其他人类活动的干扰采取管控措施。

(2) 制图实现了较精准的定位和查询功能。可识别土壤保持服务重要区的具体位置, 查询关注要素周边的土壤保持服务分布情况。既为规划决策及实施提供便利, 也有助于与其他相关生态评价规划成果关注点的比较和信息整合, 对更合理的引导政府将资金投入到更关键的区域有重要的意义^[61]。

(3) 图件较详细的呈现了规划决策需要协调的多方面的利益主体信息, 如乡镇政府、村落、景区或者工厂企业等。方便考虑需协调的多方面的利益相关者, 以便保护规划决策开展和落实。

4 结论

结合制图学理论方法及国家行业相关标准和规范, 提出了生态系统服务分布图编制的目标、原则、流程和制图方法。通过延庆区土壤保持服务分布图制图方法实践, 显示出该制图方法能够准确详细的展示生态系统服务与相关自然、人文要素的分布特征, 并通过土壤保持服务与基础地理要素的联系, 有助于识别土壤保持服务重要区, 为区域生态环境保护规划提供重要决策参考。

本研究提出的生态系统服务分布图编制方法可行性强。基于开源数据库, 使用ArcMap工具一体化制图。图式与国家、行业相关规范图式基本统一, 制图方法更简便易行, 可为相关专业人员提供制图参考, 促进生态系统服务制图的规范化和自动化。

致谢: 饶恩明老师提供2015年生态系统土壤保持服务强度数据, 特此致谢。