洛阳市生态系统服务功能价值量化及空间分异研究.pdf

收稿日期：基金项目：国家自然科学基金面上项目（）；国家自然科学基金重点项目（）作者简介：管新建（），男，河南驻马店人，教授，博士，主要从事水资源配置与管理方面的研究工作 通信作者：孟钰（），女，河南南阳人，副教授，博士，主要从事生态水文学研究工作：【水环境与水生态】洛阳市生态系统服务功能价值量化及空间分异研究管新建，石悸锋，孟 钰，符育文（郑州大学 水利与交通学院，河南 郑州）摘 要：在社会经济高速发展和生态环境日益恶化的矛盾下，分析评估生态系统服务功能价值是推进经济、社会和生态协调发展的迫切需要。采用经济学方法对洛阳市生态系统即森林、草地、水域生态系统服务功能价值进行量化，并应用莫兰指数对生态系统服务功能价值进行空间相关性分析。结果表明：洛阳市生态系统服务功能具有非常大的经济价值，其中森林和水域生态系统服务功能价值较高，水域生态系统服务功能价值系数最高；西部及中部地区生态系统服务功能价值较高，主城区及附近区域的价值较低；各行政区之间森林生态系统服务功能价值具有明显的集聚特征，草地生态系统服务功能价值空间相关性较弱，水域生态系统则表现出较强的空间差异性；森林和草地生态系统服务功能价值为正相关关系。后续发展应重视水域生态保护，同时合理规划森林和草地资源分布，促进地区社会生态经济协调发展。关键词：生态系统；服务功能价值；空间自相关分析；洛阳市中图分类号：文献标志码：引用格式：管新建，石悸锋，孟钰，等洛阳市生态系统服务功能价值量化及空间分异研究人民黄河，（）：，（，）：，：；生态系统服务功能是指生态系统与生态过程所形成及所维持的人类赖以生存的自然环境条件与效用。城镇化的快速推进导致区域气候、生物多样性、生态循环过程等发生了巨大变化，并引发了水土流失、土壤污染、动植物数量急剧减少等生态环境问题。近年来许多学者对生态系统服务功能价值开展了广泛研究，包括全球、国家、流域及地区等尺度。国内大多采用当量因子法计算整个复合生态系统服务功能的价值，忽略了地区间的差异，以及采用经济学方法对单个类型生态系统的单项或多项生态系统服务功能价值进行量化评估，对区域整体生态系统服务功能价值的评估较少。开展区域整体生态系统服务功能价值研究能够有效提高人们保护环境的意识，并且可以为区域生态系统的协调管理提供理论依据，推进生态文明建设。空间自相关分析方法作为进行地理分析的常用方第 卷第 期 人 民 黄 河 ，年 月 ，法，主要用于城市发展、规划和管理影响的定性评价，评价水资源的空间分布特征以及区域生态系统脆弱性分布和聚集特征等。采用空间自相关分析方法研究生态系统服务功能价值的空间分布特征，分析其空间相关性和聚合特征，可以为管理者提供决策依据。笔者采用经济学方法对洛阳市整体及各个行政区的生态系统（包括森林、草地和水域生态系统）服务功能价值进行量化，并利用 软件通过单变量和双变量的空间自相关分析方法分析各个类型生态系统服务功能价值的空间聚集特征以及不同类型生态系统服务功能价值的空间相关性，以期为区域生态环境建设提供科学依据。研究区概况洛阳市位于河南省西部，地处华北地台与秦岭褶皱系交界处，横跨两个一级大地构造单元，由此决定了洛阳市地形、地貌复杂，类型多样。洛阳市面积 ，辖 县（新安县、宜阳县、伊川县、汝阳县、嵩县、栾川县、洛宁县）、区（涧西区、西工区、老城区、洛龙区、孟津区、偃师区、瀍河区），总人口 万人，气候属于暖温带大陆性季风气候。洛阳市具有较为丰富的森林资源，高等植物种类繁多，森林覆盖率达，主要集中在西部和西南部地区。洛阳市水资源丰富，北部有黄河流域，洛河、伊河、涧河由西南向东北在市区交汇，还有陆浑水库和故县水库等。本文研究时间为 年，因此行政区划采用 年的行政区划，其中偃师市为现在的偃师区，现在的孟津区由原孟津县和吉利区合并而成，同时考虑到涧西区、西工区、老城区、瀍河区面积较小且位置集中，均为洛阳市核心发展区域，因此并为主城区进行研究。研究方法 生态系统服务功能价值量化）生态系统服务功能分类。目前生态系统服务功能分类大多采用联合国千年生态系统评估中提出的分类方法，将其分为供给、调节、支撑及文化功能四大类。根据该分类方法，结合不同生态系统的差异性，辨识森林、草地、水域生态系统服务类别，见图（图中表示具备该类生态服务）。供给功能指生态系统直接为人类生活及生产提供各种产品及资源。森林能够为人类提供木材、食物及药材等，草地能够为人类提供肉、蛋、奶等食物，水域在提供鱼类产品的同时为人类生活、生产提供淡水资源等，因此选取产品提供、水资源供给等指标评价生态系统服务供给功能的价值。调节功能指人类从生态系统的调节作用发挥过程中获得的收益。空气质量调节即森林、草地生态系统通过吸收空气中的污染物等调节空气质量；气体调节指植物通过光合作用固定、释放；水分调节指森林和草地通过植被截流等减小径流量、蓄积水分；水域生态系统通过水库及河道等进行水量调蓄，同时通过水体自净作用改善水质；气候调节指水域生态系统通过水面蒸发吸收热量，起到区域调温增湿的作用；废弃物降解指草地等生态系统中微生物分解家畜粪便等释放营养物质。因此，选取净化空气、固碳释氧、涵养水源、水资源调蓄、水质净化、气候调节、废弃物降解等指标评价生态系统服务调节功能的价值。支撑功能指生态系统所形成的支撑发展的条件与效用。土壤形成即森林和草地生态系统通过植被覆盖减少土壤侵蚀、保持土壤肥力；初级生产指陆生及水生植物通过光合作用吸收 并制造有机物的能力，因与固碳释氧指标重复而舍弃该指标；营养循环即生态系统内部与系统之间循环、运输、积累营养物质；生活质量改善指水域生态系统为城市生态环境如绿地、园林、喷泉等基础建设提供水源，从而维持并改善城市生态环境和人们生活质量的作用。因此，选取土壤保持、营养物质循环、生活质量改善等指标评价生态系统服务支撑功能的价值。文化功能指生态系统的美学、文化和教育等方面的非物质效益，常选取科学研究、休闲娱乐功能等指标评价生态系统服务文化功能的价值。图 生态系统服务辨识）评价指标与量化方法。综合考虑生态系统服务类型、功能特征、数据的可获取性等因素，选取森林生态系统供给功能（林业产品提供）、调节功能（涵养水源、固碳释氧、净化空气）、支撑功能（土壤保持、营养物质循环）、文化功能（休闲娱乐）项功能指标，草人 民 黄 河 年第 期地生态系统供给功能（牧业产品提供）、调节功能（涵养水源、净化空气、固碳释氧、废弃物降解）、支撑功能（土壤保持、营养物质循环）、文化功能（休闲娱乐）项功能指标，水域生态系统供给功能（水产品提供、水资源供给）、调节功能（气候调节、水资源调蓄、水质净化、固碳释氧）、支撑功能（营养物质循环、生活质量改善）、文化功能（科学研究、休闲娱乐）项功能指标，采用经济学方法对区域生态系统服务功能价值进行量化分析，其指标计算方法见表。表 生态系统服务功能价值计算方法，功能指标适用生态系统计算公式符号说明供给产品供给森林、草地、水域 为生态系统提供各类产品的价值，为第 类产品单价，为第 类产品的产量调节涵养水源水资源调蓄净化空气固碳释氧废弃物降解气候调节水质净化森林（）（）为森林涵养水源的价值，为森林面积，为降水量，为森林蒸散量，为地表径流量，为水库建设单位库容投资，为净化水质成本草地 为草地涵养水源的价值，为草地面积，为产水系数，为径流系数水域（）为水资源调蓄价值，为研究区河川水资源和水库蓄水量总和，为地下水资源总量，为水的密度，为居民用水价格森林、草地 为评估生态系统净化空气功能的价值，为生态系统的面积，为生态系统单位面积吸收第 类污染物量，为第 类污染物的治理费用森林、草地、水域（）为评估固碳释氧功能价值，为生态系统净初级生产力，为固碳系数，为固碳价格，为释氧系数，为工业制氧价格草地 为废弃物降解价值，为牲畜粪便归还草地的比率，为草地牲畜的载畜量，为不同类型牲畜平均每年粪便中所含营养物质（、）的量，为氮肥 磷肥价格水域（）为水域生态系统气候调节的价值，为水面蒸发量，为水域面积，为水的汽化热，为空调效能比，为 水蒸发耗电量，为居民用电价格水域 为水质净化的价值，为排放废水量，为废水处理成本支撑土壤保持森林、草地（）()为土壤保持价值，、分别为土壤侵蚀模数和无林地土壤侵蚀模数，为挖掘和运输单位土方费用，为土壤容重，为土壤中、的含量，为、的价格营养物质循环森林、草地、水域（）为营养物质循环的价值，、分别为不同生态系统中参与循环的平均含氮、磷、钾量，、分别为磷酸二胺化肥价格、氯化钾化肥价格生活质量改善水域 为水域生态改善生活质量价值，为年生态用水量文化科学研究、休闲娱乐森林、草地、水域 为文化功能价值，为单位面积生态系统文化娱乐功能价值 空间自相关分析空间自相关分为全局空间自相关和局部空间自相关，全局空间自相关可以探究整个研究区域变量之间的相关性，局部空间自相关可以研究某一区域与邻近区域对某属性值的集聚程度，检验局部区域的空间关联性和空间异质性。当全局空间自相关分析结果呈现出较强的空间正相关时，可以进一步进行局部空间自相关分析。单变量空间自相关分析描述的是空间某一要素同一指标的空间关系，双变量空间自相关分析描述的是某一要素的一个属性与另一个属性的空间关系。一般采用 指数（莫兰指数）来描述空间相关性，单变量莫兰指数包括单变量全局莫兰指数和单变量局部莫兰指数，计算公式为 ，（）（）（，）（）（），（）（）（）式中：、分别为单变量全局莫兰指数和单变量局部莫兰指数，、分别为空间单元、变量的观测值，为 的均值，为 的协方差，为研究对象的数量，为研究对象之间的空间连接矩阵 中的元素。双变量莫兰指数包括双变量全局莫兰指数和双变量局部莫兰指数。双变量全局莫兰指数计算公式：人 民 黄 河 年第 期 ，（）（），（）（）（）式中：为双变量全局莫兰指数，、分别为空间单元、的变量、的观测值，、分别为观测值、的均值，、分别为观测值、的方差，为、之间的空间连接矩阵 中的元素。莫兰指数大于 表示空间正相关，值越大，空间相关性越明显；莫兰指数小于 表示空间负相关，值越小，空间差异性越大；莫兰指数等于 表示空间不相关，呈现随机性。当呈现较强的正相关性时可以进一步分析其集聚特征，通常采用 指数即局部莫兰指数来衡量空间单元属性和周围单元的相近和相异程度，通过绘制 聚类图可以更加清晰表明各研究对象的聚类特征，主要有高高、高低、低高、低低、不显著几种情况，根据其不同聚类特征对空间内各单元进行分析。结果分析 生态系统服务功能价值构成分析）森林生态系统。根据各生态系统服务功能价值计算方法，计算 年洛阳市森林生态系统服务功能价值，结果见表。年洛阳市森林生态系统服务功能价值为 亿元，各指标价值排序为固碳释氧涵养水源土壤保持净化空气产品供给营养物质循环休闲娱乐。供给功能、调节功能、支撑功能、文化功能价值占比分别为、。调节功能价值远高于其他功能的，其中固碳释氧和涵养水源的贡献较为突出，表明森林对区域气体调节以及水资源调节具有重要作用。文化功能价值表 年洛阳市森林生态系统服务功能价值亿元行政区产品供给涵养水源固碳释氧净化空气土壤保持营养物质循环休闲娱乐合计主城区孟津县洛龙区新安县栾川县 嵩县 汝阳县宜阳县洛宁县伊川县偃师市合计 占比较低，说明森林生态系统为人类提供休憩娱乐等有待进一步挖掘。洛阳市西南部栾川县和嵩县森林生态系统服务价值较高，而主城区、孟津县、洛龙区的价值较低，这与洛阳市森林分布情况相吻合，后续应加大主城区、孟津县、洛龙区植树造林的力度，合理配置森林资源，缓解生态压力，促进区域可持续发展。）草地生态系统。年洛阳市草地生态系统服务功能价值见表。年洛阳市草地生态系统服务功能价值为 亿元，各指标价值排序为固碳释氧土壤保持净化空气产品提供涵养水源营养物质循环休闲娱乐废弃物降解。供给功能、调节功能、支撑功能、文化功能价值占比分别为、。调节功能价值占比最大，其中净化空气和固碳释氧贡献较大，废弃物降解价值较低。洛阳市西部洛宁县草地生态系统服务功能价值最高，中部新安县、宜阳县和嵩县价值较高，主城区、洛龙区价值最低。应合理增加主城区草地覆盖面积，改善地区环境，保持并改进洛宁县等草地分布结构，增强其稳定性和韧性。表 年洛阳市草地生态系统服务功能价值亿元行政区产品供给涵养水源净化空气固碳释氧废弃物降解土壤保持营养物质循环休闲娱乐合计主城区孟津县洛龙区新安县栾川县嵩县汝阳县宜阳县洛宁县伊川县偃师市合计 ）水域生态系统。年洛阳市水域生态系统服务功能价值见表。年洛阳市水域生态系统服人 民 黄 河 年第 期务功能总价值为 亿元，各指标价值排序为气候调节休闲娱乐水资源调蓄水资源供给生活质量改善水产品供给水质净化科学研究营养物质循环固碳释氧。供给功能、调节功能、支撑功能、文化功能价值占比分别为、。调节功能价值占比较大，其中气候调节的价值远高于其他指标的，说明水生态系统在调温增湿上起着巨大作用，对城市小气候的形成具有重要影响。洛阳市北部新安县、孟津县和南部嵩县水域生态系统服务价值较高，而栾川县、伊川县、洛龙区和偃师市的价值稍低。价值差异与气候调节功能相关，而气候调节功能价值主要与水域范围及水量有关，新安县和孟津县有小浪底水库和西霞院水库，嵩县建有陆浑水库，水资源较为充沛，且水库有较好的调蓄作用，因此水域生态系统价值较高。表 年洛阳市水域生态系统服务功能价值亿元行政区水资源供给水产品供给水资源调蓄气候调节固碳释氧水质净化营养物质循环生活质量改善科学研究休闲娱乐合计主城区孟津县洛龙区新安县栾川县嵩县汝阳县宜阳县洛宁县伊川县偃师市合计 ）生态系统服务功能构成及分布。年洛阳市生态系统服务功能价值为 亿元，其分布情况见图。图 年洛阳市生态系统服务功能价值分布情况森林、草地、水域生态系统服务功能价值占比分别为、。洛阳市嵩县生态系统服务功能价值最高，洛宁县、栾川县、孟津县和新安县的也较高，而主城区、偃师市和伊川县生态系统服务价值较低。城镇化加速了中心城市生态系统服务功能的衰退，主城区经济发展迅速，但土地利用结构不合理，因此生态系统服务功能价值低。水域面积最小，但水域生态系统服务价值占比最大，且价值系数即单位面积水域生态系统服务功能价值远高于森林和草地的，因此应加强对水生态的管理和保护，同时应对绿化覆盖面积较低的主城区等进行合理规划，使生态与经济发展相协调。空间自相关分析通过 软件对不同生态系统类型的服务功能价值进行单变量空间自相关分析，得到莫兰指数散点图见图。图 年洛阳市生态系统服务功能价值莫兰指数散点图人 民 黄 河 年第 期 洛阳市森林生态系统服务功能价值全局莫兰指数为，为强空间正相关；草地生态系统服务功能价值全局莫兰指数为，草地覆盖面积较小，分布较为离散，为较弱的空间正相关；水域生态系统服务功能价值全局莫兰指数为，为负相关，说明水域生态系统服务功能价值在行政区之间离散特征较为明显，空间差异性较强。洛阳市生态系统服务功能总价值全局莫兰指数为，为正相关。进一步对森林生态系统和整体生态系统服务功能价值进行单变量局部空间自相关分析，绘制 聚类图，见图。从图 可以看出，洛宁县、栾川县、嵩县森林生态系统服务功能价值为高高聚集，新安县、孟津县、主城区、洛龙区、偃师市为低低聚集。高高聚集区应丰富林木物种以及扩大森林覆盖范围，使生态系统更加稳定；低低聚集区需要进行植树造林，增大森林覆盖面积，使生态系统更完整有效地为人类提供各项服务。栾川县整体生态系统服务功能总价值为高高聚集，洛龙区为低低聚集，其他区域聚集程度不显著，这显示出区域整体生态系统服务功能的差异性。对于高高聚集区的后续发展需要在目前的发展策略基础上保持稳中求进，以巩固已有的服务功能价值为管理目标。低低聚集区需要对整体生态系统管理的图 年洛阳市生态系统服务功能价值 聚类图方案进行改进，寻找与区域更为匹配的发展战略，提升区域生态系统服务功能价值。同时需要改善草地结构，增强区域间水资源的调控和管理。通过 软件，采用双变量自相关分析方法对不同类型生态系统之间服务功能价值的相关性进行分析，计算得到：森林生态系统服务功能价值与草地、水域生态系统服务功能价值之间的莫兰指数分别为、，草地生态系统服务功能价值与森林、水域生态系统服务功能价值之间的莫兰指数分别为、，水域生态系统服务功能价值与森林、草地生态系统服务功能价值之间的莫兰指数分别为、。森林生态系统与草地生态系统之间服务功能价值相互影响，为较强的空间正相关关系；水域生态系统服务功能价值与森林生态系统服务功能价值为一定程度空间正相关关系，与草地生态系统服务功能价值的相关性较弱，同时水域生态系统服务功能价值的空间差异性较大，水域生态系统对森林和草地生态系统服务功能价值的影响较小。目前区域发展越来越注重生态环境保护，提高草地覆盖率相对于提高森林覆盖率更加容易，但森林的生态效益更强，因此需要提高城市森林覆盖率。森林和草地生态系统呈现较强的空间相关性，说明在同一区域草地和森林生态系统在提供生态系统服务时具有互相促进作用，与完全的草地分布相比，区域内一定林地分布能够使草地生态系统为人类提供更多的服务价值，因此可以增加森林和草地的复合生态系统来提高城市绿化质量，促进生态系统协调健康发展。针对水域生态系统的空间差异性，要加强水生态的维护，注重水资源的合理配置和调控，缓解水资源地区分布不平衡问题。结论对洛阳市 年生态系统服务功能价值进行量化，通过空间自相关分析发现，生态系统能够为人类带来巨大的经济效益，其中水域生态系统价值系数最高，但地区差异性较大，因此需要重视水资源保护以及水资源地区间的调控；森林生态系统比草地生态系统的经济效益更显著，且两者在空间上具有一定集聚特征，可以根据呈现出的聚集特征调整各行政区发展策略。同时森林和草地生态系统服务功能价值具有较强的空间正相关性，能够为城市绿化规划中森林、草地的分布设计提供一些新思考。参考文献：欧阳志云，赵同谦，王效科，等水生态服务功能分析及其人 民 黄 河 年第 期间接价值评价生态学报，（）：，：，：，：，：谢高地，鲁春霞，冷允法，等青藏高原生态资产的价值评估自然资源学报，（）：谢高地，张彩霞，张雷明，等基于单位面积价值当量因子的生态系统服务价值化方法改进自然资源学报，（）：欧阳晓，贺清云，朱翔多情景下模拟城市群土地利用变化对生态系统服务价值的影响：以长株潭城市群为例经济地理，（）：刘胜涛，高鹏，刘潘伟，等泰山森林生态系统服务功能及其价值评估生态学报，（）：赵苗苗，赵海凤，李仁强，等青海省 年草地生态系统服务功能价值评估自然资源学报，（）：梁鸿，潘晓峰，余欣繁，等深圳市水生态系统服务功能价值评估自然资源学报，（）：饶恩明，肖燚，欧阳志云，等海南岛生态系统土壤保持功能空间特征及影响因素生态学报，（）：刘景红，郑晓，樊俊美，等基于 模型的浑河中上游水源涵养服务价值评估应用生态学报，（）：田宁宁，张建军，李玉婷，等晋西黄土区退耕还林地涵养水源和保育土壤功能评价水土保持学报，（）：，：吴泽宁，黄硕俏，狄丹阳，等黄河流域农业系统水资源价值及其空间分布研究灌溉排水学报，（）：娄佩卿，付波霖，刘海新，等锡林郭勒盟草地生态系统服务功能价值动态估算生态学报，：（）：昝欣，张玉玲，贾晓宇，等永定河上游流域水生态系统服务价值评估自然资源学报，（）：侯孟阳，姚顺波，邓元杰，等格网尺度下延安市生态服务价值时空演变格局与分异特征：基于退耕还林工程的实施背景自然资源学报，（）：李小璨，阿荣，佟宝全呼和浩特市生活性服务业空间布局特征及评价干旱区地理，（）：，：施馨雨，赵筱青，普军伟，等基于斑块尺度的云南省景观生态安全时空演变及归因生态学报，（）：，：【责任编辑 吕艳梅】（上接第 页）无资料地区传统的支沟防洪分析计算过于单一，仅依靠地区水文手册及暴雨洪水图集计算，不能体现现代城区排洪渠的角色定位，本文对于受地势影响需承接管网排水、超标准城区坡面流及天然塬面汇流的支沟防洪分析计算提供了解决方法。）初步探索了水利行业及市政行业的重现期标准衔接问题，城区超标准坡面流概化计算法仅对于市政，与，相近时是可行的；对于转换计算后设计暴雨量相差较大的情况，还有待进一步研究。参考文献：，：，：张建云城市化与城市水文学面临的问题水利水运工程学报，（）：潘成忠，上官周平降雨和坡度对坡面流水动力学参数的影响应用基础与工程科学学报，（）：丛翔宇，倪广恒，惠士博，等基于 的北京市典型城区暴雨洪水模拟分析水利水电技术，（）：臧文斌，赵雪，李敏，等城市洪涝模拟技术研究进展及发展趋势中国防汛抗旱，（）：【责任编辑 许立新】人 民 黄 河 年第 期