生态系统服务功能是生态系统与生态过程所形成及维持的人类赖以生存的自然环境条件与效用^[1]。不仅包括向经济社会系统输入有用物质和能量(如食物、其他生产生活原料等)，还包括创造与维持地球生命支持系统所必需的环境条件(如洁净空气、水、自然景观等)，以及生态系统接受和转化来自经济社会系统的废弃物(污染物的吸收与降解、有害气体吸收等)^[2-3]。生态系统服务功能的定量评价，即生态系统服务价值的评估是区域生态建设、生态保护、生态功能区划、自然资产核算和生态补偿决策的前提和依据，成为生态学的研究热点^[3-5]。国内外学者就生态系统服务的框架体系^[5-7]、不同生态系统的生态服务价值化评估方法^{[4, 6, 8]}、区域生态服务的空间格局和动态变化^[9-10]、生态安全和生态环境管理^[11-13]、生态服务与区域福祉的耦合^[14-15]等方面展开研究。打造京津冀生态环境支撑区是河北省在京津冀协同发展战略中“三区一基地”定位的重要一环。为此，本研究以1980，2000和2015年Landsat TM/ETM遥感影像解译生成的土地利用数据为数据源，采用单位面积价值当量因子法，基于谢高地等改进的当量表^[3]，构建细化的土地利用二级分类当量表，计算并分析不同时段河北省各生态环境支撑区生态系统服务价值及变化。旨在通过生态系统服务价值有效发现生态系统服务特征、变化及存在的问题，有助于推进地区生态环境建设，不断提高各支撑区的生态系统服务功能，是建设经济强省和美丽河北的必然要求，也是实现京津冀地区可持续发展的迫切需要。

1 研究区概况

河北省地处华北平原，位于北纬36°03′—42°40′，东经113°27′—119°50′，居环渤海地区的中心地带，内环京、津，北与辽宁、内蒙古为邻，西与山西为邻，南与河南、山东接壤，具有独特的区位优势。河北省地貌类型复杂多样，是全国唯一具有坝上高原、山地丘陵、盆地平原、湖泊海洋的省份。河北省属北温带大陆性季风气候，春季干燥多风，夏季炎热多雨，秋季昼暖夜凉，冬季寒冷少雪。全省年日照时数为2 450~3 100 h，年平均气温0~13 ℃，多年平均降水量531.7 mm，主要集中在7—8月份。土壤类型多样，共有21个土类，55个亚类，164个土属，357个土种，主要以褐土、潮土、棕壤等为主。多数土壤质地为壤质，适宜农林牧业生产。河北省境内主要由滦河水系、海河水系、辽河水系、内陆水系4大水系组成。多年平均地表水资源量1.20×10¹⁰ m³，人均地表水量仅为全国平均水平的10%。受地形影响，河北省植被分布有明显的垂直差异，植物数量繁多。《河北省建设京津冀生态环境支撑区规划(2016—2020年)》将全省的生态环境保护建设分区划分为5个生态功能区，包括京津保中心区生态过渡带、坝上高原生态防护区、燕山—太行山生态涵养区、低平原生态修复区、海岸海域生态防护区(图 1)。各生态支撑区基本特征如表 1所示。

2 数据和方法 2.1 数据及来源

所采用的土地利用/覆被数据来源于中国科学院资源环境科学数据中心(http://www.resdc.cn)。土地利用数据以Landsat TM/ETM遥感影像为主要数据源，通过人工目视解译生成。土地利用类型包括农田、森林、草地、水域、居民地和未利用土地6个一级类型以及25个二级类型。该数据库是目前我国精度最高的土地利用遥感监测数据产品，已经在国家土地资源调查、水文、生态研究中发挥着重要作用^[16-17]。

2.2 基于当量因子计算生态系统服务价值 2.2.1 单位面积生态系统服务价值当量

采用千年生态系统评估(the millennium ecosystem assessment, MA)的方法，将生态系统服务分为供给服务、调节服务、支持服务和文化服务4大类，并进一步细分为食物生产、原料生产、水资源供给、气体调节、气候调节、净化环境、水文调节、土壤保持、维持养分循环、生物多样性和美学景观等11种服务功能(表 2)。单位面积生态系统服务功能价值的基础当量是指不同类型生态系统单位面积上各类服务功能年均价值当量(以下简称基础当量)。基础当量体现了不同生态系统及其各类生态系统服务功能在全国范围内的年均价值量，也是合理构建表征生态系统服务价值区域空间差异和时间动态变化的动态当量表的前提和基础。本研究以谢高地等^[3]对生态系统服务价值评价研究为基础，考虑到计算精度，本文采用生态系统二级分类，细化土地类型，使结果更加精确。由于缺乏对海洋生态系统服务功能及价值的系统研究数据，在本研究中尚未包括海洋生态系统。同时，在改进的生态系统服务功能中，充分考虑了水资源的供给服务，并提高了水域的水文调节服务的价值当量，这在缺水的京津冀地区是非常必要的。

2.2.2 单个标准单位面积生态系统服务价值当量因子的价值量

1个标准单位面积生态系统服务价值当量因子的价值量是指1 hm²全国平均产量的农田每年自然粮食产量的经济价值。直接量化不同类型生态系统对生态服务功能的潜在贡献能力。由于受人类活动影响，农田在自然条件下的粮食产量的经济价值很难确定。根据谢高地等的研究，单位面积农田食物生产生态服务价值D(即1个当量因子的价值)相当于当年研究区域平均粮食单产市场价值的1/7，即：

$ D = 1/7\mathit{S}/\mathit{A} $

(1)

式中：S——年农作物经济总价值；A——年农作物播种面积。

由于单位面积价值是引起生态系统服务价值变化的重要因子，以评估年份的价格为表征的价值量会出现由于价格因素引起的生态系统退化但生态功能服务价值增加的悖论。因此为了纵向对比结果的准确性，本文统一采用2015年河北省的单位面积的价值量。根据2016年《河北农村统计年鉴》的农作物播种面积和经济价值，计算可得河北省2015年单位面积生态系统服务价值当量因子的价值量为1 731.5元/hm²。

2.2.3 生态系统服务功能价值估算

针对每个栅格对应的土地单元，其生态服务价值估算：

$ \mathrm{ESV}_{i j}=Y_{i j} \cdot D \cdot A_{i j} $

(2)

式中：ESV_ij——计算栅格的生态系统总服务价值(元/km^-2·a^-1)；Y_ij——栅格对应土地利用类型的生态服务价值当量因子(表 2)；D——1个标准单位生态系统服务价值当量因子的价值量；A_ij——栅格对应土地面积与单位面积(1 hm²)的比值，本研究栅格分辨率为1 km，因此A_ij的值为100。

3 结果分析 3.1 河北省生态环境支撑区土地利用方式变化 3.1.1 河北省土地利用方式变化的总体特征

表 3为1980，2000和2015年河北省土地利用类型面积变化的总体情况，除城乡、工矿、居民等建设用地增加之外，其他生态用地都为呈现面积缩小的趋势。1980—2015年河北省城乡建设用地增加了33.2%，新增加面积中86.2%是由农田转换而来的，其次是草地占5.5%。1980—2015年河北省面积减少最多的土地利用类型为农田，转出3 610 km²，转进1 055 km²，减少面积2 555 km²。转出面积中90.8%转为城乡建设用地，转入类型中主要包括水域、草地、未利用地和森林。35 a间森林转出573 km²，78.2%和18.3%转为草地和农田，且主要在2000年之前。2000年后，森林转出面积为89 km²，而且60%以上转为工矿用地。同时，1980—2015年森林转入的面积为257 km²，山区植树造林、退耕还林、太行山绿化等工程的生态效应显现。659 km²水域面积消失，66.2%转为农田，围水造田现象严重。转出的300 km²的未利用地中1/2以上转为农田，27%转为建设用地，尤其是转为厂矿、大型工业区、油田、盐场、采石场等其他建设用地。697 km²草地转出，46.1%转为农田，30.0%转为城镇建设用地。虽然草地生态系统也有转进，但是分析发现转进的507 km²中，87.6%是由森林退化转为草地，这一过程主要发生在2000年之前。

3.1.2 不同生态环境支撑区土地利用方式变化的特点

由于不同生态环境支撑区的主要土地利用类型、社会经济发展速度等不同，导致其土地利用变化的规律又有自身的特点。平原区是城镇化速度较快的地区，1980—2015年土地利用方式整体变化是农田转为城乡建设用地，其他生态用地转为农田的逆生态化趋势。京津保生态区的农田用地转出758 km²，93.3%转为城乡建设用地，其中，46.0%为城镇用地，37.7%为农村居民点，9.1%为厂矿等建设用地。与此同时，263 km²的其他生态用地转为农田，82.1%的面积曾经是湖泊、滩地等水域生态系统，以白洋淀的减小最为明显。1980—2015年低平原生态修复区有1 542 km²(占平原区农田的4.3%)的农田转出，其中转为城乡建设用地的面积占97.9%，转出的农村居民点占1/2以上。而转进的只有89 km²，主要是植被覆盖度高的草地和河滩等生态用地。类似的，海岸海域生态防护区同样农田减少，城乡建设用地增加，但转为厂矿、大型工业区、油田、盐场等其他建设用地面积比例较高。且水库面积在过去的35 a内增加了114 km²，沙河水库、陡河水库等水域面积增加。与平原地区农田大面积转出的状况不同的是，1980—2015年坝上高原生态区却新增农田186 km²，其中60.2%是由草地转换的，29%是由荒漠开辟而来。燕山—太行山山地是河北省森林分布的主要区域，同时也是森林面积变化最大的区域，森林覆盖减少400 km²，其中太行山森林的减少最多，燕山太行山消失的森林里79.9%转为植被覆盖更低的草地，9.5%转为工矿用地，其余转为农田和居民用地。燕山—太行山森林减少主要是在2000年以前，2000年以后由于三北防护林，太行山工程等的实施，森林呈现增加的趋势。

3.2 河北省生态系统服务价值时空变化 3.2.1 河北省生态系统服务价值的总体特征

1980—2015年河北省生态系统服务价值的年平均值为3.56×10¹¹元。从生态系统类型分析，各生态系统服务的价值从大到小依次为：森林＞草地＞水域＞农田＞湿地＞荒漠，分别占总生态系统服务价值的34.2%，26.6%，19.1%，17.2%，2.9%，0.04%。从提供的生态服务功能来看，以调节服务为主，占总服务价值的62.8%，主要是水文调节和气候调节功能。其次是支持服务，包括土壤保持和保持生物多样性方面的价值。而食物生产和原料生产的价值量仅占生态系统总价值量的5.3%和3.7%，每年提供的美学景观价值占4.1%。空间上，生态系统服务价值的分布呈现西北高，东南低的总体格局。水域生态系统的生态服务价值最高，其次是燕山—太行山森林、高覆盖草地，东南部农田和坝上东部的覆盖草地，荒漠等地的生态服务价值较低(图 2)。

3.2.2 河北省生态系统服务价值的时空变化

以2015年现价计算，1980，2000，2015年河北省生态系统服务价值分别为3.60×10¹¹元、3.54×10¹¹元、3.53×10¹¹元。35 a间生态系统服务价值减少70亿元，且主要发生在2000年以前。

各生态系统和各种生态服务类型的服务价值均呈现下降的趋势(表 4)。各生态系统类型中，水域的生态系统服务价值减少的最多，其次是农田和湿地。其中，森林和水域生态系统服务价值的减少主要发生在1980—2000年的第一阶段，2000年后的第二阶段有所增加；农田和湿地的生态系统服务价值持续减少，第一阶段速率大于第二阶段；由于有大量森林在第一阶段转为草地，所以草地生态系统服务价值表现为先增加后减少。

各种生态服务类型的服务价值均为下降，下降主要发生在2000年前，调节服务减少的最多，其次是支持服务、供给服务和文化服务。空间上，生态系统服务价值减少的区域主要分布在城镇和农村居民点周围，河流、湿地周边，生态系统服务价值增加的区域主要分布在燕山—太行山草地、森林恢复区。

3.3 河北省生态环境支撑区生态系统服务价值及其变化 3.3.1 河北省生态支撑区生态系统服务价值的总体特征

1980—2015年京津保、坝上、燕山—太行山、低平原区和海岸海域地区5个生态环境支撑区的年平均生态系统服务总价值分别为1.74×10¹⁰，4.45×10¹⁰，2.36×10¹¹，3.75×10¹⁰，2.03×10¹⁰元。单位面积的生态服务价值分别为1.02×10¹⁰，2.31×10¹⁰，2.44×10¹⁰，8.76×10⁹，1.76×10¹⁰元/km²。单位面积生态系统服务价值最大的生态区为燕山—太行山、其次是坝上高原、海域海岸、京津保，最低的为低平原区，仅为燕山—太行山区的1/3。从生态系统服务功能的构成上看(图 3)，京津保和低平原区以农田提供的服务为主，占低平原区的67.0%，燕山—太行山以森林生态系统服务功能为最大，占总价值的43.6%，其次是草地和农田。坝上高原生态区各生态系统服务价值比例较为均匀，草地比例最高为34.4%。海域海岸生态功能区的生态系统服务价值以水域和湿地为主，二者占总服务价值的59.8%。从各生态支撑区的生态服务价值构成可以看出，低平原区和海域海岸区生态系统服务价值来源比较单一，坝上高原和燕山—太行山的类型较为丰富。从生态系统服务类型看，都以调节服务为主，尤其是水文调节的服务功能价值最大。其次是以土壤保持功能为主的支持服务，而食物、原料和水的供给价值在8%~24%之间，以美学景观为主的文化服务的价值量占总价值的5%以下。

3.3.2 河北省各生态支撑区生态系统服务价值的变化特征

1980年以来河北省生态系统服务价值减少的量中，京津保生态区所占比例最高占41.5%，约28.1亿元，其次是燕山—太行山生态区生态价值总量减少2.63×10⁹元所占比例为38.7%，低平原生态修复区占12.8%，海岸海域生态区占8.7%，坝上地区基本保持不变。纵向分析，京津保生态区1980—2015年生态系统服务价值降低15.0%，为生态系统服务功能降低最严重的生态区，远远高于海岸海域生态区的2.75%和低平原区的2.28%。燕山—太行山生态区尽管生态系统服务价值总量减少较大，但是由于生态区整体生态系统服务价值量高，所以减少比例仅为1.11%，坝上高原生态系统服务价值变化最小。图 4为1980—2015年不同生态支撑区生态系统服务价值量的变化。从图 4可以看出，京津保地区所有生态系统和所有生态系统服务价值都有所下降，尤其是水域和湿地面积的缩小让京津保地区生态调价服务功能价值大幅下降。相反，坝上高原的调节服务有所增加，但是由于坝上湖泊多属于季节性湖泊，不排除由于数据来源的季节性差异导致的湖泊面积提取带来的较大不确定性。燕山—太行山地除了草地生态系统服务价值有所上涨外，其他生态系统的服务价值均有不同程度下降，尤其是森林和水域，导致生态区生态调节服务大幅减小。海岸海域生态区水域生态系统稍有增加，但是远远不抵农田、湿地等的减少，导致整体生态服务功能退化。

4 讨论与结论 4.1 河北省生态系统服务价值变化及原因

目前关于河北省建设京津冀生态环境支撑区的生态系统服务价值估算研究主要集中在不同地区和不同生态系统的生态服务价值估算，如京津冀区域、环京津冀、张承地区、森林生态系统等^[18-20]。本研究的生态服务价值空间分布与变化趋势与相关研究结果较为一致，从空间上看，坝上高原、燕山—太行山丘陵区生态系统服务供需呈盈余状态，生态服务价值较高^[18]。张承地区的生态系统服务价值为1.895×10¹¹元/a，与本文的1.917×10¹¹元/a较为接近^[20]。从多年的变化情况看，无论是环京津地区还是整个京津冀地区，生态系统服务价值都呈现下降趋势^[18-19]，但由于研究时段的差异导致下降幅度有所不同。

生态系统服务价值在空间上受土地利用方式变化的直接影响。通过分析发现1980—2015年河北省土地利用方式变化的总体特征为：农田转为城乡建设用地，水域和草地等生态用地转为农田的逆生态化过程，2015年城镇和农业用地均为1980年的2.5倍。土地利用方式的逆生态化导致生态系统服务价值流由高到低，进而导致河北省整体生态系统服务价值降低。土地利用方式转变的深层原因主要有：(1)由于城镇化和工业化发展对的建设用地需求的增加，导致建设用地的扩张占用周边农田；(2)农村广大地区受人口和经济因素的驱动增加了对农田需求，导致湿地、草地等生态用地的公地悲剧发生。河北省生态系统服务价值降低的过程主要发生在2000年以前，之后有逐步恢复的趋势。

2000年，是河北省生态建设的一个重要时间节点。1999年河北省人民政府印发《河北省生态环境建设规划》后，2000年退耕还林工程在张家口市、承德市坝上地区开始试点，随后在河北全面展开，河北省生态环境建设进入一个新的阶段。2001年京津风沙源治理工程建设工程正式启动，河北省三北防护林体系建设、太行山绿化工程等生态环境建设范围越来越大。这些生态建设工程不断提升生态系统服务功能，基本扭转了生态系统服务功能恶化的趋势。但是，由于河北省生态环境基础较差，同时正处在经济转型发展，城镇化加快阶段，生态环境压力仍然较大。面临京津冀协同发展重大国家战略、规划建设河北雄安新区等重大机遇，河北省应以建设京津冀生态支撑区为指引，优化生产—生活—生态的“三生空间”布局。制定严格的土地利用制度，不断增加人们的生态环境服务的需求，提高城镇的建设、管理现代化，并进一步加强生态环境建设，提升生态系统服务功能，最终实现区域可持续发展。

4.2 当量因子法估算生态系统服务价值

当量因子法是目前计算生态系统服务功能价值的常用方法，有较为简单、易于操作和结果便于比较等优点，但短板也十分突出^[3]。首先，表现在不同时间段和空间尺度上，由于表示价值的价格基础不同，导致同一地区的生态系统服务价值在不同尺度上的计算结果存在较大差异。例如估算一个县的生态系统服务价值时，将其放在县、市、省、生态区、全国等不同尺度中计算的结果可能会存在很大差异，利用不同阶段的价格计算出的结果也不同。其次，目前的计算生态系统服务价值在空间上受土地利用方式变化的直接影响，基本上只考虑了土地利用类型的质变影响，对同一种土地利用类型中不同程度的量变差异考虑较少。植被覆盖程度、水体质量等的差异，虽然有少数研究^{[3, 10]}以NDVI，NPP，降水等作为考虑因子，但仍然不能细致地刻画生态服务的差异。下一步研究应深入考虑不同生态系统生境质量的影响，更加科学和精细地计算生态系统服务价值，以期为生态评价、生态补偿和生态保护等提供可靠参考。

4.3 结论

(1) 1980—2015年河北省土地利用类型中除城乡、工矿、居民等建设用地增加之外，其他生态用地都呈现面积缩小的趋势。各生态支撑区土地利用方式整体变化是农田转为城乡建设用地，其他生态用地转为农田的逆生态化趋势。

(2) 以2015年现价计算，1980—2015年河北省生态系统服务价值的年平均值为3.56×10¹¹元，各生态系统服务的价值从大到小依次为：森林＞草地＞水域＞农田＞湿地＞荒漠。1980—2015年河北省生态系统服务价值减少约7.00×10⁹元，且主要发生在2000年以前。

(3) 1980—2015年京津保、坝上、燕山—太行山、低平原区和海岸海域地区5个生态环境支撑区的生态系统服务总价值分别为1.74×10¹⁰，4.45×10¹⁰，2.36×10¹¹，3.75×10¹⁰，2.03×10¹⁰元。单位面积生态系统服务价值最大的生态区为燕山—太行山，其次是坝上高原、海域海岸、京津保，最低的为低平原区。

(4) 1980年以来河北省生态系统服务功能减少的价值中，京津保生态区占41.5%，其次是燕山—太行山为38.7%，低平原生态修复区占12.8%，海岸海域生态区占8.7%，坝上地区基本保持不变。35 a间京津保生态区生态系统服务价值减少最多，降低了15.0%。